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The relative importance attached to the field strengths is perhaps one of 


‘the main points of contrast between the classical and the quantum theory of 


the electromagnetic field. Whereas to MAXWELL the field strengths were 
probably of prime importance in themselves, in modern physics they have 
come to play an auxiliary role, second to that of quantities such as the cross 
sections or the shift of energy levels. 

This state of affairs, however, does not necessarily imply that the field 
strengths are devoid of any physical significance outside the range of validity 
of the classical theory. Before such a conclusion could be reached, it should 
have been shown, by the discussion of appropriate experiments, that the field 
strengths are non-measurable quantities; namely, that their measurement is 
subject to uncertainties far stronger than those accounted for by the theory. 
This point seemed to have been confirmed in a paper by LANDAU and PEIERLS (!), 
which appeared as early as 1931. In the ideal experiment described by these 
authors, the measurements of a field strength at a certain point of space- 
time is reduced to that of the momentum of a point charge used as a test body. 


(*) This review is based on a series of informal seminars delivered by the author 
to a small group of theoreticians in the N.R.C.. The author’s thanks are due to 
Prof. L. RoseNFELD of Manchester University, for introducing him to the subject. 

(**) National Research Laboratories Postdoctorate Fellow. Now at Palmer Phy- 
sical Laboratory, Princeton University, Princeton, NIJESZU:-S7A% j 

(1) L. Lanpau and R. PEIERLS: Zetts. f. Phys., 69, 56 (1931). 
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This last, being subjected to an acceleration during the measurement, emits 
a certain amount of radiation. The impossibility of controlling the energy 
and momentum of such a radiation was proved to be one (and the largest) 
source of uncertainty in the measurement. 

However, BoHR and ROSENFELD (*) (1933) pointed out that the result 
reached by LANDAU and PEIERLS was due to the atomic nature of the test 
particle, instead of being the consequence of the properties of the electromagnetic 
field. No difficulty arises if test bodies of finite extension are employed, thus 
measuring average values of field strengths over finite space-time regions, 
instead of field strengths at a certain point. Moreover, the latter case can be 
. treated as the limit of the former when the dimensions of the region, over 
which the average is performed, tend to zero. 

The question of the measurability of charge-current quantities was also 
considered by Bonr and ROSENFELD (*) (1950), who indicated the way in 
which average charge-current components can be unambiguously measured. 

In the following we shall expound the above points of view. We shall 
also indicate how methods recently developed in quantum field theory (Tomo- 
NAGA, SCHWINGER, FEYNMAN, Dyson et al.) can be advantageously used to 
improve the discussion of the more formal part of the problem of measurability. 


1. — Measurements of field strengths. 


1:1. — We shall first summarize LANDAU and PEIERLS’s discussion of the 
measurement of a component of the electric field (1). 

Let e be the charge of the particles used as a test body, # the component 
of the electric field in a certain direction, and p(t) the component, in the same 
direction, of the momentum of the test body at the time t; tbis last quantity 
is assumed to be known. The test body is accelerated by the electric field 
during a time interval 7, at the end of which the momentum has the value 
p(t + T) according to the equation 


(1) ot + T)— p(t) =eET. 


The measurement of p(t + 7) is, therefore, equivalent to that of E, while 
of course any uncertainty Ap in the value obtained for the former quantity 


(2) N. Bonr and L. Rosenrerp: Kgl. Danske Vid. Sels., Math.-fys. Medd., 12, 
No. 8 (1933). 

(Cf. also W. HeIrLER: Quantum Theory of Radiation (Oxford, 1944), Chap. II. 

(3) N. Bour and L. RoseNFELD: Phys. Rev., 78, 794 (1950). 
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will produce a similar uncertainty AH for that of the latter: 
(2) eTAE> Ap. 


There are several ways in which p(t + 7) can be measured. They all amount 
to letting the test body interact with some physical agent during a time 
interval At, in the course of which the position is altered by an unknown 
amount Av. From Heisenberg’s principle 


(3) Ap: Ar~wh, 


h 
(4) AE> STAR 
Clearly, the larger e is, the greater the accuracy. Nevertheless, it is impos- 
sible arbitrarily to choose a large value of e without augmenting another source 
of uncertainty. This would result from the fact that, during the time inter- 
val At required for the measurement of the momentum, the test body is 
accelerated and its energy altered, by an uncontrollable amount A€, due to 
the emission of radiation. A€ satisfies the inequality 


(5) Age >< 


(6) Ap> di: 


we obtain, for the components of the electric field, 


eAx 
(7) AE > GTA? 


The total uncertainty of the electric field is the sum of 4E and 4,E, and 


h eAx 
(8) e ro 


Replacing the right hand side of this inequality, considered as a function of e, 


x 
ni 
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by its minimum value, we have 


Vhe 
(9) AE + AE> TA 
A similar discussion could be carried out for the case of a component of the 
magnetic field, leading also to the uncertainty (9). 

The relation (9) is of the greatest importance. In fact, should it be taken 
literally, it would imply that the measurement of a field component is subject 
. to a limitation which is not accounted for by the quantum theory of the electro- 
magnetic field. In this theory field components are represented by operators 
defined at each point of space-time. Operators at different points are, 
in general, non-commutable; this signifying that the measurements of the 
corresponding physical quantities are complementary. No limitation, how- 
ever, occurs in the theory for the measurement of a single component at one 
specific point. 


1:2. — As was anticipated in the Introduction, BoHR and ROSENFELD (?) 
pointed out that the result of LANDAU and PEIERLS does not necessarily affect 
the consistency of the theory. In fact, as the theory of the free electromagnetic 
field is independent of the atomic structure of matter, the use of a point charge 
as a test body creates a situation which the theory is not meant to comprehend. 

Measurements of freé fields should be carried out by the use of apparatus 
consisting of continuous distributions of charge and/or currents supported by 
bodies of finite extension. So instead of the field components fay(%) (4) let 
us consider the averages 


nis 
(10) Bo) = [fuo(0) (da). 


v 


Here v indicates both a space-time region and its volume V-e7, V being an 
ordinary three-dimensional volume and 7 a time interval, and (dx) = dx;(da). 
For the sake of brevity, let us denote by # the quantity F),(v), i.e. the 
x-component of the electric field average. The measurement of H can be 
performed by placing a rigid test body with a charge distribution of uniform 
density @ so as to occupy the space region V at the beginning of the time 
interval 7, and by measuring the change of momentum of such a body: 


(11) pit + T)— p(t) =oVTE. 


(4) As customary, Greek indices assume the values 0, 1, 2, 3, while x denotes a 
point of space-time. 
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Note that, for the experiment to have meaning, the test body should not effect 
any appreciable displacement from its initial position; namely, at the time 
t + T, it should still be covering, more or less, the region V. To achieve this, 
a system of elastic forces must be made to act constantly on the test body. 
Therefore p(t + 7) in eq. (11) actually represents the final momentum plus 
the amount of momentum lost in process of overcoming the action of the 
elastic forces. 

Proceeding as in the previous section, let us consider the sources of un- 
certainty affecting the measurement. Corresponding to (4) we find 


h 
(12) AE> oda VT’ 
On the other hand, the uncertainty due to the emission of radiation is greatly 
reduced because of the replacement of the point charge e by a continuous 
distribution of density o. In fact, in the latter case, the strength of the 
radiation field emitted by the body will not exceed the order of magnitude 
o Ax, this being the order of magnitude of the current oAx/At, multiplied 
by At according to Maxwell’s equations. A rough estimate for the uncertainty 
of the momentum caused by the emission of radiation yields 


(13) A,p~ @VAxrAt, 
which corresponds to (6). Finally, in analogy with (9), we have 


he _ gdo At 
LA RAT ERE 


(14) ARSA 


Bour and ROSENFELD have shown that the right hand side of this inequality 
can be made arbitrarily small: if the first term is reduced by the choice of 
a large value for 9 Av, the second term also can be rendered negligible by 
choosing At sufficiently smaller than 7. This argument holds for any values 
of V and 7, however small. 


1:3. — From the elaborate analysis of the measurement of field components 
made in reference (2), we have selected only those aspects necessary to outline 
the contrast between the viewpoints of LANDAU and PEIERLS and BoHR and 
ROSENFELD respectively. 

We wish to re-emphasize that the discussion summarized in sect. 1°2 is 
valid only in the approximation for which the electromagnetic field produced 
by the test body is treated classically, and the atomic structure of matter is 
disregarded. 


88 E. CORINALDESI 


2. — Average field strengths in quantum electrodynamics. 


21. — In the theory of the free electromagnetic field, the commutator of 
two field strength operators, {(r) and f{{(a"), is a e-number, namely an ordi- 
nary function of @ and w#', containing % and c as sole parameters. Also a 
c-number is the commutator of the operators 3/}(v) and F(), which are 
obtained, analogously with (10), by averaging the previous ones over the two 
regions v and v' respectively. 

The physical meaning of the commutator of ae) and ye) has been 
interpreted by BoHR and RosENFELD (*) as follows. Two test bodies, set to 
cover v and v’ respectively, with the intent of measuring the field averages 
F,y(0) and F,,(0') are, during the experiment, sources of electromagnetic fields 
by which they interact with each other.. A rough estimate shows that such 
interaction induces reciprocal Hmitations in the measurements of the two field 
averages. Moreover, the product of the uncertainties AP,,(v):AF,,(v') (5) is 
seen to coincide with the absolute value of the commutator 


(15) [ so), ee!) = oop [tect ko ae Oy gk yplp— Ook yk — Onoky k 0) O(k 3) E(k )(dk) - 


waa) Jae’) exp [i(k, w—~«a’)). 


, 


In this formula, k? = k,k,, 6(€) is the Dirac function and e(k), function of 
the four-vector k, is denned by 
(16) e) ek)=—-1 k<0, 


ko being the zero component of k. 

Therefore, also in the case of the measurement of two field averages, theory 
and experiment are in complete agreement in the approximation of free electro- 
magnetic fields. 


2:2.- It was pointed out in sect. 1°2 that a system of elastic forces should 
be applied to the test body in order to prevent it from shifting appreciably 
from the initial position. Let this be achieved by the use of a system of 
latches and springs, which connect the test body to a sort of immovable cage. 
The latches Keep the test body rigidly fixed to the cage in the period preceding 
and following the measurement during which the springs remain un-stretched. 


(5) The electromagnetic fields of the test bodies, by which such uncertainties are 
estimated, are calculated on the basis of Maxwell’s equations. The possible impor- 
tance of quantum effects will presently be discussed. 
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The latches are lifted only during the measurement (i.e. during the time interval 
T), when the springs are left to act. Moreover, in order to reduce the impor- 
tance of the electromagnetic field produced by the test body during the 
measurement, let the cage bear a charge distribution capable of counterbalane- 
ing, in the rest period when the latches are on, that of the test body. 
During the time interval 7, the two charge distributions, being shifted 
with respect to each other, are equivalent to a polarization density Pipa eh 
the following this will be referred to, somewhat loosely, as «the dipole ». 
Thus, during the measurement, the test body is subject to the action of: 


a) the field to be measured, #; 


b) the springs, of the type — x & (& denoting the displacement from the 
initial position, i.e. from the cage); 


c) the electromagnetic field produced by the dipole, namely to a reaction 
of the test body on itself, ,gy- 


It can be proved that the reaction E,,, is equivalent to a force of the 
type x’é, if the field of the dipole is calculated classically and if £ is so small 
that É and the higher powers may be neglected. If, therefore, the elastic 
constant x of the springs is chosen to coincide with x’, the actions d) and ec) 
cancel each other, and can be completely disregarded throughout the discus- 
sion (6). However, the proof of this compensation has been based on the use 
of the Maxwell field produced by the dipole. As is well known from quantum 
electrodynamics, such a Maxwell field has no other meaning than that of an 
expectation value, to which fluctuations are superposed. Such fluctuations 
exist even in the absence of a Maxwell field, i.e. in vacuo. We could say, more 
intuitively than es tak that, although the expectation value of a free field 


compo ient f,,(7) is zero in the absence of radiation, the expectation value of 
the product x gas field components, say {®()-f((2'), is on the contrary dif- 


ferent from zero. This is due to the fact that alight quantum can be «emitted » 
by one of the field components in the product and «absorbed » by the other, 
this being possible even in vacuo. To estimate the order of magnitude of the 
fluctuations bearing on the measurement of a field average 7’,,,(v), let us con- 
sider the vacuum expectation value 


(17) 5 (520) ), Bool (ye => an ap | One k,+ 6, ol ko —by kK —Syokyko) d(k2) (Ak) - 
al f 4 
Lr | (dx) I (dx') exp [i(k, e—e')]. 


v v 


(5) If this.is done, the remark made in sect. 1:2 concerning the meaning of 
p(t + T) becomes superfluous. 
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Here the {} denote the anticommutator. Otherwise the symbols used in this 
formula are those of eq. (15). 

As for the physical meaning, note that, when the pairs of indices yw, v 
and o, o are identical, (17) represents the mean square deviation of 3'))(v) 
from the classical Maxwell value (here zero in absence of radiation), and, in 
the opposite case, expresses the correlation between 3/(v) and Too(0). 

23. - We may now proceed to the more practical question of evaluating 
the quantity (17) for some special region v. Let v consist of a time interval T 
and a sphere V of radius R. Then 


(18) fam) [ ae) exp riti eno (0, 


v v 


where 


(19) g(K) = | exp [i(k, x)] (dx) — = [sin (KR) — KR cos (KR)] , 


y 


depends only on K=|k|. Inserting eqs. (18) and (19) into (17), this latter 
is seen to be different from zero only when the pairs of indices 4, v and 9, o 
are identical. (Such a result depends on the spherical symmetry of V, the 
spatial part of v). Moreover, the mean square deviations of the various com- 
ponents of the electric and of the magnetic field are all equal to one another 
their common expression being the (convergent) integral (7) 


(20) S(v) = sia | sint (>) (q(K))?K AK. 

0 
Thus the mean square deviations are a known function of 7 and R. In the 
following, we shall not expand the study of such a function, but confine our- 
selves to its limit for 7 and R tending to zero. For the sake of simplicity, 
let us carry out the limit by first assuming c7 and R to be equal, cT= R= e, 
and by letting e tend to zero. It is then apparent that (8) 


(21) lim S(v) = co. 


e>0 


(?) E. CORINALDESI: Ph. D. Thesis, Manchester University (1951). 
(8) A similar result would be obtained by letting 7 or R tend to zero independently. 
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This result is of the greatest importance. Recalling the arguments of Part 1, 
we see that it affects the statement made in sect. 12, that the field strengths 
at a certaint point can be measured as the limit of a field average over a 
region v, the dimensions of which can be chosen arbitrarily small. But the 
smaller v is, the greater becomes the importance of the fluctuations, which 
mask the field to be measured. The present conclusion is not, however, 
identical with that reached by LANDAU and PEIERLS as to the non-measura- 
bility of a field component at a point of the space-times. In fact, the fluctua- 
tions are by no means foreign to the theory. On the contrary, they have 
been calculated on the basis of the theory, of which they are therefore a 
consequence. 


2:4. — There are situations which fall entirely into the domain where the 
classical electromagnetic theory is applicable, and where no quantum effects 
need to be taken into account in the discussion of measurements. This will 
be the case whenever the quantities (15) and (17), expressing the non-com- 
mutability and the fluctuations of the field components, are very small. 
Following BoHR and ROSENFELD (?), let us define a critical quantity Q as 
the square root of the absolute value of (15), assuming that v’ be of the same 
size and shape as v, but shifted with respect to v by an amount comparable 
with its dimensions. We shall presently express this concept more precisely. 
Another critical quantity G can be defined as the square root of (17). The 
classical theory is valid when the field components involved are all larger 
than Q and S. 

A rough estimate by BoHR and ROSENFELD seemed to indicate that © is 
always smaller than Q in all cases of physical interest, thus pointing to the 
latter as being the true critical quantity of the theory. This conclusion, how- 
ever, is not confirmed by a study of eqs. (15) and (17) (7), which are in a form 
particularly suitable for comparison. i 

Let v consist of the time interval 7 and a sphere V of radius R. Thus the 
non-vanishing components of (17) are equal to S(v) of eq. (20). For the sake 
of simplicity, we assume that v’ is obtained by applying to v the displa- 
cement represented by the four-vector J = (J,, J). Then 


bo 
bo 


(22) i (42) i (42) exp [i(k, @—a')] = 


v 


= exp [— ik, DI | (ae) | (ae) exp [u(k, ©—a')). 


As a consequence of this relation, 2S(v) is readily proved to be an upper limit 
to the absolute value of any component of (15). So the critical quantity Q 
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is smaller than V2 6, and the fluctuations seem to be the prevailing quantum 
feature. 

Since this point is of some importance, it can be verified for several types 
of regions. For example, if v is a four-dimensional parallelepiped of sides 
Lo, I,, Ls, L3, all correlations vanish, while the mean square deviation $, 
(i =1, 2,3) of the i-th component of the average electric field is equal to 
that of the corresponding component of the average magnetic field, although 
in general 8; 4 S;. It is also easy to prove that, if v' is obtained from v by 
a translation, the absolute value of any component of (15) is always smaller 
than S,+8,+8,. From such results it might be inferred that, while Q can 
' in some cases be much smaller than 6, in no case does it exceed an upper 
limit yS where y is a small number. 


2:5. — In the previous sections, measurements of average field components 
were discussed in their relation to the theory of the free electromagnetic field. 

The idea of a free electromagnetic field is, nevertheless, somewhat abstract. 
It seems to be generally accepted nowadays that the right theory will have 
to consider the electromagnetic field and the matter fields with which it interacts 
as a unity. Although the final theory is still to find, familiar perturbation 
methods can supply some idea as to whether the interaction with the matter 
fields bears on the problem of measurability. 

In the more general case in which the interaction with the matter fields (°) 
is taken into account, the field operators in the Heisenberg representation 
will be denoted by gothic characters but without the upper index (°) which 
characterized the operators of the free electromagnetic field. 

The first new feature to be encountered is that now the commutators of 
field components are not c-numbers, but instead operators, which are built 
out of the various electromagnetic and matter field components Occurring in 
the theory, and can be expanded in powers of the fine structure constant 
a = e?/he (19). The zero order terms of such expansions coincide with the 
corresponding commutators of free fields, as they should do. Among the 
terms of the first order in x, on the other hand, some are seen to be affected 
by tbe usual divergences of quantum electrodynamics. However, in the case 
of interaction with the electron and the scalar meson field, the special form 
of the divergent terms allows far their elimination by renormalization. For 
example we have 


(23) Lf, (0), fe) = 0), fe] + (CLAM), fee) + 2) +... 


(*) Such as the electron field, the scalar meson field, ete. 
(19) Assumed to be the same for the matter fields considered in the following. 
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where C is an infinite constant, i.e. a divergent integral, and Q is an operator 
free from divergences (7:11). Writing (23) in the form 


(23') [fo (0); Foo(@")] = (1 + oO) [fin (@), fe) +2 +... 
and carrying out the renormalization of the field strengths 


loi 
(24) We >(1 ae = ie ’ 


the commutator of the renormalized field components is free from the diver- 
gences occurring in the termsof the first order in «: 


(25) Ful); fo] = [{n,) fer + a2 +... 


It has been shown by Dyson (2) that in tbe case of interaction with the 
electron field all divergences can be eliminated by renormalization up to any 
power of «. 

In the following we shall constantly assume that the field strengths have 
been renormalized. 

We have already pointed out that £ and all the higher order terms in (25) 
are operators. Thus, in order to gain information as to the importance of 
such terms in problems of measurability, their expectation value should be 
calculated for that state of the fields involved which best corresponds to the 
conditions under which each experiment is carried out. For our purpose the 
vacuum can be assumed to be such a state. 

After these premises, we can, with the help of (25), evaluate the com- 
mutator of the (renormalized) field averages %,,(v) and %,,(v'). We find that 
the vacuum expectation value of such a commutator, <[%,,(v), Tos(0)]>, 
differs from (15) in that 6(k?) is replaced by È 


È i k2 + 4m? k? + 4m? 
(26) 0068) — pa |P LI sm fe ( 3k | 


1 ke Ay2\ 3/2 
+ 5 O(k? + 4) (— = [ex 


where 0(É) denotes the discontinuous function 
(27) ag) =0 £>O0, 0(É) =1 E<0. 


The symbols m and x denote the mass of the electron and of the scalar meson 


(11) W. PauLI: Feldquantisierung, Lecture notes (Ziirich, 1951). 
(12) F. J. Dyson: Phys. Rev., 82, 428 (1951). 
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multiplied by c/h (namely the inverses of the Compton wave lengths). We 
have confined ourselves to these two examples, although the formulae for other 
types of fields, such as the vector meson field, can easily be obtained. From 
(26) the expression (15) for the commutator of two average field components 
can be corrected to the first order in x, such a correction resulting from the 
simple process of emission and absorption of charged particles, electrons and 
mesons, in vacuo. Without entering into a more elaborate discussion, we 
shall content ourselves with stating that the above correction is negligible 
compared with the zero order term. Note that such a statement cannot be 
based only on the smallness of x; the behaviour and order of magnitude of 
‘ the integrals involved must also be carefully considered. 

As well as the commutator, the anticommutator of two average field com- 
ponents is also modified when the interaction with the matter fields is taken 
into account. It is easy to show that (1/2) <{%,,(v), (v)}> can be obtained, 
accurate to the first order in «, on replacement of 6(k?) by (26) in eq. (17). 
The correction thus produced, however, is genuinely small only when the 
region v is finite, while it can exceed the zero order term when one of the 
dimensions of v tends to zero. In the following section we shall give special 
attention to this point. 


2:6. — In sect. 2:3 the behaviour of the mean square deviations of average 
components of the free electromagnetic field was studied in the special case 
when the dimensions of the region v tend to zero. For a better understanding 
of Part 3, we shall now briefly discuss the behaviour of the first order cor- 
rections in similar circumstances. 


2°6.1. — Instantaneous space averages. — Consider first the case where V is 
a finite space region, and 7 is very small. For the x-component of the electric 
field average (the case to which we shall need to refer later), we have 


(28) lim <{Folv), Folv)}> = lim (o, log T + o = 


™>0 T>0 È 


SIZE ye 


Thus the mean square deviation of %,(v) diverges even if 7 alone tends 
to zero, while V remains finite; the divergence being logarithmic for the zero 
order term of the expansion in powers of «, and linear for the term linear in «. 
This result can be confirmed by considering the quantity 


1 he, 
= | (ds) f (02!) {fale 9, fale’, DY = Gr | (o(aey)exe ae — 


Vv Li 0 


(29) 
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che | iS a> /k? + 4m? k? + 4m? 
<A (410)? Cdk: (| nl y: ie a ) De 
0 V4 amt 
ae i k? + 4y2\92 _ \ 
Vrs au? 


V being a sphere of radius R and g(K) the expression defined in eq. (19). This 
formula gives the mean square deviation of 


(30) + | fos (2 8 (ax) 
r 

namely the space average of a component of the electric field at a certain fixed 
time, which is the same for all points of V. The first integral on the right 
hand side of (29) is logarithmically divergent, while the integrals over ky which 
occur in the «-correction are linearly divergent. Thus (28) and (29) agree 
with each other in fixing the degree of divergence of the fluctuations which 
affect the measurement of (30), a quantity referred to in Part 3. as an 
«instantaneous space average ». 


2:6.2. — Averages over. a thin space layer. — Let us assume that the spatial 
part of the region v is a thin layer of thickness b. For b very small the mean 
square deviations of %,,(v) increase logarithmically with decreasing b, even 
if the terms of the first order in x are taken into account. Such a result is dis- 
similar from that found in the previous subsection. To let the time dimension 
of v tend to zero does not produce the same result as letting one of the space 
dimensions do so —a characteristic asymmetry of the properties of space and 
time in the present problem. 


3. — Measurement of charge-current components. 


3:1. — In this and the following sections, we shall discuss measurements 
of a charge-current distribution j,(w). It has been made clear by Bour and 
ROSENFELD (*) that only the measurement of averages of the type 


(31) Tyo) ==] joa) a 


v 


is a physically well-defined procedure. As usual, v here denotes a finite space- 
time region (and also its four-dimensional volume). 
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By the use of Maxwell’s equations, the averages (31) can be expressed in 
terms of the electromagnetic fields produced by the charge-current distri- 
bution in question, 


DOSSI: 1 
2 nese ha 
(32) OE I sa (a) = 5 i fra dows 


where s is the three-dimensional boundary of v. In the special case is which v 
consists of a time interval 7 and a space region V, eqs. (32) assume the form 


1 Li 
(33a) Jo) = rr|* fe n)ds , 
(335 Jo) == Hays.) Biden 
) Whee arya fan - 77] e 


Here , indicates the (two-dimensional) boundary of V. The second term 
in (33b) represents the difference of the values assumed at the times 0 and 7 
(the beginning and the end of the time interval) by the space average, over 
the region V, of the electric field. 

The physical meaning of the formulae (32) and (33) is that the measu- 
rement of the quantities (31) can be reduced to that of the flux of the electro- 
magnetic field through the boundary of v. 

Bour and ROSENFELD (*) have shown how such a measurement can be 
performed by laying a system of test bodies on the boundary of v, and by 
controlling the momentum imparted to them by the electromagnetic field. 
For instance, in the case of (33a), a system of charges has to be placed on 
the surface S, and the increase of the component of the momentum in the 
direction n during the time interval 7 has to be observed. In the case of (335), 
things are more complicated. By a system of currents disposed on S only 
the first term will be measured. The measurement of the second term shouid 
be carried out by setting a charged test body to occupy the entire region V. 
Measurements of momentum increase should then be performed, one at the 
time 0 and the other at the time 7, to determine the two corresponding 
instantaneous values of the space average of the electric field over the region V. 
Two such measurements of momentum will take place during two short time 
intervals At, around t =0 and t = T respectively. 

There are statistical effects, however, which must be considered in the 
above measurements of the flux of the electromagnetic field through the boun- 
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dary of v. In analogy with what has been said in Pat 2, it must be stressed 
that, by the ideal experimental procedure outlined above, we measure the 
sum of (32) and of the flux of the fields which the test bodies themselves pro- 
duce. The field actions of the test bodies, on the other hand, fluctuate uarond 
their maxwellian values. Such fluctuations are connected with the elementary 
processes of emission and absorption of light quanta, and of creation and 
annihilation of electrified particles. 

It has been observed in Part 2 that the fluctuations due to the former 
of these effects exist even for a free electromagnetic field. They have no 
bearing, however, on the measurement of the field flux (32). In fact, it can 
be proved that if the quantity 


tb 
(34) 3 I {3,(0) VD) Io (v)}> a5 | do, | do; <{fux(#) » falw )}) 


v2 
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is expanded in powers of the fine structure constant «, the first term of the 
expansion is linear in x. Thus the zero order terms, which would exist even 
in the case of free electromagnetic fields, is identically zero. In the following 
we shall confine ourselves to the first non-vanishing term in the expansion, 
of (34), which we shall loosely identify with (34) itself. 

The commutator of two average charge-current operators %,(v), 3,(v') is 
also of importance, representing the complementarity of the two quantities 
involved. For the time being we wish only to point out that such a com- 
mutator vanishes in the zeroth approximation, in which the electromagnetic 
fields are considered as free. In fact, expanding the quantity 


1 A 7 "A , 
(35) [Fu(0), I,(04)] =w] ao, | dog [fua (2), Î,g(24)] , 
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in powers of «, the first non-vanishing term of the expansion is seen to be 
linear in «, for any regions v and v’. 


3:2. — Confining ourselves to the linear term in the expansion of (34) in 
powers of «, the expression 


(36) Tue = pa | (A) | (Ae) citta, OD, 


v 


can be easily obtained, where ji"! indicates the charge-current density in the 
interaction representation. 
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The calculation of $,, yields the expressions (1*) 
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The first eq. (37) refers to the case of electrons and the second to that of scalar 
mesons (mass me/f and ye/h respectively) (14). The various symbols occurring 
in (37a,b) have all been previously defined. 

The component So represents the mean square deviation of SW («charge » 
fluctuations — actually 9, has the dimensions of a charge density). It was 
first shown by HEISENBERG (!) in the electron case that the integrals by 
which So is expressed diverge as long as v is sharply defined, and that con- 
vergence may be secured by the contrivance of «smoothing out» v on its 
space boundary, say in a shell of width d. In Heisenberg’s paper, b was 
assumed to be larger than the Compton wave-length of tbe electron, h/me. 
The opposite case was considered by Jost and LUTTINGER (1%), who proved 
that, for b<h/me, Joo iS proportional to log (A/meb). A similar result can be 
obtained for scalar mesons. 

The smoothing out of the region v on its space boundary amounts to 
calculating Qo. for the time interval 7’ and for the various volumes V’ enclosed 
by the members of a family of surfaces 8S’, one of which coincides with S, 
and by taking the average of the results. Physically, it corresponds to the 
fact that the system cf charges laid on S in order to measure (33a) forms a 
shell, the thickness of which can be identified with 6. Thus when b tends to 
zero the charge fluctuations increase logarithmically with b for both electrons 
and scalar mesons. | 

We shall now consider the remaining components of $,,. For the sake 
of simplicity let us assume that v consists of a sphere and a time interval. 
By this particular choice all the non-diagonal ccmponents of %,,, become 


(13) E. CoRINALDESI: Nuovo Cimento, 8, 494 (1951). 

(44) The case of vector mesons will not be dealt with here. For it, cf. previous 
reference. 

(15) W. HEISENBERG: Leipziger Ber., 86, 317 (1934). 

(16) Unpublished; cf. reference (18). 
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equal to zero, while $1, $2 and $3, are equal to one another, their 
common value being in the following denoted by Syx. This represents the 
mean square deviation of any average current component 9; (i = 1, 2, 3) and 
will be referred to as «current fluctuations ». It is then easily seen that the 
integrals occurring in x. diverge when the region v is sharply defined. The 
present case, however, is somewhat different from that of So. In fact, the 
introduction of a smoothing out of the space boundary of v is not sufficent 
to render $x x finite (1°). It is necessary instead to smooth out the time 
interval 7, by assuming that a latitude b,/e be allowed in the fixation of the 
beginning and of the end of such an interval. Moreover Sx x is seen to increase 
linearly with c/bo when b, is much smaller than #/me or h/uc; namely more 
rapidly than Yo does for decreasing b. In order to trace the difference in the 
behaviour of Yo and Sxx, let us go back to the quantities (33a,b), and see 
whether it is possible to anticipate the behaviour of the fluctuations affecting 
their measurement, simply by inspecting their respective form. Writing the 
corresponding relations for the field operators, we have 


ue 
(39a) sx(0) = ral di | (Gn) dz, 
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It is then easily seen that, if V is assumed to be a sphere as has constantly 
been done, the first term of (390), which involves the magnetic field, gives a 
a vanishing contribution to the mean square deviation jxx. The second term, 
on the other hand, is actually responsible for the discrepancies between Joo 
and Yxx. It involves the knowledge of the instantaneous values of the space 
average of the electric field over V. In sect. 2°6.1 it was shown that the mean 
square deviation of an average component of the electric field over a region v 
diverges linearly if the time dimension of the region tends to zero (thus giving 
in the limit the instantaneous space average of the field component as it occurs 
in the second term of (39b)), and diverges only logarithmically when one of 
the space dimensions does so (thus giving in the limit the average over a 
time interval and a space region with one vanishing dimension, a shell, as it 
is the case in So and in the first term of Yxx). 

Considering the question from a more physical point of view, we recall 
that the measurement of the second term in (335) is carried out by setting a. 


(17) E. CORINALDESI: Nuovo Cimento, 9, 194 (1952). 


7 - Supplemento al Nuovo Cimento. 
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charged test body to occupy the entire region V and by determining the change 
of momentum of the test body first during a time interval At around the 
instant t = 0 and successively during a time interval At arount the instant 
t=T. It is clear that At is to be identified with the latitude b,/¢ allowed in 
the definition of the time interval 7. 

The dissimilarity between the behaviour of Yo and of Fxx might be of 
a certain importance in connection with an observation made by BoHR and 
ROSENFELD (?), who pointed out that even if the fixation of the boundary 
ot V is affected by a latitude v very small compared with the Compton wave- 
length of the electron and of the scalar meson, the charge fluctuations will 


.not have an excessive effect, as they depend only on the logarithm of the 


ratio between b and the Compton wave length. In the corresponding case, 
on the contrary, the current fluctuations might be appreciable. 

We may end by pointing out that the fluctuations are by far the pre- 
dominant quantum feature in the measurement of current averages. It can 
very easily be proved that the absolute value of the commutator (35), or 
better, of its vacuum expectation value can never exceed a quantity of the 
order of magnitude of the larger of the mean square deviations Yoo and Fxx, 
a result which is similar to that found for the electromagnetic field components. 
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Some Wave Propagation Properties of the Electron Beam 
in Vacuum and in Ionized Media (*). 


O. E. H. RyDBECK 


Research Laboratory of Electronics, Chalmers University - Gothenburg, Sweden 


(ricevuto il 29 Agosto 1952) 


1. — The properties of the electron beam as a wave propagating medium has 
attracted considerable attention recently. In the latest types of micro-wave 
tubes, the travelling wave tube (TWT) and the electron wave tube (EWT), 
the electron beam is used as a wave propagating medium which, if suitably 
connected to input and output circuits, permits a very wide frequency band 
to be transmitted. Through interaction with periodic structures, for example a 
helix as in the TWT-case, a periodic d.c. potential field etc. or with an other 
electron beam as in the EWT-case, considerable amplification or gain can be 
obtained. For further details see for example The Theory of the Travelling- 
Wave Tube, by O. E. H. RyDBECK, Reports from the Research Laboratory of 
Electronics, No 1 (RRLE 1) (1948), and On the Theory of Electron Wave Tu- 
bes I and TI by O. E. H. RyDBECK and S. K. H. ForscreNn, RRLE 15 (1951), 
RRLE 31 (1953). 

In the auroral ionosphere streams of electrons appear as in the solar corona. 
From the general wave propagation point of view the influence of the electron 
stream upon the radio wave propagation in ionized media should therefore 
be investigated. In order to study the mechanism of the outbursts of solar 
radio noise such investigations are of fundamental importance. 

New properties of the electron beam as a wave propagator even in the 
millimeter wave length range are now being studied at the Chalmers Research 
Laboratory of Electronics and at several other laboratories especially in 
the USA. 

It is the purpose of this lecture to give a concentrated survey of the fun- 
damental wave propagation properties of the electron beam, its behaviour 

(*) Conference held in Naples at the Meeting of ionospherical studies (see page 176 
of the Supplemento n. 2 del Nuovo Cimento, 9 (1952). (Editor's Note). 
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in ionized media and some of its recently discovered properties in inhomo- 
geneous media, which are of special interest to the ionospheric and micro wave 
physicists. 


2. — For the sake of symmetrical simplicity let us study a cylindrical 
electron beam of radius a travelling with an average speed v, in an evacuated 
wave guide of radius b (b > a). The beam travels in the centre of the guide 
along its axis which we denote 2. As we study the steady state the combined 


time-factor and axial propagation factor can be written 


exp [j(wt — yz)] , GN 


where is the angular signal frequency and y the axial wave propagation 
coefficient (y = 27/4., where A, is the axial wave length). In order to simplify 
matters even further, as the time for the presentation of this material is very 
limited, we assume that the electrons only oscillate in the axial or z-direction. 
This means that we have a strong axial, focussing magnetic field. 

If è is the current density of the beam and e the electronic charge, the 
electron density, N,, of the beam is N,——di,/ev,. The angular critical fre- 
quency of an ionized medium with this electron density, as we would like 
to introduce the conceptions of the ionospheric wave propagation, thus is 
On — VN,e2/mey, where m is the mass of the electron and £ the dielectric 
constant of free space. As the electrons can not oscillate in the transverse 
direction the transverse dielectric « constant » is e, and independent of the 
electron stream. It should be added in this connexion that static space charge 
effects are neglected in the beam, that is it travels with constant speed and 
cross section. 

The dielectric « constant » in the axial direction, e. is 


eal 


(0 Yo)? J 3 


(1) REZZA 


as is easily shown. This is a very important expression which, it should be 
added, only holds when the alternating velocity, v, (ac-velocity) is small com- 
pared with v (the de-velocity). We thus restrict us to the «small signal 
theory » and avoid non-linear effects. One readily sees from (1) that when the 
beam drift stops (v = 0) we obtain the regular dielectric constant of the 
ionized medium. 

If the space charge wave travelling along the beam (due to the oscillations 
of the electrons around their average drift positions) travelled with a phase 
velocity equal to v (which in practice does not occur, see (1)) the axial wave 
propagation coefficient should be w/v) which we denote by « (i.e. x = w/v). 
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As, however, the space charge waves, except when they are violently perturbed, 
almost travel with a phase velocity equal to v it is convenient to introduce 
the notation y = «— Ay, where |Ay| usually is much smaller than «. 

If E. denotes the axial electric field strength of a given mode or type of 
space charge wave, the corresponding ac-velocity is 


(2) va = jeE.Jom-a/(a—y) - 
Y 


A 
The corresponding ac-current (produced by the space charge wave) becomes 


ADE x 


8) = i%a, (dn 


where z, =: Vo/îo is the specific de-impedance of the beam (Vo is the acce- 
lerating voltage giving the beam the average velocity vo). 

When all E,-modes, and their corresponding y-values have been found, 
the magnitudes of them, or their degree of excitation are determined by the 
kind of velocity modulation impressed upon the beam. When the mode 
amplitudes have been found the ac-current is easily determined along the 
beam and its cross section by relation (3). 

Before we proceed further it should be added that « usually is much larger 
than k, = 22/4), where 4, is the vacuum wave length. If the electron beam 
travels with 1/10 of the speed of ligth x/k, = 10. Usually, therefore, y > ho. 

Introducing cylindrical coordinates 2, 0, and 9, the TM-modes (*) (trans- 
verse magnetic) travelling on the beam become (if we for convenience omit 
the time and propagation factor exp [j(ot — y2)]) 


Ez — A,I,(he) ; E, = ibi “Al, (ho), 
0 


(4) } and 


. D ko 


Le pasto) ’ (@ <4) 


where hy = V2 — Ke > ko, h= Dh, D = €:/50, Zo BV ie) eas the intrinsic 
impedance of free space and I, and I, are modified Bessel-functions of orders 0 
and 1. 


(*) Because the electrons have no transverse oscillation, the TE-modes are not 
perturbed by the electron beam and therefore are not discussed here. 
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Outside the beam the field components obviously become 


(5) E. = Ax{Ko(hoo) + CI(Mb0)}, ete. (b>o>a) 


A, and C are determined by the boundary conditions that the radial wave 
impedance, E./H, should be continuous for o =a and b. If the wave guide 
tube is a perfect conductor, then E,(0 = b) = 0. 

Because J,(v) increases exponentially with w one infers that as b tends to 
infinity (very wide wave guide), i.e. we approach the important case of a 
cylindrical electron beam in free space, C tends to zero. 

As h, is so much larger than k, the influence of the guide walls disappears 
when d is of the order of a few vacuum wave lengths. For the very wide 
tube case, or the electron beam in free space, we thus have the following 
boundary relation 


I 
(6) = ha Pay has 23 


This relation can always be satisfied, if h, i.e. if D, becomes imaginary. And 
this is easily possible if only Ay is small enough so that e,/e) becomes negative. 
Introducing D = jQ, and H=Qh,. the boundary relation (6) transforms to 


(6a) TE ZE = Ip E, 


This relation only holds when h, is large, and real, that is for waves, space 
charge waves, attached to the beam and roughly travelling with the same 
velocity. The regular guide waves, which usually travel much faster than 
the beam and thus hardly interact with its electrons, are also contained in 
the general boundary relation (C0). For them normally also h) is imaginary 
(i.e. y< ky) and the boundary relation gets a form different from (6a) and 
only involving normal Bessel- and Neumann-functions (for further details se 
introductory reference 2). As the slow electron beam only slightly perturbs 
the regular TM-guide waves we do not need consider them here. In a fol- 
lowing section, where we study the properties of the electron beam in an 
ionized medium, we will have occasion to return to them in connexion with 
the study of the coupling between space charge (or electron beam) waves and 
electromagnetic radiation. 

If the n-th solution of (6a) is H,a, one easily infers that for thick beams 
H,a = W° which is the n-th zero of Je). For Q we use the convenient 
form 


lod 3 / (01/00)? Ba / PACE 
e Sr e (35) di 
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where k, = 2af,/v) = 22/4, and A, is interpreted as the plasma wave length 
of a plasma oscillation propagating with the phase velocity v of the beam 
(fi = @,/22). It appears from (7) that the plasma propagation coefficient È, 
is a very important parameter. In fact it characterizes the behaviour of the 
electron beam as a wave propagating medium. 

For the n-th mode we thus obtain (treating ho as constant ow Va? — k?) 
from (6a) and (7) 


(8) = yan (Hay 


For the lowest order modes (n = 1, 2,3 say) we see from (8) that for a 
thick beam Ay, 93 = + ky. 


As y> ko, bn ~y =, because Ay is small. Therefore roughly 4, = 
= ho'tolto. A thick electron beam, from the wave propagation point of view, 
thus must be many axial wave lengths thick. According to what has just 
been shown, however, this does not mean that it has to be many vacuum wave 
lengths thick, as co > vo. As Ko(hog) decreases exponentially as exp [— Wo] 
the electric field connected with the space charge wave falls off about as 
exp [— ho(o — @)] = exp [- 22(0 — a)/à.] outside the beam. For a wave 
length of 3 em and a velocity ratio colto = 10 (a customary value) 27/4. = 
— 20 em-1, i.e. the amplitudes of the associated fields fall off very rapidly 
outside the beam. The electric and magnetic fields associated with the space 
charge wave (or electron wave as it is also called) thus are mainly confined 
to the beam. For interaction to take place between two electron beams, as 
in the EWT, or between an outer field and the beam, the outer beam or the 
outer field must be brought very close to the beam. This means that in the 
EWT the two interacting beams must be very well mixed. 

In order to sum up the excited modes the modulation conditions’ must be 
known, as has already been stated. If we assume, that velocity modulation 
takes place with the same amplitude across the beam, i.e. we assume the beam 
to pass through an ideal pair of non-absorbing modulating grids (which may 
be connected to a micro wave cavity in some fashion or other), the sum of 
all excited modes must yield a constant ac-velocity across the modulation 
plane, which we put at 2 = 0. If we denote the amplitude of the ac-velocity 
atz=0 by 00 (vo <%, ie. vo = AV/2V,, where AV is the amplitude of 
the alternating (micro wave) voltage across the modulating grids) the total 
alternating beam current 7; becomes 


Wr(a) 
J,(H 9a) 
: v ; © sin (Ayng) H,a 
ili = 2 -O-exp [— j(0 + 2/2 Se) . (11: DI) St? 
(9) t/t 7 0-exp [— j(0 + a] > RG o(Hn0) T5H,a) + JiH,a) 
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where 0 is the average transit angle x. The corresponding total ac-velocity, 
vi, similarly becomes 


(10) vi/o° = 2 exp [— j0] 5 cos (4y,2):J(H,0)-y(a). 


n= 1 


From the above expressions one finds, that close to the modulating grids, 
where (Ayz)? < 1, 


yo 
(9a) dalio = n ‘0-exp = j(0 = m/2)] bd [(Ayn2)? < a £ 


This is the ac-klystron current according to the purely geometrical conception. 


We note from (9) and (10), however, that the strong bunching of electrons 
according to the geometrical theory (faster electrons overtaking slower ones) 
does not appear. The mutual repulsion between the individual electrons 
makes the entire beam elastic. All electrons oscillate around their mean velo- 
city with the plasmic frequency, or rather very nearly the plasmic frequency as 
it would be for a very thick beam. This means that every mode has its planes 

located m:22/Ay, apart, where 
Diameter of beam ——— | the ac-velocity and the ac-cur- 
I rent are the same respectively. 
It is also to be noted that the 
accurrent and the ac-velocity 
are out of space phase. Finally 
it should be noted that rela- 
tions (9) and (10) very fully 
describe the elastic wave pro- 
pagating property of the elec- 
tron beam. 

The formal difference bet- 
ween the space charge wave 
and the regular wave propa- 
gated in the empty wave guide 
mainly is that the former be- 
comes so to speak amplitude 
modulated by the plasmic 
oscillations at the same time as it propagates with the phase velocity vo. 

Fig. 1 and 2 show the cross sectional variation of the axial electric field, H, . 


a=0.5em;A=3 cm; v=0.1005 220; is=1mA/mm? 


Fig. 1.- The first three #, beam wave modes as 
functions of radial distance. 


3. — In the auroral ionosphere or in the solar ionosphere the electron (or 
ion) beams travel in media with variable properties. In such media one can 
not simply predict what kind of axial propagation will take place. The 
propagation depends to a great extent upon the variation of the electron or 
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ion densities, or the beam velocity with the coordinates. If the potential 
varies, for example (and in this connexion we mean the accelerating potential), 
the plasmic beam will be 

compressed or elongated, if +. Diameter of beam __ 

we try to describe it in ele- / Vee 

mentary terms. 

Before we approach the 
problem of the propagation 
of the space charge wave in 
a medium with variable pa- 
rameters, let us study briefly 
what happens in an ionized 
medium. In order to use 
the same or similar notations 
as before we still stick to the 
cylindrical system. Let us 
further for simplicity assume 
that 1) the beam completely 
fills up the guide tube and 2) that the guide tube is filled with an ionized 
gas of electron density N.. Collisional effects are completely neglected in 
this brief study. 

The transverse dielectric constant is e, as before (we still have the strong 
axial focussing field which acts also on the electrons of the gas) but the longi- 
tudinal dielectric constant now becomes. 


Anp=3 CM; %=0.1 00; 220: H=1mMA/mmMm?, 


Fig. 2. — Variation of total axial electric strength 
across the beam. 
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where w: = VN ,e2/egm is the angular critical frequency corresponding to the 
electron density N,. 

As the beam completely fills up the guide (a = d) we have for the deter- 
mination of y, the folllwing relation, viz. 


Qha = W, 


where as before W° is the n-th zero of J,(v), and (with the new value of e.) 


ee Hones 
(11a) Q=| dog (=) as 


If we for the moment consider h, as a constant (i.e. when Ay is small) we obtain 
the following relation for Ay, viz. 


(12) ini Vi x (eye a 
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From this relation it appears that when (@./0)? > 1 + (W?/ha)?, Ay, becomes 
imaginary and amplification is possible. For a thick beam, and that is what 
we usually deal with in practice, this means that strong amplification is pos- 
sible, when (@,/m)? is slightly larger than one. For the ionospherist this is 
especially interesting to note because in the free ionized medium regular 
electromagnetic waves cannot propagate when (w/w)? is larger than 1. Ti 
we insert (12) in the expression for the axial dielectric « constant » (11) we 
find, for the various modes, 


*(11b) Ga SSS (4 ; 


hya 


This also serves as an explanation why the propagation is possible even when 
(@./m)? is very large. 

We have thus found the very important fact that an electron beam can carry 
waves and information (of course only with the average speed of the beam) 
through a heavily ionized medium where no regular electromagnetic waves 
can travel. We will see in a later section, which surveys the case when the 
focussing axial field has an arbitrary strength, that on the whole this remark- 
able property of the beam is practically independent of the field strength. 
Thus an electron bean passing through a heavily ionized auroral ionosphere 
for example can carry information through the layer at a carrier wave fre- 
quency where no regular electromagnetic wave transmission is possible. This 
equally well applies to an electron beam in the solar ionosphere. 

However, one important question remains to be solved. Does any coupling 
take place between the beam and electromagnetic waves, as it approaches 
regions where regular electromagnetic wave propagation is again possible? 
This is a very important problem if we try to explain the ejection of solar 
noise radio outbursts. 

Fortunately already relation (12), which only holds for small Ay-values, 
gives us a clue to the solution of the problem. As (@,/@)? approaches 
14 (W2/h a)? the amplitude of dy becomes very large, so large that A in (12) 
can no longer be regarded as a constant. Relation (12) then becomes a quartic 
in y, its exact form, which has many peculiar and interesting properties, which 
will be described graphically a little later. 

But, before we proceed, let us study briefly what kind of space charge 
and electromagnetic wave equations we get in the general case when we 
suspect that coupling occurs and we therefore cannot give the axial pro- 
pagation factor the simple harmonic or symbolic form exp [— jyz]. For the 
sake of simplicity we still assume that the waves are travelling in a cylindrical 
tube. One can show, that under these circumstances, for a given radial mode, 
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the transverse magnetic field gets the approximate form 


(13) H, = B-J,(H, e/a): Fe), 

where F(z) is the axial propagation function to be determined. If coupling 
appears, F(z) necessarily will be of complex nature. Now one can show, 
although the time does not permit me to demonstrate it here, that we obtain 
two coupled wave equations for the space charge wave and the electromagnetic 
wave F(z), viz. 


wave equation of electromagnetic wave — i 


(14) {i (02/0)?) ee 3) — Ei} He alee 


and 


wave equation of space charge wave 


az ke 3/1 dk\2 1 &k lee k, 
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where for the sake of simplicity ©? has been given the form wî = wo” exp [de]. 

We immediately notice that as we approach the critical level ©} = ©? 
(critical frequency level of the ionized medium) the coupling between the wave 
equations becomes terribly strong. They are in fact completely distorted. 
A more detailed analysis will show, that energy from the space charge wave 
(which is continuously growing in amplitude as we approach the critical fre- 
quency level from the deep ionization side, as for example in the solar iono- 
sphere) is transformed more or less completely into electromagnetic wave 
energy at and around this critical frequency level. The coupling is somewhat 
similar to that taking place at the ordinary critical frequency level in the 
magnetic ionosphere and there produces the triple split. 

If we regard the parameters w; and ©} as constants equations (14) and (15) 
are easily transformed to a fourth order wave equation having the roots of (12). 
At this moment it is therefore especially interesting to look at a graphical 
representation of its roots as functions of the ionization expressed as w3/@? (=X). 

Deep down in the ionized medium, X,> 1, the electron beam proceeds 
towards the critical level X, —1. The lowest order space charge mode (and 
the most important one) consists of two waves, labelled 1 and 3. Their Ay 
are complex conjugates, as long as X, is somewhat larger than 1. Wave 
number 3 continuously builds up its amplitude (it must be remembered that 
we are all the time thinking in terms of a small signal theory so the equations 
do not show how kinetic energy from the beam is converted into space charge 
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wave and later into electromagnetic wave energy) whereas wave 1 rapidly 
decreases in amplitude and looses significance. At or near the critical coupling 
or transformation level (quite close to X, == 1) wave 3 transforms over into 
an electromagnetic wave of the regular wave guide type. It is further very 


space charge waves 
domain (plasma) 


electromagnetic waves 
domain (wave guides) 


these waves travel in negative dire 
tien, that is, against the stream 


x=1.6 cm; A,=3 em; v=0.1 00; %=1MA/mm?; n=1; Jm= imaginary; Re= real. 


Fig. 3. — Waves roots of equation (12). 


interesting to note, that if deep down in the ionized medium an electro- 
magnetic wave could be excited (one notes from the field strength expressions 
that such a wave would have practically no axial electric field strength com- 
ponent for X,>> 1 and thus travels practically unaffected by the electrons 
which are not allowed to oscillate in a transverse direction in this simplified 
model) it would travel first with the speed of light, as wave number 2 in the 
diagram. At the coupling level, however, it very smoothy transforms into a 
space charge wave. As its axial electric field strength component develops 
it becomes retarded by the beam which it modulates. However, it is to early, 
in the present stage of the theory, to judge if this kind of wave is of signifi- 
cance or not in the solar or in the auroral ionosphere. Finally it should be 
stated that the complete solutions also contain an electromagnetic wave 
running against the stream and mainly affected by the stationary ionized 
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medium. This wave, the fourth solution of the quartic (12), is labelled 4 in 
the diagram (fig. 3). 

It is interesting to note from the coupled wave equations, that the coupling 
disappears when a tends to infinity, that is in the infinite plane case. A 
further study of the case indicates, that a radial inhomogeneity, in this case 
represented by the tube boundaries at 0 =a is essential to produce the 
effect of coupling. We are at present investigating if a radial inhomogeneity 
of the ionization or of the axial magnetic field would produce a similar, strong 
effect. Such an effect should be very important as a possible explanation 
of the generation of solar outbursts of radio noise. 


4. — So far, we have for the sake of simplicity limited ourselves to the case 
of the strong axial magnetic field, which prohibites the electrons to perform 
transverse oscillations. In spite of this simplification the problems have 
become complicated enough. 

However, in order to demonstrate to the ionospherist, how the magneto- 
ionic waves (in the longitudinal case to which we restrict ourselves) become 
perturbed by the electron beam, we sketch briefly, without going into the 
fairly laborious mathematical details, what happens generally. 

Because the electrons now also oscillate, more or less freely, in the trans- 
verse direction, TE-waves also are excited by the beam or propagate slowly 
along it. The axial electric field strength now gets the form 


(16) E, = AJ,(xo) + BI (YO), 
where 2 and y are related to each other by the relations 


(X, = 08/02, X2 = 3/0") 


x = g,(X,, X2, Ho) + gx(X1, X2, Ho) (H, = static magnetic field) 
and 
k2 
an) a= iq la x} 
Py 
ditt ee Vi al 
il (o — yvo)[(0—- PV) + Oy] @(@0+ wa) ki f 
P, 
oe Set aio a a Si vay 
| (0 — yvo)[(@ — yvo) — Og] o(@— Wp) ke f 
Ps 


(w, = the angular gyro frequency). 
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Because H, differs from zero in this case (both TM and TE modes are 
excited) it is a very difficult procedure to determine y from these relations 
(both E, and E, must be zero at the tube walls). If we, however, limit our- 
selves to the case of a very wide or rather infinitely wide tube, x and y are 
zero and we are reduced to the case of plane waves. 

The roots of the possible waves thus are determined from the relations 
P,= 0, P, =0, and P; = 0. First of all it is interesting to note, that the 
static magnetic field, represented by the angular magnetic gyrofrequency is 
absent in P, which thus (compare (12)) yields the same space charge wave as 
before independent of the electric field. It is further very interesting to note 
‘that P, and P; = 0, for v == 0, i.e. for stationary electrons, yield the «- and 
z-components (separated in critical frequency by the gyrofrequency) of the 
magneto-ionic theory. When v 4 0 we obtain x- and z-waves perturbed by 
the electron beam traversing the ionized medium. Further numerical studies 
of this perturbation, which have recently been done at our institute indicate 
that the perturbation normally is quite small even near the critical level 
X, = 1, except when the medium is inhomogeneous in a transverse direction. 

The TE waves thus produced in the medium form the perturbed electro- 
magnetic waves of the magneto-icnic theory. At the same time, when ana- 
lyzing P, = 0 and P3 = 0, one finds (because they are quartics in y) that 
they also have two solutions, for which y — a. These solutions represent 
transverse space charge waves, which follow the beam. There are two such 
waves associated with the #-wave and two such solutions associated with the 
z-wave. 

Thus a total of 10 different waves are produced or possible. In the posi- 
tive or drift direction: two longitudinal space charge waves (the important 
type of space charge wave we studied at the beginning), four transverse space 
charge waves (all these waves travel with about the speed of the beam), one 
perturbed x-wave practically with the speed of the magneto-ionic theory, one 
perturbed z-wave, practically with the speed of the magneto-ionic theory. 
In the negative, or anti drift direction: one very slightly perturbed x-wave 
and one very slightly perturbed z-wave, both very closely with the speeds 
of the magneto-ionic theory. 

This survey has shown the complex groups of waves produced by an electron 
beam travelling along the static magnetic field in an ionized medium. If the 
beam at a given level is assumed to have a certain amplitude versus time 
variation, as we assume in the solar ionosphere for example, the Fourier- 
transform of which is known, one can, at least theoretically, study how the 
disturbance spectrum travels through the medium and changes form as it 
approaches the critical level, where coupling into electromagnetic radiation 
takes place. 


SOME WAVE PROPAGATION PROPERTIES OF THE ELECTRON BEAM ETC. US: 


5. In conclusion a few words should be said about the behaviour of the 
electron beam or rather the space charge wave in vacuum when the acce- 
lerating potential varies along the path or direction of propagation. This is 
also a case of inhomogeneous wave propagation, because the phase speed of 
the beam varies along the axis of propagation. 

One finds in a way, similar to the deduction of the coupled wave equations 
(14) and (15), that in the case of plane wave propagation (i.e. in the case of 
a thick beam), when the beam current is constant, the axial propagation 
function gets the form 


18) E, = A II) exp [—jW] 


si 3/2 
Vo” 


where //(z) is a solution of the special space charge wave equation: 


all k2 om eal 2 : 1 2y 
(19) ; i (ii ro. 
“a Vo 


de? 1— X, dz dz 


To make the results more general we have included a constant ionization, 
represented by X,. 

If we neglect the ionization, and only study the original case of space 
charge waves in a medium of variable accelerating potential (this corresponds 
more or less to the propagation of regular waves in a medium of variable 
dielectric constant), it is convenient to introduce an exponential velocity change 
(which we assume can be obtained by suitable electron - optical arrangements). 
If we, for the sake of simplicity, put v = v exp [fz], equation (19) gets the 
simpler form 


(19a) all ee 1 pe) 77 = 6), 


From this equations one readily infers that strong amplification is possible 
if B > 2k,/3, i.e. if B > 47/34,. This means that the velocity change must be 
large within a distance shorter than one plasma wave length. In order to see 
what this means for the technical use of this new effect in the development 
of broad band micro-wave tubes we calculate the plasma wave length for a 
typical beam. If the current density is 1 mA/mm?, and the accelerating voltage 
1000 volts, the plasma wave length becomes about 12 em (01/27 & 150 MHz 
and N, x 3,35-108 electrons/em*). Thus it is technically normally quite 
simple to make use of this effect, which varies but little with the wave, 
or beam modulation, frequency. It can be shown from relation (18) as well 
as from the original wave equation (19) that the wave for maximum effect 


114 O. E. H. RYDBECK 


should run into lower velocity (6 negative), that is it should suffer a com- 
pression. This compression, if done as rapidly as required by the condition 
just discussed, immediately increases the ac-current amplitude. 

A further study, applying this principle to our original space charge wave (9) 
which has been produced by velocity modulation at a pair of plane grids, 
shows that the velocity jump (v/v),) Should be inserted at the ac-current node, 
i.e. at a point where Ay, — 2a, for maximum effect, as was originally 
shown by L. M. FreLp and co-workers (U.S.A., 1951). At such a point the 
ac-velocity is maximum (compare relation (10)) and the effect, in approximate 
elementary terms is to increase the ac-velocity from v° to 0° Volvo, The ap- 
proximate amplification thus becomes (vy/vo,)!?. If the beam runs from 1600 volts 
into a 100 volts region, the ac-current increase therefore is of the order of 
8 in amplitude. 

For many reasons it should be advantageous to raise the beam to the 
original velocity level after amplification. One could then, for example, run 
it through a new voltage step and obtain additional amplification. Further 
it is convenient to raise the beam to its original velocity for good demodulation 
at the end of the tube (by a helix or a cavity). Fortunately this can be done 
practically without loss of gain (without loss of amplification) if the voltage 
step, now bringing the beam up to its original velocity, is inserted at a place 
where the ac-velocity is zero, i.e. at a current maximum. Similar effects are 
also produced if the beam is run through properly arranged periodic poten- 
tial structures. 

At present we are testing this new effect in our tube laboratory in Gothen- 
burg. It is probably the simplest mechanism of electronic amplification in 
nature. It varies little with frequency and may prove to be very important 
for the future development of electron beam tubes for the amplification of 
millimeter and shorter waves with great bandwidths. 

I hope that this lecture has given you at least some idea of the fundamental 
properties of the electron beam and its space charge waves, which probably 
will be of vital technical importance in the telecommunication of the future 
and in the vast field of radio wave propagation in the auroral and solar iono- 
spheres. 


Finally I wish to express my warmest thanks Dr. MARIA Morzo for her 
generous help and assistance in preparing the lecture in manuscript form. 
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1. — Introduzione. 


In questo lavoro ci proponiamo di presentare un nuovo bilancio energetico 
della radiazione cosmica a latitudini molto settentrionali (intorno ai 50° N). 
Negli anni scorsi sono stati presentati calcoli dettagliati di questo tipo (!-?), 
o stime (*) intese a mettere in evidenza lo stato di bilancio; stato non molto 
soddisfacente in quanto sussisteva un notevole divario tra il valore numerico 
del flusso di potenza in ingresso nella atmosfera trasportato dalla componente 
primaria, e il flusso di potenza dissipato nella atmosfera dalle varie compo- 
nenti secondarie nei diversi processi in cui dette componenti secondarie ven- 
gono gradatamente assorbite. Vi sono ragioni che consigliano di rifare a questo 
punto un bilancio dettagliato; infatti la nostra conoscenza delle varie compo- 
nenti e dei processi da esse subiti ha assunto una discreta consistenza; vedremo 
che la situazione è notevolmente migliorata rispetto ai tentativi precedenti. 
Nel calcolo di Rossi infatti, del 1948, si ragionava su un solo tipo di mesone 
(la particella u) perchè il passaggio 7 + non era ancora conosciuto, e si con- 
siderava per la detta particella una disintegrazione binaria in elettrone e 
neutrino, u >e + y, per quanto una analisi dettagliata della componente 
molle al livello del mare suggerisse per la frazione di energia totale dei u, che 
in seguito alla disintegrazione, va in componente elettronica, un valore minore 
di 1/2. Tuttavia la scelta del valore 1/2 fatta da Rossi rappresentava allora 


(1) B. Rossi: Rev. Mod. Phys., 20, 537 (1958). 

(2) G. PupPi e N. DaLLaPORTA: Progress in Cosmic Ray Physics, Vol. I, Cap. VL 
(North Holland Publishing Co., Amsterdam, 1951). 

(3) J. A. VAN ALLEN e S. F. SINGER: Phys. Rev., 78, 819 (1950). 
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un ottimo, tenendo anche presente il fatto che, ragionando su un solo tipo 
di mesone era giustificato attribuirgli un decadimento binario. Per queste 
ragioni dunque il calcolo di RossI veniva a sottostimare sensibilmente le per- 
dite «invisibili ») in neutrini. Nel bilancio di Puppr e DALLAPORTA del 1951 
venivano apportate le dovute correzioni al bilancio di Rossi tenendo conto 
appunto del fatto che i mesoni u non vengono generati direttamente nella 
atmosfera, ma provengono dai mesoni x per disintegrazione, e che il decadi- 
mento dei u è un processo ternario, u —e + y + vy, e quindi in media solo 
1/3 della energia dei u va all’elettrone. Si otteneva così rispetto al bilancio 
di Rossi un aumento circa del 25% sul valore del flusso di potenza dissipato 
nella atmosfera dalle varie componenti (*). 

L’accennato miglioramento non era però sufficiente a sistemare il bilancio 
in quanto il flusso di potenza in ingresso, che si era andato sempre meglio 
definendo risultava notevolmente più elevato. Dalle esperienze condotte infatti 
e con razzi (8) e con palloni (45), pur con qualche riserva su un sensibile 
effetto di albedo, si giunge a definire un valore numerico per il flusso di po- 
tenza in ingresso, all’incirca doppio. Questa situazione non soddisfacente, sotto- 
lineata anche da altri autori (45), merita di essere analizzata molto accurata- 
mente, perchè, se reale, può significare una effettiva e grossa lacuna nell’at- 
tuale quadro descrittivo dei raggi cosmici e delle particelle elementari. 

Nella costruzione del bilancio intervengono infatti tre elementi sostanziali. 

Il primo di tali elementi è la conoscenza della composizione e della inten- 
sità della radiazione primaria in quanto la conoscenza dello spettro primario 
permette di definire il flusso di potenza in ingresso. 

Il secondo elemento è la conoscenza della natura e della intensità delle 
varie componenti secondarie che hanno origine dalla primaria, per potere, in 
base ad esse, calcolare il flusso di potenza dissipato dalle componenti. 

Il terzo elemento è la conoscenza delle modalità esatte dei vari processi 
in cui, da una parte la primaria genera le componenti secondarie, dall’altra 
le dette componenti secondarie si trasformano l’una nell’altra o vengono gra- 
dualmente assorbite. Poichè i primi due elementi sono di natura sperimentale, 
è ovvio che, ove si riuscisse a definirli in modo soddisfacente, il problema del 
bilancio potrebbe diventare critico nei riguardi del terzo elemento, e precisa- 


(*) Rossi ha condotto, indipendentemente da Puprr e DALLAPORTA una revisione 
del suo bilancio del 1948, e lo ha esposto in un seminario tenuto nell’estate del 1951 
a Padova. I risultati coincidono perfettamente. 

(4) J. R. WINKLER, T. Stix, K. DwiGHT e R. SABIN: Phys. Rev., 79, 656 (1950). 

(®) M. VipaLe e M. ScREIN: Nuovo Cimento, 8, 10 (1951); M. VIDALE: 
Phys. Rev., 88, 266 (1952). 

(6) H. Newer: Phys. Rev., 83, 649 (1950). 
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mente sulla più o meno completa e corretta descrizione che attualmente noi 
possediamo dei fenomeni elementari. 

Lo schema entro cui si svolge il bilancio deve quindi essere specificato per 
quel tanto che interessa, e noi ci proponiamo di parlarne tra breve. 

Vogliamo prima dire che il miglioramento notevole ottenuto nel presente 
bilancio rispetto alla situazione precedente dipende in massima parte da una 
migliore definizione del secondo punto e cioè della intensità delle varie com- 
ponenti secondarie, e da un migliorato esame dei fenomeni di bassa energia. 
Seguiranno a questo bilancio quelli relativi a latitudini più meridionali, ed è 
sperabile che da una sempre più attendibile definizione di essi si possa com- 
prendere se effettivamente abbiamo una conoscenza adeguata di tutti i pro- 
cessi che si svolgono tra particelle elementari. 

Il bilancio si riferisce al fascio verticale della radiazione cosmica, e la gran- 
dezza che si cerca di definire è la potenza dissipata da un fascio di apertura 
corrispondente ad un angolo solido unitario, in una colonnina di atmosfera 
di base unitaria; anche quando abbreviatamente parleremo di energia dissipata 
dal fascio, i valori sono da intendersi riferiti a MeV-em-?:s-1-sterad. 


2. — Schema fenomenologico. 


Assumeremo che la primaria consista solo di protoni e di nuclei nudi di 
elementi più pesanti dell’idrogeno (rispettivamente circa 85% e 15% del flusso 
totale a questa latitudine). Escluderemo un intervento sensibile di elettroni 
di grande energia (?), come pure di neutroni, od altre particelle neutre che 
abbiano notevole interazione con la materia. 

Assumeremo che l’unica fonte di generazione delle componenti secondarie 
si possa individuare nelle collisioni delle particelle primarie contro i nulcei 
della atmosfera, sia che si tratti di generazione immediata o di generazione 
indiretta attraverso particelle intermedie. In uno di tali urti in generale ver- 
ranno prodotti vari secondari, e precisamente: 


1) Nucleoni di vario tipo, o isolati (protoni e neutroni) o. associati (nuclei 
nudi) capaci di ridare nuove collisioni del tipo di quella in esame. 


2) Nucleoni di bassa energia che verranno progressivamente frenati nei 
successivi attraversamenti nucleari. 


3) Mesoni x positivi e negativi, di vita media breve 2,6-3,0-10-8 s), che 
danno origine per disintegrazione alle particelle » di vita media 2,2-10-* s nel 
processo z+ > u* + v, cui accennammo prima. 


4) Mesoni x neutri, di vita media brevissima (~ 10-14 s), che per disin- 
tegrazione danno origine a due quanti gamma, nel processo 7° > Y+ y. 
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5) Mesoni k, +, ... pesanti, generati in numero sensibile solo da primari 
di energia molto elevata. Questi mesoni non entreranno nel presente bilancio, 
perchè alla latitudine in esame il valor medio della energia dei primari non 
è molto elevato; è presumibile invece che essi abbiano un sensibile peso nei 
bilanci a latitudini più basse. 

In conclusione riteniamo che il grosso della fenomenologia si possa spiegare 
tenendo conto solo dei primi 4 punti. Nel grafico che segue è indicato come 
hanno origine o vengano alimentate nella atmosfera le note componenti meso- 

nica elettrofotonica e nucleonica, cui 

primaria bisogna risalire dai dati sperimentali. 

-—- Una collisione del tipo di quella de- 

i | scritta, in cui si generano anche me- 

soni, si può osservare con le varie 

tecniche in uso, e si possono anche 

individuare entro certi limiti i vari 

prodotti; l’analisi di tali eventi è con- 

È lai corde con quanto da noi assunto. 

In tali collisioni la perdita di energia 

del nucleone che ne è responsabile è 

molto grande, in modo che la compo- 

nente primaria viene degradata ed assorbita in pochissimi urti consecutivi, 

in media non più di 2, per cui l’assorbimento della primaria si svolge essen- 

zialmente nella alta atmosfera. Molto lunga e varia invece è la catena dei 

fenomeni di energia relativamente più modesta in cui le componenti seconi 

darie vengono progressivamente degradate in energia ed assorbite. A parte 

la componente neutrinica, che si suppone non abbia interazione sensibile d- 

nessun tipo con la materia, la componente più lentamente assorbita è quella 

delle particelle u (cosidetta componente mesonica) ed è appunto questa com- 

ponente che trasporta una piccola frazione residua della energia primaria 

(circa il 3%) sottoterra. Altri dettagli sui processi secondari verranno spe- 

cificati in sede più opportuna e cioè a proposito della discussione sulle sin- 
gole componenti. 


nucleonica 
bassa energia 


nucleonica | alta energia 


nucleonica 
ea 


bassa energla 


3. — Dati sperimentali. 


Abbiamo già detto nella introduzione che per condurre il bilancio bisogna 
conoscere le intensità delle varie componenti lungo l’atmosfera, cioè avere 
quelle che si dicono le curve di altezza. Ma, dai dati sperimentali forniti dai 
contatori non si hanno immediatamente le componenti che interessano, ma 
certi loro miscugli, noti universalmente col nome di Totale, Dura, Molle, oltre 
ai neutroni lenti che esamineremo a parte. Le curve di altezza di queste com- 
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ponenti sono rappresentate in fig. 1 e nella tabella che segue sono specificate 
le tracce relative, cioè una quantità che misura il percorso totale delle parti 
celle di ogni componente, definito, per la componente n-sima dalla nota re- 
lazione 


TABELLA I. 


| Componente | Totale | Dura | Molle 
—— ———__k —-- —_— — — — — — — | ——— oa = oo — — 
| | 
Traccia | (per differenza) 
| (g:em-4-s—-sterad™) | 155,5 | 55,5 100,0 


I, (em-2 s-' sterad-') 


0 500 1000 9 fem? 
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Come è noto la componente totale e la componente dura si possono definire 
strumentalmente, come le componenti registrate da un telescopio di contatori 
rispettivamente scoperto oppure coperto con un assegnato spessore di piombo, 
che funziona da assorbitore, e che nel caso nostro precisiamo in 12 em. Il tele- 
scopio scoperto di contatori (una tripla o una quadrupla) è capace di regi- 
strare tutte le particelle cariche del fascio verticale dei raggi cosmici che hanno 
energia sufficiente per passare le pareti dei contatori stessi (complessivamente 
circa 5 g/cm? di ottone); il telescopio coperto coi 12 cm di piombo registra 
invece solo le particelle che o passano direttamente tale spessore, o generano 
nell’assorbitore almeno un secondario capace di attraversarlo. La componente 
molle è invece definita per differenza tra la totale e la dura, e rappresenta 
quelle particelle che vengono eliminate, assieme ai loro secondari ionizzanti, 
dall’assorbitore. Si comprende immediatamente che dura e molle, definite 
strumentalmente sono un miscuglio, variabile con la quota delle effettive com- 
ponenti che ci interessano; così la componente dura consisterà di mesoni & 
veloci e di protoni veloci, mentre la componente molle sarà formata da elet- 
troni, mesoni uv lenti, protoni lenti e protoni veloci che hanno interagito nel- 
l’assorbitore senza dare secondari ionizzanti rilevabili dal telescopio, ecc. Il 
primo compito consiste quindi nel separare le componenti effettive in seno 
alla dura e alla molle e a ricostruirle in tutta la atmosfera. Il compito è abba- 
stanza agevole al livello del mare e sotto terra per il fatto che le uniche parti- 
celle cariche superstiti in quantità sensibile sono elettroni e mesoni u, dato 
che i protoni al livello del mare costituiscono (7) solo 11,5% del flusso verticale. 

Al limite della atmosfera si presenta anche una situazione molto semplifi- 
cata per il fatto che si può assumere che la totale sia costituita dai nucleoni 
primari, e la dura da quella frazione di tali nucleoni che ha passato Vassor- 
bitore di piombo; nella alta e media atmosfera la separazione è più compli- 
cata e non manca di qualche incertezza. Ad operare tale separazione si impie- 
gano pertanto delle informazioni complementari quali la crescita con l'altezza 
dei mesoni lenti, e l'assorbimento della componente nucleonica nella atmosfera. 
La crescita dei mesoni lenti nella atmosfera, in funzione della quota si può 
stabilire per via diretta, con un dispositivo a coincidenze ritardate (8) sfrut- 
tando la proprietà delle particelle u di disintegrare in elettrone con una vita 
media nota. La conoscenza della curva di altezza dei mesoni lenti (fig. 2) unita 
alla conoscenza dello spettro differenziale al livello del mare dei detti mesoni, 
ottenibile anch’esso per via diretta, è importante perchè permette di definire 
la funzione di produzione per i mesoni p e pertanto anche la legge di assor- 


(7). M. G. Mytror e J. G. Witson: Proc. Roy. Soc., A 64, 404 (1951). 
(8) a) B. Rossi, M. Sanps e R. F. SARD: Phys. Rev., 72, 120 (1947); b) M. Sanps: 
Phys. Rev., 77, 180 (1950); c) M. Conversi: Phys. Rev., 79, 749 (1950). 
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bimento della nucleonica che li produce. La legge di assorbimento che così si 
deduce è in accordo con quanto si può dedurre dalle misure dirette fatte in 
camera di Wilson o con altre tecniche. Se si conosce la legge di assorbimento 


OR 


i (gr-? s-? sterad-") 


della componente nucleonica è possibile separare in seno alla componente dura 
i mesoni u veloci dai protoni. 

Considerando che i protoni veloci vengono assorbiti con un cammino di 
assorbimento di 120 g/em2, ed un cammino d’urto di 85 g/em? (*), e che nella 
atmosfera la componente nucleonica di media ed alta energia di carattere 
secondario, è composta all’incirca per metà di neutroni, si vede che una ragio- 
nevole legge per l'assorbimento della nucleonica primaria può essere scritta 


(9) E. CLEMENTEL e G. PUPPI: Nuovo Cimento, 8, 936 (1951); vedi anche: E. Coc- 
cont, V. G. ToncrorgI e M. WIDGOFF: Phys. Rev., 79, 768 (1950); G. BERNARDINI, 
G. CORTINI e A: MANFREDINI: Nuovo Cimento, 6, 456 (1949). 
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al modo seguente: 


J = J \exp [— 1/85] + E (exp [— t/120]— exp [ 185) 5 


dove t è la profondità nella atmosfera misurata in g/cm? a partire dal limite 
della atmosfera stessa (se vogliamo, è la pressione, nella unità predetta, misu- 
rata in quel punto) e J, è il valore del flusso primario tagliato dall’assorbitore. 
Ovviamente la legge scritta rappresenta un andamento globale della nucleo- 
nica, e non ha la pretesa di essere una descrizione dettagliata del vero stato 
di cose, tra l’altro perchè v’é nella primaria il famoso 15% di nuclei nudi, il 
cui assorbimento certamente è differente da quello dei protoni; d’altra parte 
il modo di interagire dei detti nuclei con la materia non è ancora così ben noto 
da permettere una analisi differenziata. Probabilmente i cammini ora citati 
sono una opportuna media ponderata tra un cammino di assorbimento proprio 
dei protoni ed un cammino per particelle « ed altri nuclei. 

La sottrazione dalla dura dei protoni veloci permette di ricavare i » veloci 
e precisamente i pu di range superiore a 12 cm di piombo. 

Il computo dei mesoni u tagliati, di energia quindi compresa tra 0 e 200 MeV, 
sì conduce in base alla accennata legge sulla crescita dei mesoni a fine percorso, 
stabilita con un dispositivo a coincidenze ritardate e alla conoscenza della forma. 
dello spettro differenziale dei u in tale rango di energie. Questa forma varia 
molto lentamente con la quota e si può dire che l'andamento con l'altezza. 
dei mesoni lenti nel detto intervallo energetico è ben rappresentato dall’anda- 
mento dei mesoni che stanno al mezzo della banda e cioè quelli da 100 MeV. 

Il computo del contributo di questi mesoni lenti e la loro addizione ai mesoni 
veloci permette in definitiva di ricostruire la componente mesonica totale. 

Per ottenere la componente elettronica bisogna ora tentare di computare: 
la curva di assorbimento dei protoni totali, e toglierli, assieme alla mesonica 
totale alla componente totale. Per i protoni totali abbiamo il valore al limite 
della atmosfera e, dagli spettri noti alle varie quote (?) possiamo avere alcuni 
valori nella media e bassa atmosfera; abbiamo tracciato l’andamento globale, 
tenendo presente questi dati da una parte, e l'andamento dei protoni veloci 
dall’altra. Alla fine di questo lavoro di analisi possediamo le curve di altezza 
per le varie componenti fisiche, mesoni, elettroni, protoni, su cui possiamo comin- 
ciare a costruire il bilancio, riportate in fig. 3. Nella tabella che segue, infine: 
sono specificate le varie tracce. 


TABELLA II. 


| | 
Totale | | 
| 


7 ; | Mesoni Mesoni | Protoni Protoni Elet- 
DOM PONenbe,] veloci totali | veloci | totali | troni 


| 
| 
| 
| | 
| Tracce 34,2 34,7 20,0 30,0 86,6 155,5 | g-cm-4-s-1-sterad-1 
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I, (cm-? s~? sterad-1) 


Ome 


(o) 500 g /cm? 1000 


Fig. 3. 


4. — Componente mesonica. 


Come abbiamo gid accennato nella introduzione, negli eventi nucleari ven- 
gono generati i mesoni x carichi, identici a quelli che si producono quotidia- 
namente in laboratorio con le grandi macchine, ma essi, stante la breve vita 
media, non costituiscono una componente sensibile, almeno fino a che possie- 
dono energie moderate, e pertanto ogni informazione su questa componente 
dovremo dedurla dall’analisi della componente u. Solo per energie estrema- 
mente elevate, dell’ordine delle decine di GeV, la vita media di questi mesoni x 
diventa abbastanza lunga da permettere dei notevoli tragitti. D'altra parte 
i mesoni x di tale energia sono pochi, e stante la loro forte interazione nucleare 
non si potrebbero distinguere coi mezzi usuali dai protoni, così che possiamo 
in effetti concludere che la componente x non interviene mai come compo- 
nente a sè. 

Le particelle u invece, stante la vita media relativamente lunga e la presso- 
chè nulla interazione nucleare, danno luogo nella atmosfera sostanzialmente a 
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due processi e precisamente: perdite di energia per ionizzazione e disintegra- 
zione spontanea in elettrone e neutrini. Dalla energia dissipata in questi due 
processi, più l’energia residua al livello del mare, dovremo quindi risalire alla 
energia globale della componente u e da questa, tenendo conto del passaggio 
T — wu, alla energia globale convogliata nella componente x. L'energia residua 
al livello del mare si calcola facilmente a partire dalla conoscenza dello spettro 
differenziale; detto infatti f(p, t)) dp lo spettro differenziale in momento e detta 
e(p) l'energia totale che compete ad un momento p, la si può scrivere nella 
forma: 


(ce) 


| Mita. 


. 


0 


Il computo dell’integrale fornisce il valore di 37 MeV-cm~*s~ sterad~! per 
l’energia residua. Le perdite di energia per ionizzazione lungo Vatmosfera si 
ottengono poi dalla espressione 


to (cel 


[ae | paio var, 


Ome 

dove £, rappresenta la perdita di energia, per unità di percorso, di un mesone vu. 
di momento p. Per il computo dell’integrale è necessario conoscere gli spettri 
differenziali in momento dei y lungo tutta l'atmosfera. Si può però procedere 
anche, seguendo Rossi, distinguendo le perdite dei mesoni veloci, che si pos- 
sono ritenere costanti, da quelle dei mesoni lenti, variabili col momento. Per 
i mesoni veloci si ha immediatamente come espressione delle perdite il pro- 
dotto della perdita costante per la traccia, e cioè precisamente 


(LS, , 


mentre per i mesoni lenti, detta £, la energia massima del gruppo e R,, il 
«corrispondente percorso massimo e /,,, la intensità verticale del gruppo lento, 
si può scrivere immediatamente, come espressione delle perdite globali 


ty 

Ey 

— I, it) dt. 
fz i(t) € 


0 


Il computo numerico delle perdite per ionizzazione di tutti i mesoni ‘con- 
duce al valore di 


77 MeV cm= gs-}-sterad-. 
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Resta ancora da specificare un modesto contributo dovuto a quei mesoni 
che vengono completamente frenati, e che pertanto, per lo più, da fermi disin- 
tegrano al modo noto. Il contributo in parola si può precisare in 


to 


[mR f(0, t) dt, 


0) 


© anche approssimativamente, secondo Rossi, tramite la espressione seguente 


che rappresenta il prodotto della energia intrinseca della particella u, per un 
valor medio dello spettro differenziale in range dei mesoni del gruppo lento. 
Risulta in entrambi i casi che il contributo in parola vale circa 5 MeV-em™- 
-g-1-sterad™. 

Veniamo ora al calcolo delle perdite per disintegrazione in volo, che costi- 
tuiranno il contributo più vistoso. L’energia che sparisce dal fascio dei pu, ‘a 
seguito della disintegrazione spontanea di queste particelle si può scrivere 
immediatamente nella forma 


pt IV ZIONE 


essendo b = H,m,e°/ct,, con H, altezza della atmosfera omogenea, m,, Massa 
della particella u e 7, vita media a riposo; al solito t rappresenta il livello 

del mare. 
L’espressione ora scritta si giustifica, come pure tutte le precedenti, attra- 
verso la equazione unidimensionale di diffusione per i u che si può. scrivere 
of 


9 vane 


Moltiplicando ambo i membri per e e integrando su tutti gli impulsi e su 
tutta l'atmosfera si ottiene la relazione seguente, che riproduce per l’appunto 
il bilancio sulla componente in esame: 


ty (ce) (ce) to (cel 


far[ es, Oo = | ef (py to) dp + fa [si t)dp + 


0 0 0 0 0 


to to 


È. | m,¢2B(O)f(0, 1) dt + | at | = 


0 ò 
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in cui il primo membro rappresenta la energia totale spesa nella produzione 
di mesoni in tutta la atmosfera. Il secondo membro invece mostra come questa 
energia si ripartisce e precisamente in energia residua al livello del mare, 
energia spesa in ionizzazione, energia dissipata in disintegrazione sia da fermo 
che in volo. 

Introdotta una funzione ausiliaria x(t) la espressione delle perdite per di- 
sintegrazione si può riscrivere al modo seguente: 


to 


fron + ma, 
con 


Her A 5 
a(t) = | = es ap / | f(p, t) dp . 


Questa funzione, che dipende sostanzialmente dalla proporzione tra mesoni 
di bassa e di alta energia, cioè dalla forma dello spettro differenziale sulle basse 
energie è stata computata in base agli spettri noti (*) per le quote 1033, 680, 
300 g/cm?. Al disopra di 300 g/em? non abbiamo spettri e pertanto abbiamo 
calcolato la forma di f(p, 0) al limite della atmosfera ricorrendo alla equazione 
di diffusione considerando nota (vedi più avanti) la funzione sorgente S(p, t) 
e precisamente della forma: 


il 
Sip) = G,(p) TP xia Fale 


Poichè f(p, 0) = 0, si può sviluppare la funzione f e la funzione sorgente 
in serie di Mac Laurin 


t t? 
ip, t) = h (p) iy + f.(P) DIE + ‘034 

= al ij t2 

S,(p, t) = 5 Ju(P) ea ail ae CS 


per cui si ottiene che la equazione di partenza si spezza nella serie delle equa- 
zioni: 


0 
is [Ay gu 
ge ing tare ee 


p 1-35 
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Per t > 0 pertanto la forma della funzione f sarà data dalla prima equazione. 
| La funzione «(t) è riprodotta in fig. 4 assieme ad altre funzioni ausiliarie 


1 L 


DA ie: 
DI 
aN ie 


1, + ale)) (cm 2 s-? sterad*) 
I,(em-2.s~? sterad~") 


6 
t 


rly 
i LE ire atti) 
10-3 
0 500 1000 9 /em? 


Fig. 4. 


per il calcolo della disintegrazione di u. Un calcolo numerico ha fornito come 
valore delle perdite di energia per decadimento dei u nella atmosfera 


264 MeV-cm-2 s-! sterad= . 


Di questa energia 1/3 si ritrova in componente elettrofotonica e 2/3 invece 
sparisce definitivamente dal bilancio perchè trasportata dai neutrini, che pra- 


10 


10-2 
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ticamente non interagiscono con la materia. Se sommiamo le varie fonti di 
dissipazione di energia dei u giungiamo a 383 MeV-cm-? s- sterad-!. Ricor- 
dando che i p. derivano dai x e che la energia del x è circa 1,27 volte (1°) quella 
del u, si arriva in definitiva a risalire alla energia che va alla componente 7 
carica, specificata in tabella. 


TABELLA 

| Perdite in aria dei yp per ionizzazione . . . . . 77 MeV cm? 81 sterad-! 
| Perdite sotto terra dei y . . . . . sica SIN 37 

Perdite per disintegrazione in aria Hebe Ur oo 5 2.00) 

Perdite:totall der i 0.20 5G ang aan e: 
| Perdite in neutrini nel passaggio 7 + u | 105 
} | 
ar 

Energia totale che va in 7 carichi | 488 


Come altra volta accennato (?) la energia totale che va in una componente 
può essere stimata anche a partire dallo spettro e dalla equazione di diffusione. 

Peri x e uw carichi possiamo scrivere la coppia di equazioni seguenti, consi- 
derando come variabile la energia totale ed assumendo un cammino di assor- 
bimento per interazione nucleare per i x pari a quello della nucleonica 


od a; 1 
n (Ba) È ae a x - TOP [— t/L]gn(E) , 

(peer 
du lo) b 


B 
= A = Ì 7 Ù 5 
algo e yk 


È facile vedere che per energie sufficientemente elevate, la soluzione per 
i, assunto g,(£) = yA/H’** si può scrivere: 


A SO 


= BL ky re? 


Confrontando questa espressione con lo spettro differenziale dei u al livello 
del mare, quale risulta e dai dati sperimentali diretti e dall’assorbimento sotto 
terra (1), si ha la possibilità di fissare un valore assoluto per lo spettro di 
produzione dei 7 e precisamente: 


i 


0,23 
paio cm7?-s7!- sterad— 


(10°) W. BrrnBAUM, F. M. Smita e W. H. BARKAS: Phys. Rev., 83, 895 (1951). 
(11) E. P. GEORGE: oa in Cosmic Ray Physics, vol. I, Cap. VII (North Holland 
Publishing Co., Amsterdam, 1951). 
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Questo spettro andra certamente modificato verso le basse energie, ma è 
arduo effettuare questa operazione attraverso una equazione di diffusione. 
Si pud invece continuare lo spettro verso le basse energie, tenendo conto della 
forma stabilita attraverso le misure dirette (*1%13). Abbiamo eseguito arbi- 
trariamente il raccordo tra i due andamenti a 5 GeV. Il risultato è riportato 
in fig. 5 assieme allo spettro differenziale in energia dei u al livello del mare. 
Dato il limitato rango energetico in cui siamo interessati, ripetiamo che non 
riteniamo che i risultati debbano essere sensibilmente modificati per tener 
conto di quella parte delle particelle u che potrebbero derivare da mesoni più 
pesanti dei x. Effetti più fini invece, potrebbero risentirne, analizzati su parti- 
celle p di elevatissima energia, ma questo è fuori dello scopo della presente 
ricerca e verrà discusso in altra sede. 

Se ora calcoliamo, a partire dallo spettro di produzione la energia globale 
convogliata in mesoni x carichi arriviamo al valore di 515 MeV-cm-’:s?: 
‘sterad-! in ottimo accordo con il calcolato attraverso le perdite. 

I risultati ottenuti abbisognano però di una importante correzione; infatti 
il bilancio da noi effettuato sulla dura, riguarda il fascio verticale, e sarebbe 
un indizio fedele della energia dissipata dalla primaria nello stesso angolo 
solido, se i secondari fossero collimati con i primari che li producono. Potremmo 
in questo caso scrivere che la secondaria prodotta per unità di angolo solido 
alle varie quote è grosso modo proporzionale alla intensità verticale della radia- 
zione generatrice e cioè 


Doi = le OX t/L] , 


col 
supponendo che la deformazione dello spettro primario mano a mano che si 
penetra nella atmosfera non alteri sensibilmente questo andamento. L'altro 
caso limite si avrebbe invece per una distribuzione isotropa della radiazione 
secondaria, e in questo caso si può scrivere: 


B { : 1 
Lote = 5 Jo exp 41} 1-7 exp EHF | exp (— vit 
4/L . 


27 


4rrlisote = B 7(0) dQ. 


Si vede subito che la perdita nella secondaria isotropa, relativamente alla 
componente verticale è di 1/2 al limite della atmosfera, assumendo la radia- 


(2) U: CAMERINI, P. H. FowLER, W. 0. Lock e H. MUIRHEAD: Phil. Mag., 41, 


415 (1950). 
(18) R. L. CooL e O. Piccioni: Phys. Rev., 87, 531 (1952). 
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zione primaria generatrice isotropa all’ingresso ed assorbita direzionalmente 
con legge esponenziale; la perdita però raggiunge a mezza atmosfera già il 90%. 
Questo significa che, pur tenendo conto del compenso dovuto alla primaria 
inclinata, solo una frazione del numero di mesoni prodotti isotropicamente 
dalla primaria compresa in un angolo solido unitario in verticale resta nel 
medesimo angolo solido. I due casi citati sono due casi estremi, e naturalmente 
la situazione reale è molto più complessa e varia continuamente a seconda 
della energia dei primari e dei secondari che si considerano. Possiamo tuttavia 
tentare di ragionare al modo seguente: supponiamo che i mesoni 7 Si possano 
suddividere, alla generazione, in due gruppi e precisamente: un gruppo iso- 
tropo comprendente mesoni di energia inferiore ad un valore £. ed un gruppo 
collimato di energia superiore a questo limite. In questo caso la correzione 
angolare si applicherebbe al solo gruppo isotropo e si calcolerebbe facilmente 
in base alle formole precedenti. Si può giungere così immediatamente al risul- 
tato che solo il 30% circa di tutti i mesoni prodotti isotropicamente entra 
nel fascio verticale di apertura unitaria. Per passare alla correzione sulla energia 
globale convogliata in mesoni x carichi, bisogna fissare il valore di E,. Si può 
per questo ricorrere alle esperienze di generazione di mesoni x di bassa energia 
in materiale condensato. Dai risultati (12) in lastre (!*) e in carbone si ricava 
un valore di E, intorno ai 200 MeV, ma è molto probabile che tale valore nel 
nostro caso rappresenti solo un limite superiore; infatti tali mesoni isotropi 
provengono sostanzialmente nel caso delle lastre fotografiche dallo scattering 
entro il nucleo stesso in cui vengono prodotti, nel caso della esperienza di pro- 
duzione in carbone, oltre che dallo scattering entro il nucleo stesso in cui ven- 
gono prodotti, anche dallo scattering subito nell’attraversamento di nuclei 
successivi (15). In entrambi i casi c'è da attendersi una maggiore quantità di 7 
diffusi che non nel caso della atmosfera, i cui nuclei sono abbastanza piccoli 
e non vi è possibilità di interazioni successive a cagione del rapido decadi- 
mento. Vi è però ancora un fattore che è a favore del caso dell’aria, ed è lo 
sparpagliamento angolare che si introduce nella disintegrazione 7 > yu quando 
l’energia del x è relativamente bassa. Se riteniamo comunque che E, debba 
essere compreso tra 100 e 200 MeV, la correzione da apportare sulla energia 
risulta compresa tra il 15 e il 30%. Come si vede si tratta di un contributo 
che in un futuro bilancio, quando molti dati sperimentali avranno raggiunto 
una consistenza definitiva, bisognerà cercare di definire meglio di quanto non 
si faccia ora. Bisogna anche dire che le percentuali citate per la correzione 
angolare vanno rivedute per tener conto dell’albedo. Dai risultati finora noti 


(14) J. C. BARTON, E. P. GEORGE e A.C. Jason: Proc. Phys. Soc. Am., 64, 175 


(1951). 
(15) E. CLEMENTEL e G. PUPPI: Nuovo Cimento, 10, 197 (1953). 
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non è possibile precisare quale frazione di questa radiazione noi raccogliamo 
dal basso in alto nei nostri dispositivi, che non sono generalmente protetti per 
l’albedo; comunque il tenerne conto porta a ridurre la correzione angolare, e 
nel caso limite che le misure registrassero integralmente la radiazione all’in- 
dietro le cifre delle percentuali indicate andrebbero ad un dipresso dimezzate. 

Tentiamo prudenzialmente di fissare in 62 MeV-cm-? s- sterad- la cor- 
rezione, e ricordiamo che Rossi, in base a considerazioni di carattere teorico 
sulla angolazione del processo di generazione giunge a stimare correzioni per- 
centuali della stessa entità. Un bilancio condotto sul flusso omnidirezionale, 
a parte la questione dell’albedo, sarebbe esente da questa correzione; pur- 
troppo il quadro sperimentale attuale non consiglia ancora un bilancio di questo 
tipo; esso sarà però altamente istruttivo comunque, anche per le latitudini 
più meridionali. 

Apportando la correzione convenuta otteniamo come valore finale della 
energia convogliata in mesoni x carichi: 


550 MeV-cm- s- sterad+. 


Il valore ottenuto è, come già detto all’inizio, sensibilmente superiore a 
quello dei precedenti bilanci. È interessante calcolare il numero medio di me- 
soni prodotti nella colonnina di atmosfera in oggetto; questo numero risulta 
dipendere in maniera piuttosto sensibile dalla entità del gruppo isotropo di 
bassa energia, ed il suo valore numerico oscilla intorno a 1 cm-? s+ sterad—. 
Con un valore del flusso primario intorno a 0,20 cm-? s— sterad-! ($ 8) questo 
significa che vengono prodotti in media 5 mesoni 7 carichi per nucleone o 
nucleo primario, e che l’energia media dei mesoni prodotti è di poco superiore 
ad 1 GeV; questo risultato è ovviamente e non indipendentemente, in accordo 
con i risultati diretti sullo spettro di produzione. Poichè dal confronto tra 
cammino d’urto e cammino di assorbimento per la primaria risulta che il 
numero medio di urti anelastici contro i nuclei d’aria sufficienti per assorbirla 
è un poco inferiore a 2 ne viene che in media in ogni collisione anelastica di 
un nucleone energico contro un nucleo di aria si ha la produzione di circa 
3 mesoni 7 carichi. Tenendo conto di un limitato processo a cascata entro i 
nuclei di aria, non ci sembra che i valori riportati sieno in aleun modo critici 
per una forma particolare di teoria di produzione; e ciò è anche ragionevole 
se si pensa che il nostro panorama è di non elevatissima energia. 


5. — Componente elettrofotonica. 


Dal quadro fenomenologico tracciato all’inizio risulta chiaro che la com- 
ponente elettrofotonica nella atmosfera deriva da tre fonti principali e cioè: 
1) Dai processi di knock-on e di bremsstrahlung delle particelle u e 

dei mesoni 7. 
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2) Dalla disintegrazione in volo delle particelle » in elettrone e neutrini. 
3) Dalla disintegrazione dei mesoni x neutri in 2 quanti y. 

Da tutti e tre questi processi infatti hanno origine o elettroni o fotoni 
energici i quali, moltiplicandosi in cascata, formano la componente elettro- 
fotonica. Si comprende facilmente come l'entità della componente può essere 
valutata a partire dalla conoscenza degli spettri e delle curve di altezza degli 
elettroni e dei fotoni; purtroppo il limitatissimo panorama sperimentale attual- 
mente esistente sui fotoni ci costringe a fondare il presente bilancio, come i 
precedenti, solo sulla conoscenza della componente elettronica, ed anche questa 
determinata per differenza tra totale e dura. 

In altre parole, quello che conosciamo approssimativamente è solo la curva 
di altezza di tutti gli elettroni di energia superiore a circa 10 MeV, determi- 
nata dal taglio dei contatori. Il legame tra la traccia globale che dalla detta 
curva di altezza si ricava, e la energia totale dei primari che hanno prodotto 
la componente elettrofotonica è fornito dalla teoria della cascata (!); Rosst, 
ad esempio, in approssimazione unidimensionale risale dalla detta traccia degli 
elettroni secondari alla energia globale dei primari moltiplicando il valore 
numerico della traccia per 3,26 MeV/g:em-*. Per passare a delle valutazioni 
numeriche bisogna però cercare di separare i vari contributi, e si può iniziare 
dal primo, che si valuta facilmente ammettendo una condizione di equilibrio 
lungo tutta la atmosfera tra la componente mesonica e la componente elettro- 
fotonica secondaria. Osserviamo che il processo più importante è senz'altro 
quello di knock-on delle particelle u perchè il processo di bremsstrahlung è 
subìto solo da mesoni di energia estremamente elevata, che sono molto pochi, 
e d’altra parte il processo di knock-on dei mesoni x porta un contributo tra- 
scurabile, stante il fatto che, come in precedenza accennato, i detti mesoni 
non formano una componente sensibile. Dobbiamo ancora dire che la supposta 
condizione di equilibrio tra particelle u e componente elettrofotonica secon- 
daria, generata nel processo di knock-on, è realizzata con buona approssima- 
zione solo nella media e bassa atmosfera, ma non nella altissima atmosfera, 
dove la intensità della componente » varia molto rapidamente con la quota; 
il valore numerico del contributo in esame è però così modesto, agli effetti 
del bilancio, che non ha senso tentare di precisarlo ulteriormente. 

Pertanto, assumendo che al livello del mare la elettronica secondaria dei u 
per knock-on ammonti al 7% circa della componente p (') ed integrando su 
tutta la atmosfera in base alla detta condizione di equilibrio risulta una traccia 
globale di 2,6 g-em-4-s-1-sterad-1. Resta così come traccia totale relativa 
ai processi 2) e 3) circa 84 g-cm~‘s-! sterad-, che corrisponde ad una energia 
per i primari di 

275 MeV-cm- s- sterad™ . 


Questa cifra va corretta per varie ragioni, ma sostanzialmente per tener 
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conto dello sparpagliamento laterale degli elettroni di bassa energia che sono 
predominanti. Un calcolo di tale correzione appare molto complicato, e noi 
ci limiteremo quindi a valutarne la entità dal confronto, al livello del mare 
tra le distribuzioni zenitali della componente dura e della componente elettro- 
nica. Risulta, dai dati sperimentali (1%) che tale dipendenza zenitale è meno 
pronunciata per la componente elettronica che non per la componente dura; 
poichè al livello del mare le due componenti sono certamente in equilibrio, 
la differenza può essere imputata allo sparpagliamento laterale degli elettroni 
di bassa energia. Dal confronto quantitativo tra le due distribuzioni si può 
tentare di fissare nel 25% la correzione in esame e giungere così ad una energia 
globale per i primari dei processi 2) e 3) di 


344 MeV-cem-? s71 sterad- . 


Il contributo 2), della disintegrazione dei u è anche accessibile al com- 
puto diretto. Si può infatti, in base alla conoscenza degli spettri differenziali 
dei u lungo l'atmosfera, alla vita media e al modo di disintegrazione, risalire 
agli spettri degli elettroni che provengono dalla disintegrazione; in approssi- 
mazione unidimensionale e riferendoci agli spettri medi, cioè assumendo che 
la energia dell’elettrone di disintegrazione sia la sua energia media, si può 
infatti scrivere 


fle, t)de = (3, f,(Be, 1) de) 


come spettro differenziale degli elettroni di decadimento, essendo f, lo spettro 
differenziale dei mesoni. Un calcolo in tal senso è stato fatto recentemente (17), 
e permette di stabilire, alle varie quote la percentuale di elettroni che derivano 
dal decadimento dei u; per i fini del bilancio questo non ha interesse e basta 
sapere, sulla energia globale della componente elettrofotonica, quale frazione 
sia da attribuirsi al decadimento dei u. Basta pertanto ricordare che solo 1/3, 
in media, della energia che va in disintegrazione dei si ritrova in compo- 
nente elettrofotonica; e ricavare tale valore dal bilancio della componente me- 
sonica. Resta così come contributo della terza delle fonti citate all’inizio 


256 MeV-cem-? s-1 sterad-! . 


Questa energia rappresenta anche, nel nostro schema fenomenologico, la 
energia convogliata in mesoni 7 neutri, e contiene una correzione angolare 


(16) W. L. KRAUSHAAR: Phys. Rev., 76, 1045 (1949). 
(17) R. STROFFOLINI: Nuovo Cimento 10, 300 (1953); D. CaporIN: Nuovo (i- 
mento 10, 364 (1953). 
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per il processo di generazione dello stesso ordine di grandezza di quella discussa 
per i x carichi. 

Dal confronto con il bilancio della componente mesonica possiamo vedere 
che la energia globale dissipata dalla primaria per la produzione di mesoni x 
neutri è metà della energia dissipata in produzione di x carichi. Se, come i 
risultati sperimentali (12) tenderebbero ad indicare, gli spettri di produzione 
dei mesoni x carichi e neutri sono simili, il confronto si può ulteriormente pre- 
cisare dicendo che il rapporto di produzione mesoni neutri/mesoni carichi è 1/2. 
A questo proposito è opportuno notare che tutte le determinazioni sperimen- 
tali dirette di questo rapporto, effettuate analizzando le particelle emesse negli 
eventi nucleari, in camera di Wilson, o con lastre (?:15) confermano tale valore. 
Non ci si può esimere dal far notare che, nella teoria mesonica delle forze 
nucleari, in cui mesoni neutri e mesoni carichi corrispondono rispettivamente 
a campo scarico e campo carico, un rapporto di 1/2 quale quello trovato rap- 
presenta appunto il dosaggio caratteristico di quella tale versione, tra le infinite 
possibili, che va sotto il nome di teoria simmetrica. Questa particolare ver- 
sione è risultata la più adeguata a descrivere i fenomeni nucleari di bassa 
energia, quali il legame nucleare e lo scattering elastico nucleone-nucleone, 
ma nulla autorizzava a pensare che si dovesse trovare una indicazione a suo 
favore anche nei fenomeni di energia molto elevata connessi con collisioni radia- 
tive tra nucleoni. Con una immagine tolta dal metodo semiclassico di Weiz- 
sicker e Williams, potremmo dire che i mesoni reali, emessi nelle collisioni 
radiative di elevata energia, conservano lo stesso rapporto scarico/carico dei 
mesoni virtuali responsabili, nei termini della teoria di campo, dei fenomeni 
non radiativi, quali, per l'appunto il legame nucleare e lo scattering elastico 
nucleone-nucleone. 


6. — Componente nucleonica. 


Sotto questa voce raccoglieremo il bilancio della energia dissipata dalla 
componente nucleonica nella atmosfera in varie forme e precisamente: ioniz- 
zazione prodotta dai protoni e da altri nuclei, energia di legame dei nuclei 
d’aria sedi di disintegrazioni nucleari, energia cinetica di particelle di bassa 
energia evaporate in tali disintegrazioni, emissione di quanti y di bassa energia; 
esclusa la produzione dei mesoni x carichi e scarichi. Cominciamo dalle per- 
dite di energia per ionizzazione della componente nucleonica carica. Un com- 


puto esatto di tale contributo non è ancora possibile poichè conosciamo lo 
spettro energetico della componente protonica solo approssimativamente all’in- 


(38) G. Sapvint e Y. Kim: Phys. Rev.. 85, 921 (1952); J. TinLor e B. GREGORY: 
Phys. Rev., 81, 667 (1951); vedi anche (1°). 
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gresso della atmosfera e nella bassa atmosfera (7), ma ne ignoriamo comple- 
tamente la forma nella altissima atmosfera, cioè nella zona della prima o delle 
primissime collisioni, dove certamente le cose si complicano. 

Per rendersene conto, basta pensare che lo spettro della nucleonica in in- 
gresso nella atmosfera, è uno spettro tagliato verso le basse energie, e preci- 
samente, alla latitudine in esame, ad una energia di circa 2 Gev, per effetto 
del campo magnetico terrestre. Questa regione dello spettro, al disotto del 
taglio, vuota dunque nello spettro primario, viene progressivamente invasa, 
mano a mano che si scende nella atmosfera, dai protoni secondari prodotti 
nelle collisioni contro i nuclei d’aria dalla componente di energia più elevata, 
fino a che si raggiunge una specie di condizione di equilibrio, o per lo meno una 
condizione lentamente variabile con la profondità, caratteristica della media 
e bassa atmosfera, regolata e dall’assorbimento globale della componente 
nucleonica e dal gioco tra le perdite di energia per ionizzazione e il fenomeno 
dello scambio di carica tra protoni e neutroni (°) negli attraversamenti nucleari. 
Per comprendere la vistosità dell’effetto di transizione ora descritto in questa 
parte bassa dello spettro, ricordiamo che il 65% circa dei protoni presenti 
in tutto lo spettro nella bassa atmosfera appartiene a questa regione, che è 
vuota nello spettro primario. Poichè conosciamo solo l'andamento delle per- 
dite all’ingresso della atmosfera e nella media e bassa atmosfera, dovremo fare 
una ragionevole interpolazione. Le perdite di ionizzazione specifiche nello 
spettro primario si calcolano assumendo la radiazione costituita per 185% di 
protoni e per il rimanente 15%, di particelle «, ed omettendo per semplicità 
di specificare il contributo dei nuclei pesanti. Assumendo come valore del 
flusso globale della primaria (§ 8) 0,20 particelle cm-? s-! sterad-* e perdite di 
energia per ionizzazione costanti, dato il rango energetico considerato, otte- 
niamo come perdita specifics 


0,60 MeV-g— s7 sterad . 


Al livello del mare, dalla conoscenza dello spettro differenziale si pud anche 
valutare la perdita specifica in 0,44 10-? MeV-g-! s~ sterad-. 

Nella bassa e media atmosfera lo spettro dei protoni non varia sensibil- 
mente, come forma, da quello del livello del mare e la crescita globale si può 
ritenere sia rappresentabile con una esponenziale regolata da un cammino di 
assorbimento di 145 g/cm?, per cui possiamo computare il valore delle perdite 
specifiche da livello del mare fino a circa mezza atmosfera. L’interpolazione 
tra la media atmosfera e il limite della atmosfera è stato effettuato tenendo 
presente la crescita dei protoni totali. Risulta, a conti fatti, con una traccia 
di 30 g-em7 s- sterad-! una energia dissipata per ionizzazione di circa 


80 MeV-cm-? s~ sterad— . 


Resta ora da discutere il problema delle perdite di energia della componente 
nucleonica negli eventi nucleari di bassa energia, di cui ci si pud fare la seguente 
immagine. Un nucleone che attraversi un nucleo d’aria e non vi generi mesoni, 
sia che vi resti bloccato, sia che ne.riemerga, cede in generale una frazione 
della sua energia, per collisioni elastiche, contro i singoli nucleoni del nucleo. 
Tl nucleo così eccitato o si frattura completamente, o emette una o più parti- 
celle di bassa energia e quanti y trasformandosi in un nuovo nucleo. In un 
bilancio come il presente, fondato sui dati dei contatori i quanti y sono da 
considerarsi completamente perduti perchè, stante la loro bassa energia danno 
origine ad elettroni che sono al disotto del taglio dei contatori, e lo stesso vale 
per protoni e particelle x. Bisogna quindi introdurre dei dati complementari, 
ed i più completi riguardano l'andamento con la quota delle disintegrazioni 
nucleari e dei neutroni lenti. 

La frequenza delle interazioni nucleari con la quota offre effettivamente 
una possibilità in tal senso, purchè si faccia una ragionevole ipotesi sulla 
energia media dissipata in ognuna di dette interazioni nucleari per la produ- 
zione di secondari di bassa energia e di quanti y. Questo modo di attacco del 
problema, che a prima vista potrebbe apparire piuttosto semplice è in realtà, 
almeno fino ad oggi ancora molto incerto. Infatti i dati dovrebbero essere 
dedotti o dalle lastre fotografiche o dalle camere di ionizzazione o dalle camere 
di Wilson, ma tutti si riferiscono ad eventi ionizzanti avvenuti in mezzi con- 
densati differenti dall’aria; la conversione dei dati al caso dell’aria e più che 
tutto il computo della frazione di eventi non registrati perchè in un modo 0 
nell’altro discriminati dal dispositivo, non permette che delle stime di ordine 
di grandezza. 

Riferendoci infatti ai dati desunti dalle lastre fotografiche, che per questa 
necessità sono i più completi, dall'andamento delle disintegrazioni nucleari con 
la quota si può stimare (*) che nella detta colonnina di atmosfera di 1 em? 
di base vengano prodotte circa 2 stelle al secondo. Tali stelle sono, per i criteri 
di discriminazione usati nella esplorazione delle lastre nucleari, da cui tale 
numero è dedotto, con almeno tre rami ionizzanti. Il computo delle disinte- 
grazioni nucleari con meno di tre rami ionizzanti si può fare solo in via molto 
approssimata, data la scarsità di informazioni sperimentali a tale riguardo; 
nella media atmosfera (1%) le stelle con due rami ionizzanti sembrano essere 
altrettanto frequenti di quelle da tre rami in su, e lo stesso potrebbe assumersi 
per stelle a un ramo e a zero rami ionizzanti. Si sarebbe così portati tentati- 
vamente da stabilire che circa 8 stelle al secondo vengano generate in una colon- 
nina unitaria di atmosfera, ma giova ripetere che questo numero può essere 
sbagliato del 100% in quanto le stelle prodotte negli elementi pesanti sono 


(19) M. MERLIN e O. PIERUCCI: Nuovo Cimento, 8, 528 (1951). 
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differenti da quelle prodotte negli elementi leggeri e quindi la nucleonica che 
produce stelle da tre o più rami negli elementi pesanti può generare stelle a 
meno di tre rami in quelli leggeri, e inoltre il modo di disintegrare degli ele- 
menti leggeri è anche qualitativamente differente da quello degli elementi 
pesanti essendo preferite le fratture in particelle x. Se assumiamo che ogni 
disintegrazione nucleare rappresenti una perdita di energia media di un centi- 
naio di MeV da parte della nucleonica che la genera, otteniamo un contributo 


per il fascio verticale di 
800/x MeV em? s- sterad* . 


La cifra ottenuta deve considerarsi come puramente indicativa, ed è oppor- 
tuno quindi condurre la previsione per altra via e precisamente attraverso i 
dati riguardanti le catture dei neutroni lenti. Il dato sperimentale è qui rap- 
presentato dalle curve di altezza delle catture di neutroni lenti (?%?!) ottenute 
con contatori a BF,. È possibile dal numero assoluto delle catture nel B del 
contatore calcolare il numero di catture nella atmosfera per grammo di aria 
dovute sostanzialmente alla reazione N(n, p)C. Un calcolo in tal senso (?°) 
porta a precisare che il numero di neutroni termici catturati per secondo nella 
solita colonnina di atmosfera di base unitaria ammonta a 3,5. Questo numero 
non è direttamente significativo e bisogna tentare di risalire da esso al numero 
di neutroni prodotti. Uno schema semplificato potrebbe essere il seguente (?°). 

I neutroni di cui ci occupiamo sono emessi nelle disintegrazioni dei nuclei 
di aria da parte della nucleonica, e in piccola parte della componente fotonica. 
I neutroni emessi in tali disintegrazioni avranno uno spettro energetico con- 
tinuo ma con un massimo per energie piuttosto modeste, dell’ordine di pochi 
MeV. Supponiamo di sostituire a questo spettro di generazione uno spettro 
monoenergetico corrispondente ad una energia E, = 4 MeV, e rappresentiamo 
la sorgente di questi neutroni monoenergetici con la funzione 


fy(t)(E — E). 


Si tratta di vedere, tenendo conto del rallentamento, dell’assorbimento e 
della diffusione della atmosfera e dentro la superficie terrestre, quale frazione 
raggiunge la energia termica per poter essere confrontata col dato sperimentale. 
Si giunge cosi a stabilire che per rendere conto del detto numero assoluto di 
catture bisogna ammettere che i neutroni partiti a 4 MeV sieno circa 22 em? s~. 


(°°) LC. L. YUAN: Phys. Rev., 81, 175 (1941); H. J. Kours e L. Yuan. Phys- 
Rev. 86, 129, (1952). 

(21) J. A. Simpson: Phys. Rev., 88, 1175 (1951). 

(22) M. GaLLi e G. Puppi: Relazioni al Congresso di Bergamo della S.I.F. (1952); 
M. Gari: Nuovo Cimento (in corso di pubblicazione). 
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La liberazione di tale quantitativo di neutroni richiede una energia minima 
pari a 264 MeV-cm-?s-. La frazione da addebitare alla primaria compresa 
nell’angolo solido unitario in verticale è pertanto 


84 MeV-em= s“! sterad-!. 


Più difficile è il computo relativo ai protoni, alle particelle x e ai quanti y 
per la mancanza di un quadro sperimentale soddisfacente. Possiamo tentare 
di raddoppiare all'incirca la cifra relativa ai neutroni per tener conto di tutti 
questi contributi. Si giunge così conglobando le perdite di energia dei protoni 
a calcolare in 


248 MeV-em-? s- sterad— , 


il contributo totale delle perdite di energia della componente nucleonica per 
processi non radiativi. Vogliamo ancora notare che il confronto tra i risultati 
ottenuti a partire dalla frequenza delle interazioni nucleari e quelli ottenuti a 
partire dalla frequenza delle catture dei neutroni lenti è soddisfacente; il con- 
fronto infatti porterebbe a concludere che in media vengono generati un po” 
meno di 3 neutroni per disintegrazione nucleare, cifra ragionevole. Il proble- 
ma merita una ulteriore dettagliata analisi sia dal punto di vista sperimentale, 
sia dal punto di vista teorico. È forse questo effettivamente il termine del bi- 
lancio che può subire ancora qualche significativa alterazione. 


7. Bilancio sulle camere di ionizzazione. 


Il bilancio presentato nei paragrafi precedenti era fondato sulle curve di 
altezza delle componenti fatte con i telescopi di contatori, più alcune curve 
di altezza ausiliarie come la crescita delle interazioni nucleari e quella delle 
catture neutroniche; un evidente svantaggio di questo metodo consiste nel 
fatto che gli eventi di bassa energia vanno stimati 0 calcolati in base alla 
conoscenza delle varie componenti fisiche; un vantaggio consiste nel fatto che 
appunto attraverso la conoscenza delle componenti fisiche e dei loro processi 
è possibile stimare le dissipazioni di energia in neutrini, inaccessibili diretta- 
mente. Un bilancio collaterale in un certo senso si può fare con le camere di 
ionizzazione; qui il dato fondamentale è la curva di altezza della ionizzazione 
totale prodotta nella camera omnidirezionalmente da tutte le particelle cariche 
che la attraversano o che vengono prodotte in essa. È facile convincersi che 
in tal modo si tiene conto completamente della componente elettrofotonica, 
delle perdite di energia per ionizzazione delle componenti mesonica e nucleo- 
nica carica, nonchè della energia cinetica delle particelle nucleari provenienti 
dalle disintegrazioni, e dei quanti y di origine nucleare. Si può considerare che 
non sieno comprese nel dato delle camere di ionizzazione le perdite in neutrini 
e l'energia di legame delle particelle di bassa energia evaporate dai nuclei 
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d’aria, nonchè la energia residua al livello del mare. La normalizzazione dei 
dati della camera si fa agevolmente al livello del mare, sotto un opportuno 
assorbitore di piombo, dato che in tali condizioni l’unica componente sensibile 
è la componente u di cui si conosce bene la intensità assoluta e le perdite di 
energia per ionizzazione per unità di percorso. 

Dalla correlazione tra i dati delle camere e quelle dei contatori (2%) è pos- 
sibile stabilire quale sia la ionizzazione totale da far corrispondere alla pri- 
maria verticale. Il nostro bilancio si riferisce approssimativamente ad una 
latitudine di 50° Nord; prendiamo come termine di paragone i valori calcolati 
da NEHER per 51° (Ohama) e notiamo che il valore per la traccia totale di 
157 g-cm-? s7 sterad-!, cui si riferisce il bilancio di Neher, coincide molto 
bene col nostro valore di 155,5 g-cm74s7 sterad-!. Orbene a tale latitudine 
egli calcola per la ionizzazione totale una potenza specifica di 


740 MeV-em= s— sterad= . 


Possiamo senz’altro ricavare dal nostro bilancio la potenza dissipata in 
neutrini, la potenza residua al livello del mare e il contributo dovuto alla 
energia di legame dei nuclei di aria. Sommando questi contributi il bilancio 
delle camere si precisa al modo seguente 


TINMONIZZAZION Re 740 MeV-cm- s71 sterad-! 
RARE, 5 2 5 dano GAS 
In energia di legame . . . . 110 
Residuo al livello del mare . 37 
Motale = RR e2.00 MG Vice sais vera da 


Anche qui ovviamente il dato più incerto è rappresentato dalla energia 
di legame dei nuclei. La cifra da noi introdotta rappresenta il contributo dei 
neutroni computato assumendo che ad ogni neutrone liberato corrispondano 
8 MeV; con questo si ha un contributo di 56; il resto, all’incirca altrettanto, 
è il contributo dei protoni. 

Questa stima non è l’unica, nè probabilmente la più felice, ma rappresenta 
certamente un corretto ordine di grandezza. 


8. — Flusso primario. 


Per poter trarre delle conclusioni sul bilancio delineato nei capitoli che 
precedono, bisogna discutere ancora il problema del flusso primario, da cui 


(2) H. Neuer: Progress in Cosmic Ray Physics, vol. I, Cap. V (North Holland 
Publishing Co., Amsterdam, 1951). 
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dobbiamo dedurre la potenza specifica in ingresso nella atmosfera alla lati- 
tudine in esame. 

Siamo ancora lontani dal poter dare a questo problema una risposta sod- 
disfacente e dovremo quindi accontentarci di una ragionevole stima. Come è 
noto lo spettro primario in ingresso nella atmosfera consta di protoni e di 
particelle x nonchè di una piccola frazione di nuclei di elementi più pesanti 
dell’elio, completamente privi di elettroni. Il calcolo del flusso di potenza pri- 
mario è basato sulla conoscenza dello spettro differenziale della primaria; se 
diciamo infatti f(Z) dE lo spettro differenziale in energia cinetica dei protoni, 
il flusso di potenza trasportato da questa componente si serive immediatamente 


Analoghe espressioni si debbono scrivere per le particelle x e per la com- 
ponente pesante. 

La maggior parte delle determinazioni sperimentali sulla primaria, effet- 
tuate con telescopi di contatori installati in pallone o in razzo, fornisce invece 
un valore integrale del flusso ad una certa latitudine. Poichè però ad ogni 
latitudine corrisponde, in direzione verticale, un ben determinato taglio del 
campo magnetico terrestre, si ha la possibilità, con le dette misure del flusso 
integrale, di costruire uno spettro integrale della primaria. Nella tabella che 
segue abbiamo riassunto alcune delle più recenti serie di determinazioni. 


TABELLA III. 


Taglio | 


Autore Latitudine — (in GeV) Flusso 
| J. A. vAN ALLEN e 8S. F. SINGER (8) 580 0;560 | 0.290 + 0,03 
| 50° 1,8 0,180 + 0,02 
| 410 4,0 0,073 4 0,006 
| 00 14,0 0,028 + 0,004 
| MOVIDA) OSO O a fos 550 1,0 0,296 £ 0,02 | 
410 4,0 0,102 + 0,01 | 
| 280 | 8,1 0,059 + 0,01 
t LS | 
| | | 
| WINCKLER, Stix, DWIGHT e SABIN (€) 510 | 1,6 0,220 + 0,01 
| 419 | 4,0 0,082 + 0,002 
| 310 | 70 | 0,046 +0,0015 
200 IL O 7 0,031 + 0,001 
00 | 14,0 0,026 + 0,001 
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I valori riportati in tabella si riferiscono alla componente primaria globale 
e sono in sufficiente accordo tra di loro, le piccole differenze dovendosi attri- 
buire alle diverse condizioni sperimentali e in particolare alla diversa protezione 
di vari dispositivi per gli effetti spuri. Come avvertito da molti autori, i valori 
del flusso riportati in tabella si debbono intendere come limiti superiori del 
valore effettivo non essendo corretti in alcun modo per l’effetto di albedo, cioè 
per quelle particelle cariche che entrano da sotto in su nei telescopi. Alcune 
misure recenti permettono di stimare il valore percentuale di questo effetto 
di albedo nonchè di valutare il rapporto tra protoni e nuclei più pesanti. A 
latitudine di 55° infatti Ney e THON (?4) trovano i seguenti valori a mezzo 
_di un telescopio contenente un contatore a scintillazione 


Elusso totale™ 2 2 = 2 = « (0:250h cms? sa sterads- 
Flusso protoni . =: .. - (0,190 » 
Elussozogitst Su een, O28 » 


Il valore assoluto del flusso protonico è confortato da un valore isolato 
ottenuto da M. A. CLARK, pure a 55°, di 0,19 cm? s7 sterad — (25). 

Il rapporto protoni/x ‘di 6,8 che si ha dai dati precedenti è il più elevato 
che si sia ottenuto. Una analisi completa della primaria a 41° è stata fatta 
recentemente da PERLOW e coll. (26) con un dispositivo di contatori normali 
e contatori proporzionali. Assieme ad un rapporto protoni/x di 5,3 si hanno 
i valori del flusso riportati sotto, nonchè una valutazione numerica dell’effetto 
di albedo: 


Flusso totale (range > 7 g/em?) . . . 0,079 + 0,005 
Albedo (range tra 7 e 21 g/em?) . . . 0,012 + 0,002 
Flusso Primaria (estrapolata a 0). . . 0,070 + 0,005 
Flusso Proton = e a ee 030082 0/005 
Hluss0r aio os actges Osa eee mates wee Os OL ev OsO.Oi! 
Elusso As 20 See ca 0:001 400001 


Valori assoluti sul flusso delle particelle x e dei primari più pesanti si sono 
ottenuti anche a mezzo di lastre fotografiche inviate in pallone al limite della 
atmosfera; i risultati, discussi da PETERS e coll. (27) sono in generale accordo 


(24) E: P. Ney e D. M. THox: Phys. Rev., 81, 1069 (1951). 

(25) M. A. CLARK: Phys. Rev., 87, 87 (1952). 

(28) C. J. PerLOW, L. R. Davis, C. W. KISSINGER e J. D. SHIPMAN Jr.: Phys. Rev.,. 
88. 321 (1952). 

(2?) M. F. Kapton, B. PerERs, H. L. ReynoLDS e D. M. Rirson: Phys. Rev., 85, 
295 (1952); vedi anche B. PeTERS: Progress in Cosmie Ray Physics, Vol. I, Cap. IV 
(Amsterdam, 1951). i 
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con i dati dei contatori e permettono di definire in modo ragionevole lo spettro 
delle « in ingresso, fino al taglio che corrisponde all’equatore. Per i protoni, 
ui A . 3 ; 
e quindi per la totale, resta molta incertezza dovuta al fatto che non essendoci 


Jon! 5 1,0 5 > = (Gev) 5 10° 
Fig. 6. 


misure dirette è difficile stimare quale sia la correzione da apportare per l’ef- 
fetto di albedo all'equatore. È verosimile, anche dall’esame dell’effetto di 
transizione della totale nella atmosfera che l’effetto di albedo vada percen- 
tualmente crescendo dai poli all'equatore e pertanto ritenere che se esso vale 
il 15% a 41° sia praticamente trascurabile a 55° e sia dell’ordine del 25% 
all'equatore. In tale ipotesi si può tentare di rappresentare lo spettro integrale 
della primaria (fig. 6), senza distinzione tra protoni e nuclei più pesanti, almeno 
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fino a 14 GeV, con un andamento che segue quello delle misure citate in tabella, 
ma sta costantemente al disotto per la correzione di albedo. Per avere infor- 
mazioni sullo spettro della primaria (che corrisponde presso a poco a uno 
spettro di protoni) per energie superiori a 14 GeV, non possiamo utilizzare le 
misure di flusso inclinate al limite della atmosfera, perchè esse sono talmente 
deformate dall’effetto di albedo, variabile con l’angolo di inclinazione, che non 
si saprebbe attualmente come effettuare la correzione. Esistono però alcune 
stime per primari di grande energia, fondate su vari metodi. Tralasciando 
le informazioni che si desumono dai grandi sciami (?8), ricordiamo soltanto 


. il valore stimato da KAPLON e coll. (#9), in base all’analisi di eventi nucleari 


di grande energia, di 4-10-*particelle-cm~? s-! sterad-! per primari di energia. 
intorno a 4,5:108 GeV. 

Si è spesso cercato di rappresentare lo spettro integrale, o gli spettri inte- 
grali delle varie componenti con leggi semplici (*?). Per quanto è lo scopo del 
presente lavoro basta la rappresentazione seguente che consta di due leggi 
di potenza in due intervalli energetici differenti, che non ha ovviamente nessun 
significato fisico, ma si presta ad una rapida valutazione della potenza (*): 

Per 

H<14 GeV N= 0,28/E W= 0,28 (In #— ln E), 
E > 14 GeV N= 1,06/5455 W=3(EN— BN). 


Con lo spettro così come lo abbiamo costruito, il flusso di potenza in ingresso 
che si calcola subito dalle formule della terza colonna risulta in cifra tonda 


1400 MeV-cm-? s- sterad-+ 


per primari protonici di energia superiore a 2 GeV che rappresenta il taglio 
alla latitudine media cui si riferisce il bilancio. Non abbiamo tentato una valu- 
tazione dettagliata della potenza trasportata dalle varie componenti in attesa 
che il quadro sperimentale si chiarisca ancora; tuttavia il valore della energia 
globale in ingresso non può variare di molto. 


9. — Conclusioni. 


Nella tabella che segue è riassunto il bilancio attuale ed è riportato per 
confronto il bilancio precedente: 


(28) P. H. Barrer, L. M. BoLLINGER, G. Cocconi, Y. EISENBERG e K. GREISEN: 
Rev. Mod. Phys., 24, 133 (1952). 

(29) M. F. KapPLON, D. M. Rirson e E. P. Wooprurr: Phys. Rev., 85, 933 (1952). 

(*) Una ragionevole rappresentazione con una unica espressione analitica può essere 
N= 1,2/(E + 2)15 da cui W= 1500 MeV. 
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TABELLA IV. 


Perdite 


PuPPI e DALLAPORTA | Presente 

Dai dati dei contatori | 
| Tonizzazione dei p in aria . . 62 | it 
Sottoterra delw....+..- 37 | 37 
| Decadimento in aria dei 4 . . 190 | 269 
Totali dei y ... + 289 | 383 
In neutrini passaggio T-u. . . 78 | 105 
Gorrezione angolare . ... . | 37 | 62 
Totali in x carichi . | 404 | 550 
ih 978 SORIA Bg Gah oe Boo = 244 | 256 

| | 
i In disintegrazioni nucleari . . 160 | 168 
|  Ionizzazione protoni e x . . . 57 | 80 
Totale generale... | 865 | 1020 


| MeV -cm-2-s-!-sterad—! | MeV-cm_2-s71-sterad 


Dai dati delle camere di ioniz- 

zazione | 
lonizzazione di tutte le compo- 

FONICO eci — | 740 
INENCOLONIA Re ee g — 313 
Inenergia dilegame nuclei di aria | 110 
SOLO LORIA oC ow Os G oc | n SI 

| Totale generale. . . | — | 1200 


| | = MeV -cm-2-s~2-sterad—! 


Dallo spettro primario 


| Flusso di potenza in ingresso . — 1400 
| | MeV-cm~?-s-1-sterad+ 


Come si vede, tutte le voci del bilancio hanno subito qualche aumento, e 
cid, come dicemmo, @ dovuto in buona parte alla diversità delle curve di 
altezza usate. Malgrado si possa essere influenzati dalla circostanza che si sa 
che il flusso in ingresso è ancora più elevato del flusso di potenza dissipato nella 
atmosfera, non riteniamo, nella ricostruzione di questo ultimo, di aver abbon- 
dato per eccesso. La cifra globale è ancora inferiore a quella in ingresso, ma la 
discrepanza resta limitata a qualche cosa come un 20%, che stante le incer- 
tezze da cui risultano affetti i totali, e che non è possibile precisare, potrebbe 
anche non essere significativo. Non si ravvisa quindi la necessità, almeno in 
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base al quadro sperimentale attuale e limitatamente alle latitudini setten- 
trionali, di pensare che la vera fenomenologia sia sostanzialmente diversa da 
quella che è alla base del presente bilancio e che abbiamo tratteggiata all’inizio. 
In particolare, se la generazione dei mesoni k ha luogo sensibilmente (8°) solo 
per primari di energia superiore a 50 GeV, possiamo facilmente stimare, in 
base agli spettri costruiti dianzi, che la energia trasportata dai primari che 
hanno energia superiore a questo limiti è un po’ inferiore a 1/3 della energia 
totale. Ammettendo ancora che all’incirca 1/3 di questa energia vada in k 
carichi, si vede che la sostituzione dei 7 con i k riguarda una energia complessiva 
che è circa 1/10 della energia globale; la relativa correzione pertanto, a favore 
‘del bilancio nella atmosfera, rientra nella incertezza del bilancio stesso e non 
la specificheremo oltre. Ben diversa è la situazione all’equatore dove i primari 
di energia superiore a 50 GeV trasportano più della metà della energia globale 
in ingresso; è possibile che nel bilancio all'equatore il peso dei mesoni pesanti 
sia notevole. Un bilancio sensibilmente equilibrato a latitudini settentrionali 
significa una conoscenza corretta del panorama degli eventi di media energia, 
e riteniamo che in tal senso il risultato raggiunto dal presente bilancio sia 
soddisfacente. Molti punti tuttavia meriterebbero una maggiore attenzione di 
quanta non sia stata loro dedicata nel presente lavoro; e tra l’altro il punto 
che potrebbe ancora subire delle variazioni sensibili è, come già accennato, 
quello che tratta delle perdite di energia per disintegrazioni nucleari, ed in 
genere fenomeni nucleari di bassa energia, dato che il computo di questo con- 
tributo avviene per via indiretta ed è fondato su calibrazioni che possono 
anche non essere esatte. 
Anche una migliore definizione del flusso in entrata dovrà essere tentata, 
nonchè una migliore selezione dei dati sperimentali, nei riguardi della latitu- 
dine geomagnetica, di quanto non sia stato possibile in questo lavoro. 


Desidero ringraziare i proff. B. Rossi e S. F. SINGER per alcune interes- 
santi discussioni sull’argomento e G. BATTISTI per l’aiuto prestato nella pre- 
parazione delle curve di altezza e nel computo delle tracce. 


(3°) R. R. DANIEL, J. H. Davies, J. H. MuLver e D. H. PERKINS: Phil. Mag., 
43, 753 (1952). 
(1) M. A. Pomerantz: Phys. Rev., 77, 830 (1950). 


SUPPLEMENTO AL VOLUME X, SERIE IX DEL NUOVO CIMENTO N. 2, 1953 


Rivestimenti assorbenti per microonde. 


G. LATMIRAL 


Istituto di Fisica Tecnica dell’ Universita, Centro Studi di Radiopropagazione - Napoii 


(ricevuto il 5 Febbraio 1953) 


4. — Analogia tra le equazioni dell’onda elettromagnetica piana e le equazioni 
dei telegrafisti. Conseguente estensione alle onde piane di alcuni concetti 
propri delle linee. Riflessione metallica. 


Sia un’onda e.m. piana propagantesi nel verso positivo dell’asse Z di una 
terna orientata come in fig. 1. 

Sia E, il campo elettrico, polarizzato secondo x, ed H, il campo magne- 
tico, polarizzato secondo y. Le prime due equazioni di Maxwell, che nel caso 
di variazione sinusoidale nel tempo, possono essere 
scritte nella forma (unità Giorgi): 


(1) rotH =oE + joeE = jo (e + 2) E= joe'E, 


(2) rot E = — jouH, 
si trasformano in tal caso nelle equazioni cartesiane: 
oH, 
| STAI: 
i ci ra Fig. 1. 
OH, ‘ Terna destrorsa. 
(4) a jouH, . 


Le (3), (4), sono in tutto analoghe alle equazioni dei telegrafisti: 


dI an ee G 
(5) > oi ee JE (0 oh È) V=— GV—joCV , 
OV ae 3 Ii 4 
(6) da — —jol/1=—jo(L+3;)I-—RI- jotl, 


10 - Supplemento al Nuovo Cimento. 
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scritte con la convenzione che Z sia progressivo, vale a dire misuri le distanze 
a partire dal generatore e non dal carico, come spesso si usa fare nelle linee 
(in quest’ultimo caso i segni delle (5), (6) dovrebbero essere invertiti). 

Se i campi elettrico e magnetico variano nel tempo e nello spazio nel modo 
descritto dalla relazione: 


(7) a= A exp |jw (5) 5 
si ha 

é jo 
(8) è ry. KS 


Il simbolo y designa Ja grandezza 
(9) y= joven, 


che nel caso dielettrico imperfetto diviene 


(10) y = joven = jo (e— 2) w. 
a) 
Le (3), (4), in base alla (7), possono essere scritte nella nuova forma: 
(11) VIEL ot p= E, 


Se il campo elettrico fosse polarizzato secondo y e quello magnetico se- 
condo « si otterrebbero, in luogo delle (11), le relazioni: 


(12) yH, =—joe'E, ; yH, = — jouH,. 
Indicando infine con Z, la grandezza: 
uf» 

13 ARIA PE 

( ) 0 VE | ge jo] b 


è possibile scrivere le (11) e (12) nella forma abbreviata: 


E, Lu 
(14) pg 40 
E, n 
(15) i. È E. IL Zi 
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Nel caso delle linee dotate di costanti uniformemente distribuite L, R, G, C, 
la costante di propagazione y è definita dalla nota relazione (cfr. (5), (6)): 


(16) y=at jp=V(R + joL)(G + jo) = jovL'0'. 


Parimenti Z, equivale all’impedenza caratteristica delle linee: 


R + joL rg 
a GER 


(17) Zp 


L’unica differenza formale tra le equazioni dei telegrafisti e quelle dell’onda 
piana (quando si consideri la tensione corrispondente al campo elettrico, e la 
corrente al campo magnetico) è che nelle seconde manca l’equivalente della 
resistenza serie R (per unità di lunghezza della linea). In realtà nella propa- 
gazione dell’onda piana mancano gli effetti dispersivi dovuti alla resistenza 
serie dei conduttori che, nel caso delle linee, «guidano » l’onda attraverso 
il dielettrico. 

Si noti ancora che l’impedenza caratteristica è positiva quando £ ed H si 
seguono in senso orario (cfr. fig. 1 e formula (14)) e negativa invece quando 
per sovrapporre E ad H occorre effettuare una rotazione di 90° in senso an- 
tiorario sempre guardando nel verso della propagazione (cfr. fig. 1 e for- 
mula (15)). 

L’estensione del concetto di impedenza alle onde e.m. che si propagano 
con e senza la guida di conduttori (1), estensione che nel caso dell’onda piana 
trova immediata giustificazione nelle considerazioni 
sopra esposte, consente di trattare in forma spe- Ei hi Gi 
ditiva un certo numero di problemi e partico- 
larmente quelli che riguardano la riflessione e a 
l’assorbimento delle onde piane in corrispondenza onda piflessa 
di superfici (anch’esse piane) di discontinuità. 

Sia, ad esempio, il caso di un’onda e.m. piana 
che incontri la superficie piana di discontinuità Fig. 2. — Onda piana alla su- 
che separa un mezzo omogeneo di caratteristiche perficie di separazione tra 


3 due mezzi omogenei. Inci- 
€ o,, da un secondo mezzo e£ Oz. Se 5 
AI PIO denza normale. 


E2; Mei 07 


onda trasmessa 


'E2; Ha 


l’incidenza è normale, la condizione che siano con- 
tinue, attraverso la superficie di discontinuità, le 
componenti tangenziali di E e di H, si esprime scrivendo (vedi la fig. 2): 


(18) ee =e AH 


(1) S. A. ScHELKUNOFF: B.S.T.J., 17, 17 (1938). 
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Indichiamo con Z,, impedenza caratteristica del primo mezzo, il rapporto 
tra le componenti tangenziali di H, e di H, nel primo mezzo. Se l’incidenza 
è normale, i valori di dette componenti tangenziali coincidono con i valori 
totali di E, ed H,. Ragionando analogamente per il secondo mezzo si ottiene: 


E E; E 
Z = aa = ss = n= È A 
(19) ol H, Wee. ’ Zoe 


Si noti che se H, ed H, sono assunti di egual segno nell’onda incidente, 
non possono conservare egual segno nell’onda riflessa. Se il vettore di Poynting 
ENH inverte il suo segno bisogna che H, ed Hi siano di segno opposto (19). 

Dalle (18) e (19) si deduce immediatamente, che il coefficiente di riflessione 
del campo elettrico, R,, vale: 


(20) R, === 


Analogamente si ha per il campo magnetico: 


HL, Zu-4 
9 19 1% ol 102 22,28 Pk 
se TE 


Le (20) e (21) sono identiche alle espressioni che, nel caso della propaga- 
zione lungo linee, forniscono i coefficienti di riflessione dell’onda di tensione 
e dell’onda di corrente allorchè una linea omogenea di impedenza caratteri- 
stica Zo viene chiusa su un carico di impedenza Z, ovvero su una linea illi- 
mitata avente impedenza caratteristica Zo. 

Dalle (20) e (21) è agevole dedurre il valore dei coefficienti di trasmis- 
sione A, del campo elettrico e K, del campo magnetico. Si trovano le espres- 
sioni seguenti: 

__ E. 2Zo 


H, 2Zu 
00) Lan, SMS ee A 
Cole ta pie saio: cai Teel 


Il coefficiente di riflessione in potenza vale ovviamente: 
R,=|Rs|:|Rx|=|Rx|}=|Rx}?. 


Analogamente il coefficiente di trasmissione in potenza vale |K,|:|Ky|, che 
però non equivale nè a | K,|? nè a | Ky|?. 


Riflessione metallica. — A titolo di applicazione di quanto sopra detto, esa- 
miniamo il caso di un’onda e.m. piana propagantesi nel vuoto che incida nor- 
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malmente contro una lastra « metallica » (ovvero « ohmica ») di spessore «elet- 
tromagneticamente grande ), vale a dire sufficiente perchè il campo e.m. si 
attenui in misura pressochè totale entro di essa. In una lastra « metallica » 
è o> we e quindi la costante di propagazione y vale (unità Giorgi): 


(23) = jov e'u = iol (« = 5 pa Viouo = Vouo , 
L 


‘— 45° 


Posto Vjouo = « + jf, quadrando ed eguagliando separatamente tra loro 
i coefficienti delle parti reale e immaginaria, si ottiene: 


Perchè la lastra sia tale che nel traversarla un’onda e.m. di pulsazione © 
si attenui, ad esempio 10 volte in campo e quindi 100 volte in potenza, occorre 
che il suo spessore d soddisfi alla relazione exp[— ad] = 1/10, dove « ha il 
valore dato dalla (24). 

Applichiamo ora al calcolo dei coefficienti di riflessione e di trasmissione 
le formule (20), (21), (22), sostituendo a Zu @ Zo i valori delle impedenze 
caratteristiche dei due mezzi. 

Se il primo mezzo è il vuoto e il secondo è caratterizzato dai parametri 
e, u, o, e la pulsazione è tale che risulti o > we, si ha: 


e VA / 
= Wed EA | aa \ Jou 
Vani Ries & | a a | O Zao 


La Z, del metallo può essere scritta nella forma: 
È da Vena ok 
(26) z= | 


L’impedenza caratteristica di un metallo è quindi una grandezza complessa 
che ha eguali i coefficienti della parte reale e della parte immaginaria; ossia 
è un operatore vettoriale che trasforma i campi magnetici in campi elettrici 
sfasati di 45° in anticipo rispetto ai magnetici. 

A questo punto viene naturale la domanda: può la impedenza caratteri 
stica di un metallo essere tale da non provocare la riflessione di un'onda pro- 
veniente dallo spazio vuoto; e, se la riflessione è inevitabile, qual’e il modulo 
della impedenza del metallo che la rende minima? 

Sia go la impedenza caratteristica dello spazio vuoto e Z,= Vje,Kk quella 
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del metallo. Il coefficiente di riflessione del campo elettrico è: 
Viesk — 04 

00 sa Vj OK 

Tenendo presente che: Vjo, E = (00K/V2) + j(00K/V2), si ha: 


(27) R, = 


(28) I abst oa dally 
ee VON sa 
Annullando la derivata rispetto a K si trova che R, risulta minimo per 
K= 1. In corrispondenza di tali valori, quando cioè il modulo dell’impedenza 
caratteristica del metallo è uguale a quella del vuoto pari a 0 = VUE 
= 377 ohm, si ha: 


(29) |R,| =- v2 - > 0,414;  R,=|Re|?= 0,17. 
Dea '/ 2 

La minima riflessione al passaggio vuoto-metallo si ha quindi allorchè 
|Zo|= 00; anche verificandosi questa condizione, la riflessione non si annulla 
perchè la impedenza caratteristica del metallo è sfasata di 45° in anticipo; 
si ha pur sempre una riflessione in potenza pari al 17% circa. Si noti inoltre 
che affinchè risulti |Z|= V(@7/0) = @ ed in pari tempo si verifichi essere 
o > we occorrerebbe usare sostanze ferromagnetiche aventi 4 > yo. 

Se poi o, come avviene nei metalli ad elevata conduttività, assume valori 
tanto grandi da rendere Z, trascurabile rispetto a 00, la riflessione (vedi for- 
mule (20), (21), (22)) diviene totale. 


2. Artifici atti a minimizzare le riflessioni metalliche. Metodo selettivo del 
quarto d’onda. 


Le considerazioni fin qui svolte mostrano come nel passaggio diretto vuoto 
(od aria)-metallo, la riflessione sia inevitabile, qualunque sia la impedenza 
caratteristica del metallo, e risulti praticamente totale quando la condutti- 
vità del metallo è elevata. 

Quali artifici potrebbero ridurre od addirittura eliminare tale riflessione? 

Un ingegnoso artificio per evitare le riflessioni metalliche (?) consiste nel 


(2) S. Ramo e J. R. Wurnnbry: Fields and waves in modern radio (New York 
e Londra, 1944), 

Per evitare le riflessioni « dielettriche » basta interporre tra i due dielettrici da 
«adattare » uno strato spesso 4/4 di un terzo dielettrico avente resistenza caratteristica 
media geometrica tra quelle dei precedenti due (si ricordino le proprietà trasformatrici 
delle linee lunghe un quarto d’onda). 


RIVESTIMENTI ASSORBENTI PER MICROONDE 153 


fare uso della proprietà delle linee in quarto d’onda, di presentare impedenza 
infinita all'ingresso quando sono cortocircuitate all'altro estremo. In sostanza 
il procedimento di cui si tratta è in tutto simile a quello che si dovrebbe 
seguire per far sì che un osservatore, posto all'ingresso di una linea di impe- 
denza caratteristica 0, = 377 ohm (il vuoto), non si avvedesse che la linea 
è stata cortocircuitata all’altro estre- 


mo. ma la credesse infinitamente lunga. a et 2 i SR 
Il non rilevare alcuna onda riflessa 
. DA 2, È . ‘ ae z°Z, 
equivale a supporre che la linea sla g SOT: 4 
infinitamente lunga; nel caso spaziale, dineatHettente linea non riflettente 
a non permettere all’osservatore di*@ve- “ar | nb 3 
Fig. 3. — Artificio per rendere non TI- 


dere» una lastra perfettamente con-  gettente una linea cortocircuitata all’e- 
duttrice disposta normalmente alla di- stremo. 
rezione del suo sguardo. 

Per raggiungere tale fine si può, nel caso della linea, disporre una impe- 
denza Z pari alla impedenza caratteristica Z, della linea alla distanza di 4/4 metri 
dall’estremo cortocircuitato (vedi fig. 3). 

L’impedenza della linea CMND essendo infinita, la linea, vista da AD, 
risulta chiusa sulla sua impedenza  carat- 
teristica. 

La propagazione spaziale avviene con 
modalità analoghe a quelle che si verifiche- 
i rebbero in una linea formata da due nastri 
Pera perfettamente conduttori larghi un metro e 
distanti un metro. La capacità per unità di 


Fig. 4. — Linea «equivalente » : : 3 
allo spazio vuoto, formata da lunghezza di una tale linea risulta infatti 
due nastri conduttori larghi 1 m pari a e, Farad/metro (trascurando l’effetto 
e distanti 1 m. La linea è stata dei bordi) e quindi, supponendo la velocità 


cortocircuitata all’estremo ed a 
4,/4 metri da detto estremo è 
stata interposta tra i conduttori 


di propagazione eguale a quella nel vuoto, 
anche la induttanza per unità di lunghezza 


una lastra puramente ohmica di sarà pari a x, Henry/metro. 
piccolo spessore d e di totale re- Per evitare le riflessioni dovute al cor- 
sistenza R= g/d=3770hm. tocircuito MNOP basterà disporre a Ao/4 


metri dal cortocircuito stesso una lastra pu- 
ramente ohmica avente resistività o e spessore d tali da interporre tra i due 
astri, in totale, la resistenza di 377 ohm. Dovrà quindi essere (vedi fig. 4): 
Z, = 377 ohm = o(1/1-d). 

Lo spessore della lastra ohmica dovrà essere molto piccolo rispetto alla 
lunghezza d’onda, in modo che risulti possibile considerare la lastra stessa 
come una pura resistenza posta a 20/4 metri dall’estremo cortocircuitato (fig. 3). 

Per eliminare le riflessioni per incidenza normale dovute ad uno « specchio 
perfetto » si potrà analogamente, nel caso della propagazione spaziale, disporre 
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a distanza A,/4 da detto specchio una sottile lastra puramente ohmica di spes- 
sore d e resistività o tali che risulti 9/d = 377 ohm= %. Se d è trascura- 
bile rispetto a Z, e se o — 1/0 è > we (e costante dielettrica assoluta della 
lastra sottile), risulterà eliminata la riflessione per incidenza normale delle 
onde e.m. piane aventi quella determinata frequenza. 

Per una più corretta impostazione del problema conviene non trascurare 
lo spessore della Jastra ohmica, ma disporre la faccia posteriore di questa a 
distanza A,/4 dallo specchio e poi considerare la lastra stessa come una breve 
linea resistiva, aperta al suo estremo (ovvero chiusa su una linea cortocircui- 
tata all’estremo e lunga 4/4). 
.  Rammentando l’analoga relazione valida per le linee, si può affermare che 
l’impedenza di entrata Z che l’onda e.m. proveniente dallo spazio vuoto in- 
contra, allorchè perviene alla faccia anteriore della lastra è: 


(30) Z=Z,Ctghyd, 


dove Z, è la impedenza caratteristica della lastra e y la sua costante di pro- 
pagazione. 

Nel caso di una sostanza puramente « ohmica » ovverossia « metallica » nel 
senso che le proprietà dielettriche siano trascurabili rispetto a quelle deri- 
vanti dalla conduttività (o > we), si ha (23, 25): 


IR E pene, 
(31) a V joo a Z= | DIOR f 
lo) 


Se inoltre yd è molto inferiore all’unità si può porre Ctgh yd = 1/yd e 
quindi si ottiene: 


(32) AREAS 


L’impedenza di entrata Z deve d’altra parte essere eguale all’impedenza. 
caratteristica (« impedenza d’onda ») del vuoto. Deve quindi essere: 


(33 ) EU 


espressione in tutto identica a quella già ricavata ragionando sulla linea for- 
mata da due nastri paralleli perfettamente conduttori, larghi e distanti 1 metro. 

Il metodo sopra illustrato per la eliminazione delle riflessioni da parte di 
conduttori perfetti o quasi perfetti può essere applicato anche disponendo 
tra il conduttore perfetto (di cui si vogliono eliminare le riflessioni) e la sottile 
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lastra puramente ohmica un dielettrico di elevata permettività e di bassis- 
sime perdite e ciò allo scopo di ridurre la distanza 4/4. Nel dielettrico ad ele- 
vata e la velocità di fase sarà Ve/s, volte più piccola che nel vuoto, e nello 
stesso rapporto si ridurrà anche la lunghezza d’onda. Se, ad esempio, e fosse 
pari a 4e,, la distanza diverrebbe metà di quella che si avrebbe nel vuoto. 

In definitiva, per eliminare le riflessioni, ad esempio sull’onda di 3 cm, 
da parte di uno specchio perfetto (ad esempio, una lastra di rame), sempre 
supponendo l’incidenza normale, converrebbe, per 
ridurre l’ingombro e il peso del dispositivo assorbente, 
disporre sulla lastra di rame uno strato 4/4 di una 
sostanza ad elevato potere dielettrico. Supposto, ad 


onda piana incidente 


IE 


esempio, e = 4, lo spessore dovrebbe essere di 

(3-10-2/4)-1/2 = 3;75-10* m. ESE 
Sopra detto strato si dovrebbe poi depositare un 407 LIPSIT 

sottile strato di materiale prevalentemente ohmico e fice, ge 

(fig. 5). Posto, ad esempio, d = 0,1 mm si avrebbe AI 24 DA) 


o = 377-107? ohm-m, ovvero 26,5 mho/m. 
Questa conduttivita, realizzabile con adeguate mi- 
scele di polveri conduttrici e dielettriche, risulta tale 
da soddisfare la condizione o > we, semprechè la £ 
del materiale ohmico abbia valori molto elevati. 


de 


è lastra riflettente SS 

Fig. 5. — Metodo « del 

quarto d’onda » per ren- 

dere assorbente una la- 
stra riflettente. 


Critica del metodo del quarto d’onda. 


1) Il sopra menzionato metodo « del quarto d’onda » presuppone la per- 
fetta conduttività del piano metallico di cui si vogliono eliminare le rifles- 
sioni; per detto metallo |Z,|= v(©2/0) deve risultare estremamente piccolo 
e del tutto trascurabile rispetto a 00, così da rappresentare effettivamente un 
cortocircuito. Per una lastra di rame (0 — 1,7:10-* ohm-m) e per l’onda di 
3 cm, posto 2 = wo, si avrebbe: |Z, = V(04/0) = 3,7:10-? ohm, valere asso- 
lutamente trascurabile rispetto a 377 ohm. 

Per lastre di ferro (ad esempio, di ferro-nickel, che ha una resistività circa 
60 volte maggiore del rame) e per frequenze più basse (ad esempio, per 
A, = 10 em), il cortocircuito può ancora considerarsi perfetto, a condizione 
che la permeabilità magnetica non assuma valori molto superiori a fo, vale 
a dire che la suscettività magnetica y sia quasi nulla. 

Questa ipotesi è in generale lecita perchè, salvo che si operi su film di ferro 
elettrolitico 0 comunque su sostanze ferromagnetiche trattate in modo da 
ridurre le dimensioni dei domini elementari di Weiss molto al disotto delle 
ordinarie dimensioni di alcuni microns, la permeabilità relativa declina, nel 
passare da 300 a 30000 MHz, da alcune centinaia di unità all’unità circa (*). 


(3) Vedasi, ad esempio: M. J. 0. StrUTT: Ferromagneties at ultrahigh frequencies. 
in Atti del Congresso Marconiano (Roma. 1948), p. 448. 
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Naturalmente per onde più lunghe e sostanze con 4 > 1 il metodo cade 
in difetto, a meno che si provveda a rendere più conduttiva (ad esempio ra- 
mandola) la superficie in questione, prima di coprirla col rivestimento assor- 
bente già descritto. 


2) Il peso del rivestimento assorbente, per onde centimetriche, può ri- 
sultare dell’ordine di alcuni kg per m?. 
Il metodo può essere ancora valido per onde decimetriche, ma non per 
onde più lunghe, per le quali occorrerebbero spessori e quindi pesi inaccet- 
tabili. 


i 3) Il « film» resistivo assorbente deve essere estremamente sottile (del- 
l’ordine del decimo di mm). 

Questa condizione non risulterà di facile realizzazione, potendo sorgere 
difficoltà nel depositare film così sottili di spessore esattamente tarato e nel- 
l’evitare che siano corrosi o alterati. Appare probabile che sia questo l’incon- 
veniente principale del metodo. 


4) Il metodo è selettivo, vale cioè per un intorno di frequenze limitato, 
in quanto che per dissintonie anche lievi la «linea 7/4» non funziona più da 
« isolatore », ma presenta una reattanza finita. Detta reattanza, a pari dis- 
sintonia, è tanto più bassa quanto più bassa è la impedenza caratteristica della 
linea. L'uso di linee ad alta costante dielettrica consente quindi di operare 
con spessori ridotti ma aumenta la selettività. 


Un vantaggio del metodo « del quarto d’onda » consiste nel fatto che un 
rivestimento costruito con questo criterio, che abbia uno strato dielettrico 
spesso / metri, assorbe le onde piane la cui lunghezza d’onda in quel dielet- 
trico sia non soltanto A= 41 metri, ma anche 4 == 41/3, 41/5, 41/T,..., 
40/(2n + 1) metri. Poco felice appare, sotto questo aspetto, la scelta, per il 
funzionamento dei Radar centimetrici, delle bande di frequenza intorno ai 
3000 ed ai 10000 MHz; un rivestimento di questo tipo che assorba Je onde 
di 10 cm può infatti presentare ancora un considerevole potere assorbente per 
le onde di 3 em. 


3. — Nuovo metodo selettivo basato sull’uso di uno strato assorbente e di uno 
strato adattatore-rifasatore ad esso anteposto (*). 


Il nuovo artificio, di cui si propone l’uso per minimizzare le riflessioni 
metalliche, si basa sulle proprietà trastormatrici delle linee e sull’osservazione 
già fatta (vedasi n. 2, Riflessione metallica, formule (23-26)) che l’impedenza 


(4) « Rivestimento assorbente per microonde », Domanda di brevetto n. 60/77 depo- 
sitata il 21 Giugno 1952 a Roma da G. LarmiRraL e D. DE ANDREIS. 
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caratteristica Z, (o «impedenza d’onda ») di un «metallo », cioè di un mezzo 
per cui o > we, è, così come la costante di propagazione y, una grandezza 
sfasata di 45° anticipo. 

In particolare, è bensì possibile fare |Z|= V(@4/0) = 0, ma la rifles- 
sione ha luogo ugualmente nella misura del 17% circa in potenza, a causa 
dello sfasamento di 45°. Il rivestire lo «specchio » riflettente con uno strato 
« elettromagneticamente spesso » di un metallo avente 1Zo|= 0 non elimina 
quindi le riflessioni. Se lo strato assorbente ha spessore d ridotto, la sua impe- 
denza di ingresso Z = Z, Tgh yd = Z, Tgh (« + ja)d (cfr. formula (24)) non 
si identifica con la impedenza caratteristica Z,; ma dal punto di vista del- 
l'adattamento col vuoto la situazione peggiora ulteriormente perchè lo sfasa- 
mento di Z rispetto all’asse reale aumenta col diminuire di xd, tendendo a 90° 
per ad tendente a zero. 

Per «d molto piccolo si ha infatti Tgh (x + ja)d — (a+ je)d e quindi Z 
risulta sfasato di 45° rispetto a Z,, che lo è gia di 45° rispetto all’asse reale. 

Il problema che occorre risolvere è pertanto quello di trasformare la impe- 
denza puramente resistiva del vuoto in una impedenza reattiva sfasata da 45° 
a 90° in anticipo rispetto all’asse reale. 

Per trasformare una impedenza reattiva in una impedenza puramente resi- 
stiva, o viceversa, conviene, in base alla già menzionata analogia tra onde 
piane e linee, far ricorso alle proprietà trasformatrici delle linee. 

Siano V’ ed J’ le tensioni e correnti d’entrata nella linea di impedenza 
caratteristica Z;; V ed I quelle di uscita ovverossia del carico; sia / la distanza 
del punto considerato dall’estremità della linea a cui è connesso il carico; risol- 
vendo le equazioni dei telegrafisti in regime sinusoidale si ha: 


(34) V'= V Cosh yl + Z'I Sinh yl , 
Vases 

(35) I'= I Cosh yl + 7 Sinh yl. 
0 


Dalle (34) e (35) si deduce immediatamente il valore della impedenza di 
ingresso della linea: 
1 44+Z,Tghyl 


ZIE IA IE 
(36) °Zi+ZTghyl 


Se la linea è priva di perdite, In sua impedenza caratteristica Z, è reale 
(resistiva) e la costante di propagazione y = « + jB diviene pari a y — jb = 
= j(2/2); sicchè indicando con Z == a + jb Vimpedenza di ingresso del ca- 
rico si ha: 


Z- + 4Z, tg 21/7 , (a + jb) + jZo tg 2al/Aa 
ri Vi Reed (Ns ha aie ae = Ds pe Ze — 
(37) 4, Z, + 92 tg 2al/A ° Zi + j(a + jb) tg 2al/A 


= 
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Trattazione del caso più generale. — Prendiamo in esame la formula (37) e 
discutiamone il significato, con riferimento al caso di un’onda piana incidente 
normalmente su una lastra riflettente rivestita da uno strato ohmico assor- 
bente con sovrapposto uno strato dielettrico adattatore-rifasatore. Ricordiamo 
e specifichiamo il significato dei simboli, tenendo presente che le unità usate 
sono sempre quelle del sistema Giorgi. 


Z,= Vp[e è la impedenza caratteristica dello strato dielettrico adattatore- 
rifasatore di costanti e, u= 4%, 0 = 0. 


Z' è la impedenza di ingresso di detto strato, spesso / nietri, quando è chiuso 
sul carico rappresentato dallo strato. assorbente. 


Z=-a+jb=|Z|exp[jg] (dove g= arctg b/a) = Z, Tgh yd è Vimpedenza di 
ingresso dello strato assorbente spesso 4 metri, contraddistinto dalle 
costanti e, u= fo, o (diverse da quello dello strato dielettrico, per le 
quali abbiamo tuttavia usato gli stessi simboli dato che nelle nostre 
formule non ricorrono mai promiscuamente). 


Z, è la impedenza caratteristica dello strato assorbente ed il suo valore è 


Z,= Vul(e— jo/@). Se si verifica essere o > we lo strato assorbente 
ha proprietà «ohmiche » o « metalliche » e si può porre: 


ge pics Ver ne Vee 
/ 7 | ) = ° 
a 20 20 
Analogamente la costante di propagazione y diviene: 
y = Viano = YH! + 5 / 2H = a +i 


Di conseguenza si ha: 


Tgh ad + j tg ad 
1-+jTgh ad-tg xd 


(38) Z= a+ jb = Z, Tgh yd= 


Per ad > 2,3 si può porre Z ~ Z,: Si noti che la condizione xd =: 2,3, che 
assumeremo come condizione limite da soddisfare affinchè lo spessore sia « e- 
lettromagneticamente grande», equivale a quella per cui Vonda e.m. si at- 
tenua di 10 volte in campo nell’attraversare una sola volta lo strato assor- 
bente. 

Per-ad = 7/4 si ha Z= Zy/ 


239.5. 


Ae À rappresentano le lunghezze d’onda nello strato dielettrico adattatore- 
rifasatore e nel vuoto. 
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La fig. 6 indica la disposizione degli strati e delle linee ad essi equivalente 
La formula (37) può essere scritta nella forma: 


,Z+jZotg 2al/i _ 7, a+tib+Ztg 22/4) 
° Z,+jZ tg 2r1/A °(Z, — b tg 2al/A)+ ja tg 2al/d 


(39) =F 


Il nuovo procedimento proposto consiste nell’imporre che Z’ risulti reale 
ed eguale a 0) = V(4,/,)- Si scrivono così 
due equazioni tra i parametri in gioco. 

Per risolvere il problema si può: 


7) 


onda piana incidente 


Dey ae 


©9/)S/Ja eed 
ezuapadwi 


1) noto il rapporto b/a (ovvero noto 
g= arctg b/a), fissare il valore di 1/4 e 


o 
( 


. ita ‘i A : 
quindi ricavare a, b, Z, ’ N E strato adattatore e rifasatore 
oqey.e ! dm 
2) prestabiliti b/a e Z, calcolare al LY 7 
iN Q SUL, strato assorbente”, / 
60s, UA; Vaz 
A 58 WSs riflettente 
3) partire da uno strato o comunque 3 SÒ NN 
. È 3 
da un sistema assorbente dato, vale a dire : : : 
Fig. 6. — Assorbimento delle mi- 


Ae i Sy. Ly oe foes eee Dette croonde con strato assorbente de- 
vare i valori di Z, e dil/A corrisponden-  positato sulla lastra riflettente da 
temente occorrenti per « trasformare » occultare e con strato dielettrico a- 
Zin 0o- dattatore-rifasatore anteposto (nuo- 

Il procedimento 1) fornisce un metodo vo metodo proposto). 
speditivo per risolvere il problema, metten- 
dosi a priori in una data situazione di compromesso tra spessore totale e lar- 
ghezza di banda. 

IN procedimento 2) è affine al procedimento 1) e può in certo qual modo 
essere considerato una variante di esso. 

Il procedimento 3) ha notevole interesse pratico e può essere sviluppato 
analiticamente, oppure, con minore esattezza ma con grande speditezza, grafi- 
camente, mediante l’uso di diagrammi circolari in coordinate cartesiane 0 
curvilinee (carta di Smith). 

Dando per nota la derivazione teorica dei diagrammi indichiamo sempli- 
cemente la soluzione grafica del problema mediante una facile costruzione (5) 
che permette di evitare le forti inesattezze in cui si incorre interpolando su gra- 
fici già preparati (fig. 6). 

In base al grafico di fig. 7, la resistenza caratteristica Zj e la lunghezza I 
di una linea priva di perdite atta «a trasformare » 0, in Z e viceversa pos- 
sono essere trovate operando come segue. 


(5) J. C. SLATER: Microwave transmission (New York e Londra, 1942), p. 29 e segg. 
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Si riportano 0 e Z in un diagramma cartesiano RK, X. Si traccia l’asse del 
segmento 007; sia C il punto d’incontro di detto asse con l’asse R. Si traccia 
Ja circonferenza avente centro in C e raggio CZ = Cop. 

Le intercette di questa circonferenza con l’asse R sono o ed r. 

Si ha Z = vro,. Pertanto, se Z, deve essere reale positivo, » deve anche 
esso essere positivo. È facile provare graficamente che tale condizione non 

è soddisfatta per Z qualsiasi, ma soltanto 

È 3 a, Per |Z|> e@osin|y|/2, dove @= arctgb/a è 

Zr cts F (ovy. I= fr arecte 73) Vangolo di fase di Z. Per g = 45°, ad esempio, 

x deve essere |Z|> 0,382 00. 

; Individuato il valore di Za con Je riserve di 
cui sopra, si congiunge Zj con Z e si traccia 
l’asse del segmento Z;Z. Detta x Vordinata del 
punto di intersezione di detto asse con l’asse 
X si ha 2/Z,= ctg 471/27, donde 1= Af4n- 
. are ctg LIVAR Si noti che, se il punto Z è nel 
Fig. 7. — Diagramma circolare primo quadrante, vale a dire se Z è induttivo, 
per la ricerca della resistenza per |Z|>Z, deve risultare (per costruzione) x 


x| x = Vr 


|° 


caratteristica di Z; e dellalun- positivo; di conseguenza deve essere 47/7 < 7/2 

ghezza 1 di una linea priva E, Mg x y 

di perdite atta a «trasfor- cd SI va ae di x na . ro i 
È = è VE a X negatlv p ez . 

mare» la resistenza Z' = © /8; per |< Cas ce ee / 


nella impedenza Z= a + jb Condizioni opposte si verificherebbero se il pun- 
e viceversa. to Z si trovasse nel quarto quadrante, vale a. 
dire se Z fosse capacitivo. 

È pertanto vero che, dalla (39), posto Z'= Qo e noti due dei quattro para- 
metri VA (reale), 1/A, a, d (ovvero anche a, b/a come da noi supposto), gli altri 
due possono sempre essere ricavati; affinchè le soluzioni trovate siano fisi- 
‘amente possibili occorre però che la scelta dei parametri da prefissare sia. 
fatta coerentemente con le sopra esposte regole e riserve. 


Esempi di calcolo. — Ci proponiamo ora di trattare, a modo di esempio 
alcuni casi specifici, applicando il procedimento 1). 

Nel primo caso supporremo che sia |Z|= Zi; di conseguenza, giusto quanto: 
si è detto alla fine del precedente numero, dovremmo porre 1 = 4/8, vale a. 
dire tg (27//Z)= 1. Supponiamo inoltre che lo spessore dello strato ohmico 
assorbente sia elettromagneticamente grande (vale a dire che sia ad = 2,3) 
e di conseguenza sia Z = a + jb = Z,. Poichè la Z, del metallo è sfasata di 
45° anticipo rispetto all’asse reale, sarà b/a = 1. Noti 1/4 e b/a, caleoleremo in 
base alla (39) i valori della Z, = Vuj]e del dielettrico e della « a » dello strato 
« metallico », dalla quale il valore di |Z,|= V@©y/o è subito deducibile. Posto 
jt == fo ed essendo noto ©, ricaveremo poi i valori della e del dielettrico e 
della o del metallo e quindi gli spessori | e d dei rispettivi strati. 
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Nel secondo caso opereremo nuovamente con xd = 2,3 e quindi con 
bla =1, ma porremo 1 = 0,1764 (con che tg 2z1/A = 2). Risulterà quindi 
LiL. 

Nel terzo caso ammetteremo di operare ancora con |Z|= Z, e quindi con 
1 = 7/8; ma anzichè scegliere xd = 2,3 porremo ad = 2/4; a questa posizione 
corrisponde l’altra (cfr. formula (39)): Z = Z 


0 / 2305 OVVETO 


b sin 689,5 
DOA. 


& cos 689,5 


Riportiamo qui appresso in forma tabulare, i risultati numerici ottenuti, 
confrontandoli con quelli (quarto caso) ottenibili col metodo classico dello: 
strato ohmico di spessore molto piccolo rispetto a A (in pratica porremo 
sempre d = 10 ‘m) depositato su di uno strato dielettrico spesso 7/4 che 
riveste la lastra riflettente. Distingueremo tre sottocasi a secondo che la e 
dello strato dielettrico sia eguale ad eo, 480, 1680- 

La tabella si riferisce a due diverse frequenze: fo = 3000 MHz (A, = 10 em) 
ed fo = 10000 (Ap = 3 em). I dati calcolati sono quelli di cui qui appresso 
riassumiamo i significati. 


IZ, e |Z| sono i moduli delle impedenze caratteristica e di ingresso dello 
strato assorbente (ohmico) di spessore d e di conduttività o > we; g è 
Vangolo di fase di z= a + jb (vy = arctg b/a). 


Zi è la impedenza caratteristica (resistiva) dello strato dielettrico adattatore- 
rifasatore, di spessore J e di costante dielettrica relativa &/€. 


s è lo spessore complessivo dei due strati (s = d + 1). 


| R| è il modulo del coefficiente di riflessione in campo (elettrico 0 magnetico) 
per una frequenza superiore del 10% a quella in esame, a cui teorica- 
mente corrisponde riflessione nulla. Il coefficiente di riflessione in po- 
tenza è ovviamente |R!?. 


Dall'esame della tabella si deducono le osservazioni seguenti. 


a) Nel primo caso la o dello strato assorbente risulta molto bassa e non 
può quindi essere considerata molto maggiore di we (vedi nota (*)); occor- 
rerebbe procedere ad un calcolo di seconda approssimazione, sostituendo a p 
il suo vero valore (che non è di 45° ma un poco inferiore) e quindi ponendo 
ancora 1 = 4/8 ma b/a<1. Il calcolo è stato omesso per brevità atteso che 
i risultati non differiscono in modo sostanziale da quelli tabulati, che a titolo 
orientativo possono essere considerati atti a descrivere questo caso, in cui 
ponendo xd = 2,3 ed 1 = 4/8 si ha modesta selettività ma notevole spessore. 
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b) Lo spessore totale s = d + { può essere ridotto adottando 1 > 7/8; 
così è stato fatto nel secondo caso, in cui si è posto ad = 2,3; 1 = 0,1764. 
Alla diminuzione di s corrisponde però un aumento di |R|. 


c) Una ancora maggiore diminuzione di spessore si ha operando (terzo 
caso) con strato assorbente non elettromagneticamente spesso; vale a dire 
operando non con xd = 2,3 ma, ad esempio, con xd = z/4. In tal caso Z è 
ancora eguale in modulo a Z ma è sfasato di 239,5 in anticipo rispetto a Z. 
e quindi di 230,5 + 45° = 689,5 rispetto all’asse reale. Naturalmente l’essere 
anche d esattamente commisurato a Z, causa un forte aumento della selet- 
tività e quindi di | R|. 

d) Usando il metodo classico (quarto caso) e ponendo e = &; 40; 160, 
si ottengono risultati che, dal punto di vista del compromesso tra spessore 
totale e selettività, sono comparabili con quelli dei tre casi precedenti, otte- 
nuti applicando il nuovo metodo. Il nuovo metodo è però a nostro avviso 
preferibile perchè consente i seguenti vantaggi: 


— Lo strato assorbente è di spessore non decimillimetrico ma millimetrico; 
ad un pari errore assoluto in d corrisponde quindi un minore coefficiente di 
riflessione. Inoltre ded 1 sono paragonabili, laddove con il vecchio metodo d < 1. 

— Lo strato assorbente, che è il più delicato, è posposto e non anteposto 
a quello dielettrico, che meglio si presta a resistere agli agenti atmosferici. 

— Il nuovo metodo consente una varietà di applicazioni molto maggiore 
dell’antico, nel quale praticamente l’unica variante consiste nell’operare con 
strati dielettrici di diversa costante. 


4. — Osservazioni e note conclusive. 


Il metodo di assorbimento proposto, che consiste in ultima analisi in una 
più generale applicazione della analogia tra la propagazione delle onde piane e 
quella lungo le linee, dovrebbe permettere di realizzare rivestimenti assorbenti 
di ridotto spessore in una varietà di tipi molto più vasta di quella consentita 
dall’applicare il metodo classico « del quarto d’onda » inizialmente descritto. 
Questo metodo classico opera con uno strato assorbente, che deve necessa- 
riamente essere di spessore decimillimetrico, anteposto allo strato dielettrico 
adattatore-rifasatore, là dove è generalmente possibile operare con strati assor- 
benti di spessore meno critico, interposti o posposti allo strato dielettrico 
adattatore-rifasatore, pur senza aumentare lo spessore complessivo del rive- 
stimento. 

Tutti questi procedimenti sono selettivi e si può affermare che la riduzione 
dello spessore ha sempre come contropartita un aumento della selettività. 


11 — Supplemento al Nuovo Cimento 
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Raccomandabile, quando si desideri uno spessore limitato, appare essere l’uso 
di uno strato dielettrico spesso Z= 4/8 e di uno strato assorbente « metal- 
lico » con xd = 7/4. 

Alcuni metodi non selettivi sono descritti nella stampa tecnica (978), ma 
richiedono spessi rivestimenti formati da più sottili strati assorbenti a con- 
duttività crescente con progressione geometrica, separati da strati dielettrici 
a bassa permittività. Ovvero implicano l’uso di sostanze con uw + fy (il che è 
fonte di notevoli complicazioni dal punto di vista della preparazione e tara- 
tura delle sostanze stesse), presentanti fenomeni di isteresi magnetica e co- 
munque tali che 4 ed e risultino complesse (senza di che non potrebbe esservi 
assorbimento) e che la radice del loro rapporto sia reale, praticamente indi- 
pendente dalla frequenza ed eguale alla resistenza caratteristica 0, dello spazio 
vuoto in tutta la gamma considerata. Deve quindi essere e = e'— je"; 
u= u'— ju" ed essere inoltre e'/u' = e"/u" = é/uo: L'attuazione di queste 
sostanze risulta ardua, atteso che quelle che possono essere ottenute con pro- 
cedimenti relativamente semplici, come quello di usare gomma impregnata 
con polveri ferromagnetiche, hanno impedenza d’onda relativamente bassa e 
molto inferiore a Qo. 

La e della sostanza risulta infatti aumentata dalla presenza delle parti- 
celle metalliche ferromagnetiche (destinate ad aumentare w ed a creare pu’) 
che si polarizzano sotto l’azione del campo elettrico agente. 

L’uso di sistemi di adattamento toglie d’altra parte a questi rivestimenti 
assorbenti il carattere della non selettività che dovrebbe permettere loro di 
comportarsi come corpi neri ed infirma altresì la loro attitudine ad operare 
con angoli di incidenza qualsiasi. 

I metodi selettivi si applicano soltanto al caso di angolo di incidenza noto 
e costante. 

Nel caso del mascheramento radar, questo angolo può essere assunto pari 
a 90° (incidenza normale). È infatti praticamente soltanto l’incidenza nor- 
male quella che, quando si opera con superfici di dimensioni notevoli rispetto 
a A,, € quindi a reirradiazione molto direttiva e quasi speculare, può causare 
una notevole reirradiazione di potenza nella direzione di provenienza dell’onda 
incidente. 

Se la superficie equivalente del bersaglio Radar è in massima parte dovuta 
ad «aree riflettenti », l’efficacia dei metodi selettivi può pertanto essere note- 


(9) S. A: SCHELKUNOFF : Electromagnetic Waves (New York, 1943). 

(7) L. N. RIDENOUR: Radar Sistem Engineering M.I.T. Radiation Laboratory Se- 
ries, vol. I (New York e Londra, 1947), p. 69 e seg. 

(8) C. G. MoxrGoMmERrY, R. H. DickE e E. M. PurcELL: Principles of microwave 
circuits, M.I.T. Radiation Laboratory Series, vol. VIII (New York e Londra, 1948), 
p. 396 e segg. 
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vole, nella gamma di frequenze per cui sono predisposti. Se invece le « aree » 
diffrangenti » (dovute a corpi piccoli rispetto a Z, od a reirradiazione di lobi 
secondari) assumono notevole importanza, soltanto l’uso di rivestimenti assor- 
benti di spessore sufficiente e di impedenza caratteristica pari a 0, può dare 
risultati soddisfacenti. 

In bersagli caratterizzati da superfici piane od a curvatura grande rispetto 
a A, 1° area riflettente » prevale per solito notevolmente sull’« area diffrangente » 
sicchè l’applicazione di procedimenti selettivi dovrebbe poter essere tentata 
con buone probabilità di successo. 

È bensì vero che una diminuzione dell’area equivalente del bersaglio di 
n volte implica una riduzione della portata Radar di sole Wn volte; ma è 
anche vero che per riottenere la stessa portata occorre aumentare di n volte 
la potenza di punta del Radar. 

Operando su superfici di tipo riflettente, una riduzione, ad esempio, ad 
un decimo del coefficiente di riflessione in campo, che comporta una riduzione 
dell’area equivalente ad un centesimo del valore iniziale, può essere contro- 
bilanciata soltanto moltiplicando per 100 la potenza di punta del Radar il 
che, in generale, non risulterà compatibile con le limitazioni di consumo e di 
ingombro a cui è soggetto l'impianto. 

Segue da quanto si è detto che, tra le contromisure Radar, anche quella 
del mascheramento con rivestimenti assorbenti non può essere sottovalutata 
e che data la relativa facilità con cui sottili rivestimenti assorbenti selettivi 
possono essere realizzati su onde centimetriche, deve essere considerata buona 
norma quella di usare Radar operanti in gamme diverse, o quanto meno atti 
ad essere sintonizzati in gamme di frequenze molto vaste. 

Deve inoltre essere tenuto presente che lo spessore complessivo dei rive- 
stimenti assorbenti è funzione crescente della lunghezza d’onda sicchè il masche- 
ramento per assorbimento risulta in pratica tanto più agevole quanto più breve 
è la lunghezza d’onda. 
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1. — Introduction. 


La technique moderne est caractérisée par la place toujours plus grande 
qu’y prennent les systèmes asservis ou servomécanismes. Apparus d’abord 
sporadiquement comme des solutions ingénieuses de certains problèmes parti- 
culiers, ces systèmes sont devenus suffisamment nombreux et variés pour qu'on 
ait senti la nécessité d’en étudier les caractères généraux; la théorie ainsi formée 
est la Cybernétique, nom autrefois donnée par AMPÈRE a l’art de gouverner 
(les hommes), et repris récemment par N. WIENER en donnant au verbe 
gouverner son sens concret. Le rapprochement des deux acceptions n’est pas 
un simple calembour: en effet, l’étude des servomécanismes a montré que ceux-ci 
étaient capables d’un comportement si complexe que des parallèles troublants 
(et féconds) s’établissaient naturellement avec des activités biologiques et 
sociales des animaux et de Vhomme. Un abus de langage a fait parfois donner 
aux plus complexes des servomécanismes le nom de « cerveaux »; mais il reste 
que le mécanique mord sur le vivant, et que, réservée la part de l’intelligence 
eréatrice, beaucoup d’activités « serviles » de Vhomme apparaissent comme 
justiciables de la théorie des mécanismes asservis. 

Le présent exposé est consacré à l’examen de quelques problemes fonda- 
mentaux de la Cybernétique, envisagés du point de vue du physicien; on n’y 
abordera done pas les implications biologiques et sociales auxquelles WIENER 
attache une importance primordiale; mais il sera commode de poser les prin- 
cipes en ayant présente a l’esprit cette extension a Phumain dont il vient 
détre question. Et, pour commencer, la définition la plus féconde que l’on 
puisse donner d’un servomécanisme sera, croyons-nous, la suivante: un système 
matériel capable, dans un domaine défini, du méme comportement qu’un exé- 
cutant humain, c’est-a-dire qu’un homme a qui il est demandé d’exercer son 
jugement, mais non une activité intellectuelle créatrice. 


(E 
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Analysons le comportement de cet homme: il recoit d’abord des consignes 
plus ou moins précises, relatives & ce que doit étre son activité; puis il est 
confronté A un donné (circonstances extérieures) dont la nature générale est 
connue, mais qui peut varier dans le temps d’une facon imprévisible méme si 
elle est restreinte. Il extrait de ce donné ce qui peut se rapporter aux consignes, 
autrement dit, il sépare l’information recue en utile et en parasite (*). La pre- 
mière se présente, à tout instant, comme un tableau de prémisses, extrait 
d’un ensemble prévisible de tableaux possibles; les consignes sont des règles 
gouvernant les opérations de l’esprit moyennant lesquelles è chaque tableau 
possible de données correspond un tableau de décisions ou actions, extrait 
lui aussi d’un ensemble prévisible de tableaux possibles. La correspondance 
entre le donné et les décisions doit étre univoque si l’on veut que le fonction- 
nement du systéme soit satisfaisant; mais il n’est nullement nécessaire qu'elle 
soit biunivoque: c’est ainsi, par exemple, que l’homme pourra avoir à choisir 
entre les deux branches d’un dilemme, compte tenu de nombreux facteurs. 

Il est essentiel d’observer que l’existence de seuils différentiels de perception 
pour les sens de l’homme impose au donné une structure discrete (cette pro- 
priété reste vraie méme pour un donné purement affectif, mais peut-étre n’est-il 
pas utile de considérer ici cette question). L’ensemble des tableaux de données 
sera donc dénombrable, et en fait il sera fini, A cause de l’existence d’une borne 
supérieure correspondant a la «saturation » des sens. Tl en sera de méme de 
Vensemble des tableaux d’actions commandées par le donné; et a plus forte 
raison cela sera-t-il vrai s’agissant d’un servomécanisme, c’est-à-dire d’un 
système où les données comme les actions sont des grandeurs physiques effecti- 
vement mesurées. 

Du point de vue mathématique, on a done affaire, pour les données comme 
pour les actions, à des tableaux de nombres, c’est-à-dire à des vecteurs, dont 
les coordonnées sont des fonctions du temps; le système qui assure la cor- 
respondance entre le vecteur «action» y et le vecteur « donné » x est un 
opérateur a permettant de représenter l’espace du donné dans Vespace de 
l’action. On a la relation vectorielle fondamentale 


(1) y= ax. 


Pour le physicien, (1) représente un système de transmission multipdle, 
où les coordonnées de x sont les variables d’entrée, celles de y les variables 
de sortie, et x est la fonction de transfert du système. On pourra encore ap- 
peler x et y les «signaux » d’entrée et de sortie; la suite des x dans le temps 
constitue alors le « message » recu, celle des y, le message transmis. D’après 


(1) La notion, ici intuitive, d’information sera précisée plus loin. 
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ce que l’on vient de dire, ces messages sont lun et l’autre a structure discréte. 

On appellera information tout ce qui est susceptible de subir l’opération 
de transfert, ou de constituer le résultat d’une telle opération: les messages 
sont essentiellement des porteurs d’information. Il n’est pas indispensable 
pour l’instant, mais il sera commode, et satisfaisant pour Vesprit, de définir 
la quantité d’information contenue dans un message. Nous postulerons que 
cette quantité est additive par juxtaposition de messages: si H, est l’information 
fournie pas un message A, H, celle d’un message B, le message formé par la 
suite des signaux de B succédant a la suite des signaux de A sera porteur 
dune information H, + H,. Or, si A est choisi parmi m possibles et B parmi 
n possibles, leur juxtaposition peut prendre mn déterminations; la quantité 
d’information satisfera done a l’équation fonctionnelle 


(2) H(mn) = H(m) + H(n), 
dont la solution est 
A(n) = k-log n. 


On définira de la méme facon le débit d’information, c’est-à-dire la quantité 
d’information recue par unité de temps: la quantification s’applique aux 
intervalles de temps tout comme aux amplitudes, et si, pendant l’unité de 
temps, il est possible de mesurer le signal r fois, et que le résultat de la mesure 
soit l’un de n nombres possibles, le débit d’information sera 


(3) kr-log n. 


Enfin, le message pourra étre troublé par des parasites, c’est-à-dire par 
la superposition au signal émis dans l’intention qu’il soit transmis de signaux 
aléatoires, dus au mauvais fonctionnement ou aux imperfections des canaux 
de transmission situés en amont du mécanisme de transfert, ou encore a ce 
mécanisme lui-méme; les deux parasites types sont, en prenant Vexemple 
dun canal de télécommunication électrique, la diaphonie entre circuits (ot il 
y a corrélation entre signaux utiles et parasites), et le bruit de fond (ot les 
parasites sont incohérents). Les parasites sont, comme tout autre message, 
porteurs d’une information parasite, qui devrait, autant que possible, étre 
écartée avant d’influencer le message transmis: c’est le probleme du filtrage 
de l’information. 

Pour important que soit ce probleme, il passe cependant après celui de 
la transmission correcte du message utile, supposé seul regu. Parmi les opé- 
rations élémentaires que comporte le transfert de l’information figurent en 
général des dérivations ou des intégrations par rapport au temps, de sorte 
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que le signal transmis est fonction, non seulement du signal recu, mais encore 
de la facon dont celui-ci varie dans le temps; un mécanisme réglé pour opérer 
une certaine transformation d’un signal statique ou a variation stationnaire 
ne répondra en général pas de la méme facon si la variation est differente de 
celle prévue. Mais si le fonctionnement pour un signal quelconque ne peut étre 
parfait, encore devra-t-il présenter certaines caractéristiques minima de qua- 
lité: le système devra étre stable, c’est-à-dire tendre vers une annulation de 
la variation du signal transmis après que le signal regu aura cessé de varier; 
fidèle, c’est-à-dire fournir, une fois les variations amorties, un signal transmis 
indépendant du mode d’établissement du signal recu; et rapide, c’est-à-dire 
retrouver cet état d’équilibre au bout d’un délai aussi court que possible. 

Ces desiderata sont ceux que doit satisfaire tout système de télécommu- 
nications, et la Cybernétique trouve, dans la théorie de ces systemes, une mine 
de résultats qu'il n°y a qu’A transposer. Mais il faut bien noter que si les 
questions sont qualitativement les mémes, les bornes numériques sont, le plus 
souvent, bien différentes: car les derniers bénéficiaires d’un canal de com- 
munication sont en général l’ouie ou la vue, dont la sensibilité est logarithmique: 
alors qu'un servomécanisme commandera en général un appareil de mesure 
ou un organe plus ou moins linéaire. Il en résulte qu’on exigera du premier 
une précision de l’ordre du néper, du second une précision de l’ordre du 
pour-cent, c’est-à-dire du décinéper. Par suite, seules les techniques les plus 
précises des télécommunications pourront entrer en ligne de compte pour la 
transposition dans le domaine des servomécanismes; en particulier, on opérera 
exclusivement, dans celui-ci, sur les écarts entre signaux d’entrée et de sortie, 
convenablement transformés de maniére que l’écart soit nul dans un fonction- 
nement parfait, et non sur le signal d’entrée brut: c'est le principe bien connu 
de la réaction en radioélectricité. Autrement dit, on transformera le vecteur y 
par un opérateur f qui en assure la représentation dans l’espace des « (en lan- 
gage physique: on transformera les grandeurs de sortie en grandeurs ayant 
respectivement les mémes natures physiques que les grandeurs d’entrée); cet 
opérateur sera de plus tel que fy soit identique a x pour un certain régime, 
statique ou stationnaire, considéré comme le plus probable (« régime normal »): 
alors le vecteur z = x — fy restera constamment petit dans un fonctionnement 
satisfaisant du mécanisme. C'est sur z qu’opérera désormais x, de sorte que 
l’on aura, au lieu de (1), la relation 


(4) har 


Un mécanisme représenté par (4) sera caractérisé, du point de vue de sa 
constitution, par une succession d’élements ou d’étages qui contribuent a 
transformer le signal x en un signal y (chaîne « principale » ou « directe »), et 
par une autre succession parcourue en sens inverse, transformant y en un 
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signal de méme nature que x (chaine «secondaire » ou «de réaction »). La 
formation du signal « écart » x pourra d’ailleurs avoir lieu tout è Pamont de 
la chaîne directe (fig. la), ou en cours de transformation (fig. 1b), et cette 
formation pourra étre immédiate, comme dans les figures 1a et 1b, ou exiger 
un organe spécial discriminateur (fig. 1c). Enfin, chaque organe pourra étre 
lui-méme un servomécanisme (fig. 2). 


Fig. 1 


Les servomécanismes qui ont la constitution que Von vient d’examiner 
sont de la catégorie générale des transmetteurs ou des régulateurs: opérant 
sur des quantités qui doivent rester petites au cours d’un fonctionnement 
satisfaisant, ils sont régis par des équations qui peuvent, dans la plupart des 
cas, étre considérées comme linéaires avec une approximation suffisante. Grace 
à cette simplification, l’étude des problèmes fondamentaux énoncés plus haut 
est, pour eux, assez avancée. C’est cette catégorie seule qui fera done l’objet 
des deux paragraphes suivants; un dernier paragraphe sera consacré a quelques 
considérations sommaires sur les systèmes non linéaires, dont la théorie est 
encore à l’état d’ébauche. 


2. — Transmission de l’information utile dans un système linéaire. 


Si les opérateurs « et f sont linéaires, l’équation vectorielle (4) peut étre 
décomposée en un certain nombre d’équations scalaires, reliant une compo- 
sante de y à une composante de x; ces équations auront toutes la forme 


y = a(e— py); 


où désormais x et y sont des scalaires, et les opérateurs sont des expressions 
linéaires intégro-différentielles, représentant respectivement le gain de la chaine 
principale et le taux de contre-réaction. La correspondance assurée est désormais 
biunivoque, et en vertu d’un théorème connu l’opération (1 + «f) admet un 
inverse; on peut done écrire 


(le op) ei — aa, 


es 
OL 
= 
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et si de plus on suppose les opérateurs « et f(1 + «6)-! commutables, on aura 
pour l’écart 2 l’équation différentielle 


(6) <=e—By=(1+ fa =(1+ 7). 


L’hypothèse supplémentaire que l’on vient de faire est satisfaite, en parti- 
culier, par tous les servomécanismes à structure discontinue, où les opérateurs 
différentiels sont a coefficients constants. 

Le probléme que nous nous proposons d’examiner est, en somme, celui des 
régimes transitoires obéissant 4 l’équation (6): on voit par l’exemple de la 
théorie des circuits électriques que l’outil mathématique le plus adéquat est 
ici la transformation de Laplace ou de Fourier. Rappelons que cette trans- 
formation fait correspondre a une fonction du temps dite « original », soit f(t), 
la fonction «image » 


F(p) = | e-**f(t) dt, 


0 


ou p est une variable complexe, de la nature physique d’une pulsation (la 
partie réelle correspondant è un décrement). L’affixe de p dans le plan com- 
plexe est assujetti 4 se trouver 4 droite d’une certaine parallèle è l’axe ima- 
ginaire, de manière è assurer la convergence de l’intégrale; lorsque p est pure- 
ment imaginaire, soit iw, on a affaire a la transformation de Fourier, et l’image 
F(iw) est le «spectre» de f(t). Inversement, lorsque l’image est connue, on 
obtient Voriginal en calculant l’intégrale 


ia 
f(t) = Smet er! F(p) dp , 
a 


prise le long d’un contour qui enferme les singularités de F, dans le cas de 
l’intégrale de Laplace, ou le demi-plan de droite, origine évitée, dans celui de 
l’intégrale de Fourier. 

Les opérateurs intégro-différentiels linéaires opérant sur les originaux de- 
viennent, lorsqu’on passe aux images, des expressions algébriques en p, et, 
dans le cas des équations (5) et (6), il devient possible d’écrire 


: dp be 
p) = X(p) = I'p)X(p) , 
(7) I) I X(p) = I'(p)X(p) 
1 1 
= X(p) = X 
do 0) agi Bo) I o 


172 TH. VOGEL 


En vertu du principe de superposition applicable aux solutions des équations 
linéaires, toute variation de x(t) peut étre décomposée en variations élémen- 
taires, et il suffira, pour connaître la réponse du système a un transitoire 
quelconque, de connaitre sa réponse a une variation élémentaire. En parti- 
culier, si les coefficients des opérateurs sont constants, comme on l’a supposé 
plus haut, on peut utiliser l’intégrale de Duhamel, 


t 


y(t) = [net du , 


0 


où h(t) est la réponse du système a une impulsion-unité; celle-ci est définie 
(d’une facon peu rigoureuse) comme la limite d’un créneau dont l’amplitude 
croîtrait indéfiniment et la durée tendrait en méme temps vers 0 de manière 
que l’aire tende vers 1; et peut étre représentée, a la satisfaction des mathé- 
maticiens, par une distribution de Dirac-Schwartz d(v) telle que 


pour toute fonction continue f(u). Il en résulte immédiatement, si lon prend 
f(w) = e“, que l'image de 6(é) est 1. 

Supposons done que le signal d’entrée soit une telle impulsion: le servo- 
mécanisme sera stable si le signal de sortie est, de son còté, une fonction du 
temps dont l’amplitude soit négligeable a partir d’un certain instant. L’idéal 
serait évidemment qu'elle soit négligeable dès la disparition de la perturbation 
d’entrée, c’est-à-dire dès l’instant tr + 0; mais si elle l’est, du moins, a partir 
de l’instant t, + t,, on dira que la stabilité esc «assurée a 1/t, pròs ». 

Si lon fait X(p)= 1 dans Véquation (8), l’intégrale d’inversion donne 2(t) 
sous la forme 


a(t) 1 | er! 


_ dp. 

Day) et 1 
Admettons que l’integrandum n’ait d’autres singularités que des pòles p=p,, 
hypothèse qui est toujours vérifiée si, comme il a déjà été supposé, le servo- 


mécanisme est a structure discrète: alors, si R, sont les résidus correspondants, 
on aura 


a(t) = 2ri > Ree, 


de sorte que l’écart z s’amortira à la seule condition qwaucun p, n'ait de partie 
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réelle non négative. Le taux d’amortissement sera d’ailleurs sensiblement donné 
par l’abscisse 1/t, du péle le plus voisin de l’axe imaginaire, dans le demi-plan 
de gauche. Le systéme sera done simplement stable si tous les poles de 
1/(1 + 7) sont situés dans ce demi-plan, et l’étude de la stabilité assurée 
à 1/t, près reviendra è celle de la stabilité simple, pour un servomécanisme 


x 


à opérateur retardé 


1/1 + Tp—1/t)]- 


L’énoneé qui précède constitue, théoriquement, une solution complete de 
la question de la stabilité. En fait, il sera le plus souvent d’application diffi- 
cile, et se prétera mal è la discussion: d’une part, un servomécanisme réel sera 
en général composé de nombreux éléments ou étages, dont certains possè- 
deront leur chaîne de réaction individuelle, et dont certains, rebelles a une 
étude théorique, ne seront connus que par les caractéristiques, relevées en 
laboratoire, du module et de l’argument de la fonction de transfert en régime 
harmonique p = im; de l’autre, on désirera en général connaître influence 
sur la stabilité des valeurs attribuées à certains paramétres de dimensionnement, 
pour une structure qualitativement la méme. La détermination numérique des 
poles, et, plus encore, leur expression littérale en fonction des paramétres, 
seront en général des opérations inextricables. Aussi a-t-on imaginé divers 
critéres plus commodes, déduits du résultat fondamental précédent, et repo- 
sant tous sur le théorème connu de Cauchy: Si le lieu de Vaffixe de p dans 
le plan complexe enferme n pdles et m zéros de la fonction analytique f(p), 
le lieu de l’affixe de cette fonction fait au total (m — n) tours autour de l’ori- 
gine. (Il faut entendre que m —» est la somme algébrique des révolutions 
accomplies dans le sens + et dans le sens —). C'est HURWITZ qui semble 
avoir formulé le premier critére fondé sur cette proposition, dans un travail 
dont on n'a généralement retenu que le còté algébrique (°); la forme la plus 
communément utilisée aujourd’hui est celle de Nyquist: on prend pour 
contour fermé parcouru par l’affixe de p celui qui délimite le demi-plan de 
droite, c’est-A-dire Vaxe imaginaire, moins l'origine contournée par un cercle 
de rayon infiniment petit, et un demi-cercle de rayon infini; ce contour ne 
doit pas contenir de pòles de la fonction de transfert, l’affixe de cette fonction 
devra donc faire, autour de l’origine, autant de révolutions dans un sens que 
dans l’autre. En écrivant la fonction sous la forme 1 + 7(p), on trouve que 
l’affixe de 7(p) devra faire au total un nombre de révolutions nul autour du 
« point critique » — 1 + 70. 


(2) A. Hurwitz: Math. Ann., 46, 273 (1895). Nous devons cette remarque à M. 
Luc GAUTHIER. 
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Sans entrer dans les détails d’application, indiquons comment on peut 
tracer le lieu de 7(p): au petit demi-cercle autour de Vorigine correspondent 
les points è l’infini, d’argument — mz/2 si m est Vordre du pole a Vorigine 
de T(p), reliés par m demi-cercles de rayon oo; le grand demi-cercle de rayon 
infini donne l’origine, avec une tangente dans la direction — m'7/2 si m' est 
Vordre du zéro a Vinfini de 7(p); enfin, l’axe imaginaire a pour image la 
caractéristique de transfert en régime harmonique, 7'(iw), complétée par sa 
symétrique pour les valeurs négatives de la pulsation. Cette caractéristique 
se calcule aisément pour les éléments simples, ou bien peut étre relevée expé- 
rimentalement, modules et phases séparés, pour des organes plus complexes (3). 
La caractéristique totale s’obtient alors en observant que la fonction de trans- 
fert est le produit de celles des étages composants connectés en série; et lorsqu’il 
y a des chaines de réaction, qui forment une connexion en paralléle, on aura 
recours 4 l’expression de 1/I(p), 


Fy = BO) + aq = BO) + Fe) 


Ip) A(p) T(p) 


qui permet de séparer les opérateurs des deux chaînes. C'est ainsi, par exemple, 
que l’on aurait, pour le circuit de la fig. 2, 


al! alt al 1 
sani DA TI; Bs (Der oa (Ba aa, ’ 
Tw + apy EX) + aca +75) 
d’où 
1 il 1 f sl 1 
SE L.= 3 SMR aR 2 
Tip) Bip) DM) Ap) | 0) + am |) aa] 


x 


Le critere de Nyquist est particulicrement adapté a la discussion de ce 
qu’on pourrait appeler la permanence de la stabilité: si le fonctionnement du 
servomécanisme se trouve à quelque moment étre défectueux par vieillissement 
des organes, ou par tarissement de la source d’énergie, les modules des fonctions 
de transfert ne peuvent que décroître; le comportement du système pourra 
alors revétir deux aspects bien différents, suivant la position du point critique 
par rapport au lieu normal de 7(p): si ce point était extériewr au cercle de 
rayon vecteur maximum tangent A cette courbe, il le restera, et le système 
affaibli sera encore stable; mais si l’annulation du nombre de tours total autour 
de ce point n’était obtenue qu’en équilibrant des révolutions de sens contraires, 
une déformation du lieu de 7(p) pourrait détruire cet équilibre, et la stabilité 


(*) Toutefois, lorsque l’élément n’est justiciable que d’une étude expérimentale, il 
faudra aussi la faire porter sur sa lingarité. 
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d’un tel système serait quelque peu précaire. On comprend alors que, se 
placant 4 un point de vue pratique, certains auteurs aient pris pour indication 
du degré de stabilité les « marges » de module (ou gain) et d’argument (ou 
phase), exprimées par l’abscisse du point où la courbe de Nyquist coupe le 
demi-axe réel négatif et par sa pente en ce point (fig. 3). 


T(W)=x+iy 


Si |T|\>\A\|en CO E,F,G,H, 
F(MI)=-1 


Fig. 3. Fig. 4 (LEFÈVRE). 


On aura observé, d’autre part, que la stabilité étant déterminée par le 
nombre de révolutions que fait le rayon vecteur, l’argument de celui-ci parait 
jouer un ròle privilégié, plutòt que son module; il semble d’autant plus sou- 
haitable de conduire la discussion de ce point de vue que l’on sait réaliser, 
dans plusieurs domaines de la physique, des organes correcteurs de phase, qui 
ne modifient pas le gain d’une transmission de facon sensible. Or, la chose 
est a priori possible, puisque d’après le theorème de Bayard (ultérieurement 
repris par BoDE) il existe des relations entre les modules et les arguments d’une 
fonction de transfert, pour un système à structure discrète. 

M. DEMONTVIGNIER et P. LEFEVRE ont été ainsi conduits à formuler un 
critere extrémement intéressant, qui ne suppose que la connaissance de la 
caractéristique arg T(im), et celle des modules pour les seules pulsations ou 
Vargument prend une valeur dite « remarquable ». Si l’on met l’équation aux 
valeurs propres sous la forme 


p -—- ap’ —- =. 


p? + bps +... “ela 


ce qui est toujours possible avec telle approximation que l'on veut, la suite 
des valeurs remarquables est {2kz} si A>0, et {(2% + 1)x} si A<0, avec 
k entier. 

Dans ces conditions, la forme définitive du critère est la suivante: on 
compte pour + 2 tout point remarquable où Vargument décroît, pour — 2 
celui où il décroit, si on y a |T|>|A|; dans les mémes conditions, les points 
remarquables extrémes (@ = 0 et @ = co) ne comptent que pour +1, et 


» 
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ceux où |7|<|A|, pour 0. Le total des points ainsi comptés est égal &4 m— n, 
où m est le nombre de valeurs propres dans le demi-plan interdit (done: celui 
des solutions instables), et n est celui des pòles de 7 qui s’y trouvent. A chaque 
valeur de |A| correspond une valeur de m—vn calculée d’après cette régle, 
soit f(|4]); mais le nombre n est fixe: le système sera donc stable si m = 0, 
c’est-à-dire si 


n=—f(JA]). 


On trouvera dans l’étude de LEFEVRE citée in fine des détails de discussion 
sur quoi nous ne pouvons nous arréter ici; contentons-nous de faire remarquer 
‘que ce sont les parties de caractéristique a pente positive qui jouent un ròle 
stabilisateur, en diminuant m—n; si un système se révélait instable, on 
pourrait le stabiliser (sans rien changer aux modules |7(iw)|) en ajoutant 
des étages correcteurs de phase, qui introduiraient des pentes positives con- 
venables (4). 

La fig. 4 illustre l’application de la règle ci-dessus. 

Passons a la deuxieme question soulevée par le fonctionnement du servo- 
mécanisme en régime transitoire, celle de la fidélité. Il sera commode, pour 
la discuter, de considérer une modification élémentaire brusque, mais per- 
manente, du signal d’entrée, par exemple de prendre pour w(t) la fonetion- 
échelon 


2 1 + sent 
rat, 
dont la transformée de Laplace est 1/p. 
Le système sera fidéle si l’écart tend vers 0 pout t + co. Or, un théorème 
classique du calcul opérationnel donne la relation 


lim w(t) = lim pU(p), 


t>o DAI 


done ici 


SPE eo 
+ FO) 


A 


système sera donc fidéle si ède un pdle è l’origine; et si dle a 
Le syst sera done fidéle T possede l l’origine; et si ce pdle 


(*) Ainsi, un système stable peut fort bien comporter des parties qui, prises iso- 
lément, ne le seraient pas. Remarquons a ce sujet que lorsqu’on a affaire à un système 
pluridimensionnel, un couplage entre les chaines peut permettre de stabiliser l'ensemble, 
alors que des chaines séparées seraient instables: ce probléme a fait l’objet d’une étude 
récente de GoLomB et USDIN. 
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la multiplicité s, la transmission sera encore fidéle pour les dérivées succes- 
sives du signal d’entrée, jusqu’a celle d’ordre (s— 1). 

Il faut enfin que la réponse du mécanisme aux variations soit rapide: mais 
ici, la discussion est plus délicate, car le signal de sortie présente toujours, 
par rapport au signal d’entrée, une distorsion de forme qui dépend du signal 
d’entrée, ce qui rend quelque peu arbitraire la définition du retard. On a 
souvent fait le raisonnement suivant, qui est un peu simpliste, mais présente 
Vavantage de conduire 4 un résultat donnant une idée vive de la facon dont 
les choses se présentent qualitativement: 

Tous les mécanismes présentent le caractère commun de mieux transmettre 
les signaux harmoniques dont la pulsation est comprise entre certaines bornes 
ce sont donc en somme des filtres plus ou moins efficaces; de sorte qu’on 
pourra utilement rechercher ce qui se passe pour un filtre idéal. Supposons 
d’abord qu’on ait (cas du filtre passe-bas) 


T(io = exp [— iwty] |\0|< @ 
| Pim) = 0 |0|> a 
la réponse & la fonction-échelon est alors donnée par Vintégrale de Fourier 


(0) @o(t— to) 


al 1 [si t— Magi sinw 
y(t) == 4 LON to) dw =-—+- == dhe 
Dh” ai 0) Py ah 


U 
0 


On pourra prendre pour critére de la rapidité la pente maximum de la courbe 


y(t), puisque celle de x(t) était infinie au temps 0: soit 


dy @) sinu Wy 
v= | = Max meee "Oo = Wee == 
dt max dao(t a to) TT U IT 


Il est visible qu’une translation de l’origine des © ne modifie en rien le raison- 
nement: les résultats précédents valent done également pour un filtre passe- 
bande de largeur de bande ,. Ainsi, plus un transmetteur est sélectif, et plus 
il est long & rendre la valeur correcte du signal d’entrée: on verra plus loin 
l’intérét de la sélectivité du point de vue du filtrage de l’information parasite, 
il ya donc là une grave antinomie, qui affecte tous les servomécanismes linéaires. 
Par contre, certains types de systèmes non linéaires en sont affranchis, comme 
nous avons montré ailleurs (systèmes déferlants). 

Il est bien évident que les considérations qui précédent sont trop schéma- 
tiques pour servir & autre chose qu’à une orientation générale. Le problème 
du choix correct des paramétres, dans un mécanisme de structure donnée è 
Vavance, peut étre abordé plus sérieusement en se proposant de rendre mi- 


178 TH. VOGEL 


nimum l’écart quadratique moyen entre les signaux d’entrée et de sortie: cette 
facon de procéder, inspirée des travaux de WIENER dont il sera question au 
paragraphe suivant, a été adoptée par AIGRAIN et WILLIAMS dans une étude 
des amplificateurs, qui s’applique telle quelle aux servomécanismes linéaires 
transmetteurs en général. 

Admettons que le mécanisme idéal transmette le signal d’entrée en repro- 
duisant exactement sa forme en fonction du temps, avec seulement un certain 
retard #,; l’écart quadratique intégré entre le signal de sortie y(t) du méca- 
nisme réel considéré et celui du mécanisme idéal sera alors 


(ce) 


= | [y(t) — a(t ty) at ; 


0 


il suffit d’ailleurs, comme précédemment, de considérer le cas où w(t) est la 
fonction-échelon. Si l’on pose 


y(t) = y(co) + y(t) , 


la fonction y,, nulle è Vinfini, pourra commodément, suivant une méthode 
classique de la physique mathématique, étre développée en série de fonctions 
orthogonales formant une suite fermée ou « base » dans l’intervalle 0 < t < 00; 
si {g,} est une telle suite normée, on pourra poser 


x(t =5, ty) =1 Sa > EnPn(t) ’ 
Yi(t) = > UnPn(t) o) 


et puisque par hypothèse 


foo} 


[rupe di = Onn ’ 


0 


on aura 


I= | [Eto —Eadpalt)lf at = mt Ea). 


0 


Les coefficients €, sont des fonctions du retard #,, encore indéterminé; les 7a, 
de leur còté, dépendent de la fonction de transfert, où nous supposons que 
figurent j paramétres de dimensionnement ¢, ..., ¢;: les valeurs minimisantes 
de ces paramétres et de £, seront les solutions du système de j + 1 équations 
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à j +1 inconnues 


Explicitons ce qui précède: on connaît, dans l’intervalle (0, co), la base des 
polynomes de Laguerre 


de fonction génératrice 


exp 


pi 
= = Tha i n» 
ils sont orthogonaux relativement au poids e*, de sorte que nous pourrons 
prendre pour base les 


orthogonales dans (0, co) relativement au poids 1. 

La transformée de Laplace de ,(t) est 2/(2p + 1)[(2p —1)/(2p + 1)]", ou 
encore (1— s)s", en posant 2p = (1 + s)/(1—s); de sorte que si Y,(p) est 
Vimage de y,(t), on aura 


Mais 


de sorte que 
2 11+ 9p 
6 resa 


et 


(10) 


12 - Supplemento al Nuovo Cimento. 
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Les n, et les 07,/0c; sont donc connus par les développements (9) et (10) 
lorsqu’on connaît l’expression littérale de /"(p). Quant aux É,(t), ils peuvent 
de leur còté étre calculés une fois pour toutes, pour la fonction-échelon; on en 
trouvera une table dans le mémoire d’AIGRAIN et WILLIAMS. 


38. — Filtrage de l’information parasite dans un système linéaire. 


Le deuxiéme probleme fondamental de la théorie des servomécanis mes est 
celui de leur comportement en présence d’informations parasites. Il est bien 
évident que lorsqu’on se trouve en présence de deux messages concomitants 
dont lun est a rejeter, la discrimination n’est possible que s’il existe entre eux 
quelque différence caractéristique qui soit stationnaire au cours du temps, au 
moins statistiquement: c’est ainsi, par exemple, qu’on pourra supposer aux 
parasites une amplitude moyenne très inférieure a celle des signaux utiles, 
ou un spectre moyen plus étalé, etc. Cette hypothèse de stationnarité ne sera 
pas toujours vérifiée, et par suite les considérations qui suivent ne sont pas 
nécessairement applicables a tel probleme particulier de servomécanisme; mais 
les cas où elle est valable sont encore suffisamment nombreux et importants 
pour que les conclusions auxquelles elle conduit méritent attention. 

Nous commencerons par nous placer au point de vue informationnel, pour 
examiner un cas schématique dans Vhypothése où le rapport du signal para- 
site au signal utile est borné supérieurement en amplitude. Comme il a été 
dit dans l’Introduction, le nombre des valeurs possibles du signal d’entrée a 
chaque instant est limité par les possibilités de discrimination de Vorgane 
mécanique ou humain qui les apprécie: si une grandeur d’entrée est le résultat 
dune mesure faite à 1% près, elle s’exprimera par l’un des nombres entiers 
compris entre 1 et 100. Trés souvent, la grandeur de sortie sera choisie parmi 
un nombre de possibles encore plus restreint. 

Supposons alors que le servomécanisme se comporte comme un filtre passe- 
bas |w@|<a,: le spectre du signal y(t) sera alors limité è 


Da tof: 
Mi (e DI I y(t) exp [— tot] dt, la|< Wi, 


— Wy, 


et si on le développe en série de Fourier dans cet intervalle, ses coefficients 
de Fourier seront 


il il nina 
Crp = —— I Y(@) exp [—inzw/o,|do = an dl = | : 
sl 


(OX 


di 
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le signal y(t) sera done entièrement déterminé si l’on connaît sa valeur @,/z 
fois par unité de temps. On ne fournirait pas plus d’information en faisant 
des observations plus rapprochées (qui seraient « redondantes »). 

Si l’on veut pouvoir choisir entre n décisions possibles (transmettre une 
grandeur avec une précision de 1/n, par exemple), il faudra disposer, en vertu 
de (3), d’une quantité d’information 


Wi 
HErRolozini 
IT 


par unité de temps. Or, le signal d’entrée peut étre considéré comme transmis 
au mécanisme par un canal amont qui, de son còté, ne laisse passer que les 
pulsations inférieures à w,, et où le message utile est mélangé a des signaux 
parasites dont l’amplitude peut atteindre au plus 1/q-eme de celle des signaux 
utiles: l’information fournie au mécanisme sera donc 


We 
leh = ee ‘logq , 


par unité de temps, et comme le mécanisme ne peut évidemment pas créer 
d’information, on aura 


HH, soln! io 


On peut tirer de cette importante formule diverses conséquences, dont cer- 
taines n’étaient pas évidentes a priori; sans insister ici sur cette discussion, 
notons que si le rapport signal/bruit de fond est donné, on peut arriver a 
transmettre l’information désirée, soit en élargissant la bande passante en 
amont des sources de bruit, ou en la rétrécissant 4 l’aval; ou encore: l’augmen- 
tation du rapport signal/bruit de fond permet de rétrécir la bande du signal 
recu, à information constante. 

C'est encore A la théorie de l’information que l’on recourra pour examiner 
la possibilité de compléter éventuellement un message insuffisamment infor- 
matif è cause des parasites qui l’altèrent, au moyen d’une chaine de cor- 
rection; mais il faut auparavant étendre la notion de quantité d’information 
au cas de messages aléatoires non indépendants entre eux (ceux envisagés 
dans l’établissement de la formule (2) étaient essentiellement indépendants). 

Si un message est composé de n signaux, dont chacun est choisi parmi m 
possibles inégalement probables, de probabilités respectives pi, ...: Pm; la pro- 
babilité totale d’un message de durée 7° sera 


= mA, PT NPT 
P= PP Di Bin™ 


i 2 
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et le débit d’information par unité de temps sera 
(11) —n Y p; log p;. 


De méme, si l’on a affaire a deux espèces de signaux possibles, x et y, et 
que p(i, 7) soit la probabilité composée d’un état x,, y;, le message x, y sera 
porteur d’une quantité d’information 


(12) — ¥ pli, j) log pli, j) 


Enfin, si les deux espèces de signaux ne sont pas indépendantes, et que, 
x prenant la détermination 2,, y ait la probabilité p;(j) de prendre la déter- 
mination y;, on aura 


(13) p.(j) = ptt, i) / > vp, J). 


7 


En moyenne, l’information sur y pour chaque état de x sera, compte tenu 
de la probabilité de cet état, 


H:(y) = — > pli, j) log pj) = H(2, y) — H(a), 


i,f 


en vertu de (13). 
De méme, 


H,(v) = H(«, y)— Hy) . 


L’information conditionnelle H,(v) est ce que SHANNON appelle Vambiquité 
du message recu. 

Supposons alors un message altéré, suivi d’un message correcteur codé, 
dont chaque signal indique si le signal de méme rang du message principal 
doit étre corrigé: ce deuxiéme message sera porteur d’une information con- 
ditionnée par l’état de réception du précédent, et lon peut conclure que l’in- 
formation totale pourra étre ainsi restituée, a la seule condition que le canal 
de transmission soit capable d’un débit d’information égal a H,(2). 

Par exemple, si le message consiste en une suite de signaux binaires dont 
chacun a une probabilité p pour étre altéré (done une probabilité q = 1— p 
pour étre exact), les erreurs seront sans effet si ’on fait connaitre les signaux 
qu’elles affectent: or ceci suppose l’émission de binaires a probabilités p pour 1 
et q pour 0, done un canal débitant — (p-log p + q-logq), ce qui est bien 
la valeur de l’« ambiguité » du message primitif regu. 
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Nous allons maintenant reprendre la question du filtrage des parasites, 
comme nous l’avions fait pour celle de la transmission correct> du message 
utile, en partant de la considération de erreur quadratique moyenne; l’exposé 
très sommaire qui suit est essentiellement celui des résultats de N. WIENER 
et de ses collaborateurs L’hypothése de base étant celle du caractère statisti- 
quement stationnaire des messages utile et parasite, il convient de rappeler 
d’abord quelques notions usuelles de la théorie des fonctions aléatoires. 

Si p(x) est la densité de probabilité d’une fonction aléatoire x(t), telle que 


| pP@) de = 1, on appelle fonction de répartition de x Vintégrale 


F(x) = | pla) de, 


—_ 
et moment centré d’ordre n 


(co) 


M,() =|, xp(x) dx ; 


—-@ 


pour n= 1, on a l’espérance mathématique (ou valeur moyenne) de x, que 
Yon désignera par E[x(t)]; Vesperance mathématique du produit de deux 
fonctions aléatoires est égale au produit de leurs espérances mathématiques 
individuelles : 


E(c-y) = E(w): BYy) - 
On appelle covariance de x(t) Vexpression 


C(t, t') = Elx(t)-a*(@)], 


où l’astérisque désigne les quantités imaginaires conjuguées. Enfin, la caracté- 
ristique de x(t) est la transformée de Fourier de p(x), 
4 È ; (co) 4” "M x 
P.(0)= B(exp fioz])= | esplior]p() de = Vie 
5 ! 


-@ 
On a évidemment 


foo} 


it 
p(a) = =| E(exp [iwa]) exp [— ior] doa , 


_ o 


et la loi de composition des caractéristiques pour les éyénements aléatoires 
concomitants 


E(exp [iw(w + y)]) = E(exp [iw] exp [twy]) = E(exp [iwx])- E(exp [twy]) - 
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On pourra encore écrire sous forme polaire, 
Pruyy(@1, 03) = E(exp [im,2])- E(exp [iw,y]) . 


La fonction aléatoire x(t) est dite stationnaire d’ordre 2 (qualification que 
nous sous-entendrons désormais, car il ne sera pas question d’autres modes 
de stationnarité) si E[x(t)] = 0, et si la covariance C(t, t+h) est une fonction 
de h seul, soit 7(h). Cette fonction est dite fonction d’auto-corrélation de x(t), 
car il est aisé de voir que son expression se confond avec celle, usuelle, du 
coefficient d’auto-corrélation lorsque la fonction est parfaitement discon- 
tinue (5). Nous supposons (hypothèse ergodique) que 


Si Pon fait la transformation de Fourier 


A) = | Ro exp [tmh|da , 


—_ 0 


1 ey 
R(w) = dl r(h) exp [— twmh] dh , 


_ 0 


R(w) est la densité d’énergie spectrale de x(t) (Théorème de Bochner). 
Nous définirons de méme, pour deux fonctions aléatoires stationnaires x(t) 
et y(t), les fonctions de corrélation mutuelle 


r(2,Y; h) = [amore +h)dt, 


(col 


(5) Toutes les définitions qui précèdent s’appliquent aussi bien aux variables con- 

tinues qu’aux discontinues, en entendant l’intégrale au sens de StIeLTJES. Notons 
(ce) 

aussi, pour n’y plus revenir, que le symbole | Adt désigne ici de fagon abrégée l’opé- 


-_ 0 


ration 
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et de densité d’énergie spectrale mutuelle 


(ce) 


if 
Re, y;@) =s= r(z, y; h) exp[—imh]dh. 


On remarquera que si l’on sépare le module et l’argument d’une fonction 
aléatoire complexe x(t), sous la forme x = Xe‘, et si ne change pas au cours 
du temps, cet argument disparaît de l’expression de r(x; h), où n’intervient 
dans l’integrandum que vz* = X(t)X(t + h). Corrélativement, le transfert des 
énergies spectrales 4 travers un mécanisme d’opérateur y obéit a une loi 
remarquablement simple: si lon écrit le signal transmis y(t) sous la forme 
de Vintégrale de Duhamel, 


foe} 


y(t) = | Mat) du, 


—_ 2 
on voit que 


co [co] o (ce) 


r(y;t)= | y(f)y(t+7)dt = | du | de HR | x(t—u)x(t+T—v)dt| = 
= | | h(u)h(v)r(a2; u—v +7) dude. 


Prenons les transformées de Fourier des deux membres: observant que h est 
la réponse à la distribution de Dirac-Schwartz, d’image 1, il vient 


Ry; @)=|IF(iw)|?-E(c; @) . 
L’image de l’opérateur n’intervient plus que par son module. 
Un calcul tout semblable donne, pour le cas où a(t) est la superposition 
de deux signaux (l’un utile, l’autre parasite) 
a(t) = u(t) + bit), 
Rly; o) =|T(iw)|:-[R(u; @) + RO; ©) + 2R(u, bi; ©] - 
Le terme R(u, b; @) est nul dans les nombreux cas où le signal utile et le 


parasite sont incohérents. 
Nous pouvons maintenant exprimer Vécart quadratique intégré du signal 
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transmis et du signal utile retardé, en faisant intervenir les fonctions de cor- 
rélation que l’on vient de définir: - 
En écrivant y sous forme d’intégrale de Duhamel de x = wu + b, il vient 
. pour l’écart quadratique intégré 


foo} 


r= [ly utt at = 
[| h(c)h(x'[u(t—z)u(t—t') + dle elet) + 2u(t— r)b(¢—1')] ac de di — 
—2 [| acre t)u(t—ty) + b(t—rT)ult ty) | drdt + oon: = 


=|{ h(t)h(')[ (uw; T_T) ro)? Ds ts ee ee 


aa PD A(t) [r(w; t—t) + r(u, db; T—ty)] dt + r(u; 0) . 


_ o 


Le minimum de cette expression en fonction du choix de h (c’est-à-dire de 
l’opérateur du servomécanisme) s’obtient par le procédé classique du calcul 
des variations, en remplacant h par h + kAh, dot I + AI, et en écrivant 
que Von a 


AI + AI)] _ 
erga 


k=0 


Ces calculs sont sans difficulté, et conduisent è l’équation intégrale suivante 
dite de Wiener-Hopf, 


(co) 


(14) r(uyb ; tT — ty) —| h(t):-r(u; 7 —t)dt = 0 pourt> 0. 


aS 
Lorsque h vérifie cette équation, la valeur du minimum de J est 


[co] 


r(u; 0) — I h(t)[r(u; t—t,) + r(u, bit —b)] dr. 


_ 2 
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On remarquera que rien ne restreint le signe du décalage t,, qui peut, 
dans ce qui précède, étre une avance aussi bien qu’un retard; la connaissance 
des corrélations actuelles permet done, grace aux expressions obtenues, de 
prévoir la valeur probable du signal x(t) dans l’avenir; c'est Vamorce de la 
théorie de la prédiction de WIENER, que nous ne pouvons songer a 
développer ici; notons seulement que, du point de vue mathématique, cette 
partie est méme exempte de certaines difficultés qui grévent le cas du fil- 
trage (t) > 0). 

Nous n’insisterons pas ici sur la facon de résoudre l’équation de Wiener- 
Hopf; il suffira de savoir qu’on a recours, pour le faire, à la transformation 
de Fourier et 4 la décomposition des densités d’énergie spectrales en facteurs 
dont les pòles et les zéros sont convenablement répartis dans le plan des pul- 
sations complexes. 

Les considérations qui précèdent ont été développées pour un canal de 
téelécommunication (donc, & chaîne ouverte), mais sont directement appli- 
cables au cas des servomécanismes, conformément a notre équation (7), lorsque 
les signaux parasites arrivent par la méme voie que le message utile; elles 
ont done trait, sous cette forme, aux altérations dont le siòge est en amont 
du mécanisme. Une modification relativement légère permet de les étendre 
au cas d’un message altéré au cours de sa transformation, donc aux parasites 
qui naissent è l’intérieur du mécanisme, entre deux étages de celui-ci, par 
exemple. On aboutit alors 4 une équation intégrale de nature analogue a celle 
de (14), pour la discussion de laquelle nous renverrons & la thése de M. PE- 
LEGRIN. 

En résumé, la théorie du filtrage permet de dégager un message utile des 
perturbations qui peuvent l’accompagner, si l’un et les autres peuvent étre 
considérés comme stationnaires et ergodiques; il suffit alors de connaître pour 
le passé les fonctions de corrélation correspondantes, et le filtrage sera d’autant 
meilleur dans le présent que ces fonctions auront été mieux déterminées pour 
le passé. Toutefois, et c’est la une caractéristique précieuse de la méthode, 
méme une approximation grossière sur les fonctions de correlation permet 
d’améliorer notablement Vinformation transmise. 

Pour obtenir expérimentalement une fonction de corrélation, on peut en- 
registrer les messages et sommer les produits des ordonnées décalées de 7; 
graphiquement (fig. 5), on remplace l’intégrale par 


N-m 


Foo ae eT YnYn mo 
venti we 


A l'Institut de Mécanique des Fluides de Marseille, A. FAvRE enregistre sur 
rubans magnétiques, et fait les produits par difference des carrés de y(t) + 
+ y(t + t); l’appareil du Massachusetts Institute of Technology (LEE, CHEAT- 
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HAM et WIESNER) procède par «échantillonnage » au fur et à mesure de la 
réception du message: on fait la moyenne des produits d’ordonnées prises par 
couples a des intervalles 7 pour chaque couple, et 4 des intervalles 7 suffi- 
samment grands, entre couples 
successifs: il suffit d’émettre une 
suite de signaux pulsés conve- 
nablement écartés, commandant 
des créneaux dont l’amplitude est 
proportionnelle à Vun des élé- 
ments du couple, et la durée au 
second. 


i——m divisions ——~_———n divisions ———» 


‘ 
1 
' 
i 
' 
' 
' 
4 
4 
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' 
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i 
' 
' 
‘ 


Les courbes qui représentent 
r(y;t) en fonction de 7 peuvent 
Fig. 5. — (LEE). étre ainsi tracées. La considéra- 
tion de l’intégrale qui exprime 7 
permet de voir qu’une chaîne aléatoire (sans périodicités cachées) donnera 
une courbe présentant un maximum pour t = 0 et tendant vers 0 pour |r| 
croissant; par contre, la fonction 
d’auto-corrélation d’une fonction pé- @- signal (crenéaux de durée variable) 
riodique est elle-méme périodique. Il Foe fede) 
en résulte que dans le cas important 
où le message utile est périodique, le @-fonction d' autocorrélation du signal @ 
bruit parasite aléatoire, et où les 
deux sont incohérents, la séparation 
est aisée à partir de la courbe i i ia toe 
d’auto-corrélation du message total 
recu (fig. 6). 

Nous empruntons à la thèse de NANA 
J. P. Cosras un exemple de filtrage 
réalisé de la facon qui vient d’étre Os ii eee 
indiquée (fig. 7). Le signal « utile » A 
était ici une suite aléatoire de cré- . 
neaux, les parasites étant constitués EER AEE 
par un bruit de fond uniforme (bruit 
blanc); Pécart quadratique moyen 4 la sortie était de 0,517 (valeur calculée: 
0,4727, valeur pour un circuit parfait: 0,4657). 

Une propriété fondamentale des circuits linéaires, dont il a été fait usage 
dans le paragraphe 2, est la possibilité de séparer les variables lorsque le méca- 
nisme est à chaines multiples (ce qui correspond, du point de vue mathéma- 
tique, à réduire la matrice de transfert 4 ses axes). On peut prévoir que cette 
facilité ne subsistera pas si les signaux sont affectés de perturbations cohé- 
rentes (de la nature d’une diaphonie), qui doivent amener un couplage entre 


N divisions 


©-composante périodigue 
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les canaux: c'est ce que montre effectivement un calcul dù a Costas. Sup- 
posons qu’on ait affaire 4 n chaînes transportant des messages utiles w; et des 
bruits 6;, un calcul tout semblable a celui que nous avons exposé pour le cas 
unidimensionnel donne, pour l’écart quadratique moyen dans le canal i, l’ex- 
pression 


(co) Oo oO 


r(u;; 0 25 | hole orta + X| | Mhae — Marat, 


-— 0 - 0 - 


où Von a posé, pour alléger l’écriture, 


© AU 
©. sie = 


01s (t) = ult) + r(b;; #5 1) + (bi, By5 1). @- message regu v(4b 0) 


© - message aprés filtrage 


Qu(t) = TU, U;; Tt) + r(b;, U;; T) , 


©-message utile expédié 
Cet écart sera minimum lorsque les h vérifieront le 


systéme d’équations intégrales simultanées Fig. 7. — (Costas). 
ew 
Tdi (O | i 
> pro gar eu 0) [b=H1, 256, % 5 


—-o 


et la valeur du minimum sera 
r(u;; 0 3 [ht ii(T)o;(Tt—-t) dr. 


Ainsi, pour dégager la composante uw, du bruit qui l’altère, il faudra agir 
simultanément sur les opérateurs de tous les canaux. 


4. — Servomécanismes non linéaires. 


On connaît les difficultés que présente l’étude des systèmes régis par des 
équations différentielles non linéaires; elles sont particulicrement grandes 
s’agissant des servomécanismes, où les organes non linéaires qui interviennent 
le plus souvent sont de la nature des relais «par tout ou rien », done a très 
forte non-linéarité, non justiciables des méthodes de développement en série 
qui sont la ressource habituelle lorsque la non-linéarité est faible. On ne s’éton- 
nera donc pas si nous ne sommes pas encore en possession d’une théorie satis- 
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faisante des servomécanismes non linéaires; un certain nombre de résultats, 
fragmentaires mais déja utilisables, ont cependant pu étre obtenus récemment 
dans ce domaine, qui est en pleine évolution; ils procèdent, tout comme dans 
le cas linéaire, des deux points de vue fondamentaux et complémentaires du 
déterminisme et de la statistique. 

Nous nous bornerons ici à quelques indications relatives aux mécanismes 
à chaîne directe unique, tout en observant que le probleme pluridimensionnel 
ne peut plus se ramener a celui-ci par superposition, comme lorsque l’équation 
différentielle est linéaire. 

Les questions de stabilité sont extrémement difficiles 4 traiter d’une ma- 
nière vraiment correcte du point de vue mathématique; et la difficulté devient 
en général insurmontable lorsque l’ordre de l’équation est supérieur au second. 
Mais si Pon accepte de faire quelques approxi- 
mations et hypothèses qui ne paraissent pas dé- 

organe raisonnables a des physiciens, on peut arriver a 
vole Lngetce traiter, au moins dans ses grandes lignes, toute 

une catégorie, pratiquement importante, de pro- 
blèmes (DUTILA, KOCHENBURGER, LOEB). 

Supposons l’organe non linéaire inséré dans 


Fig. 8. 


une chaîne dont les autres éléments soient linéaires, et aient une fonction 
de transfert connue (fig. 8): nous aurons affaire è une équation différentielle 
du type 

Iu=—v, 


où L est un opérateur linéaire, et v = f(w) est la réponse fonctionnelle connue 
de Vorgane non linéaire. Cette réponse peut étre considérée comme résultat 
dune transformation du signal d’entrée w par un opérateur non linéaire o, 
de sorte que l’on aura, par exemple, au leu de (6), 


CÀ 
Pe; 


I 


(1+T+0)!x. 


La fagon la plus commode de traiter 7, on Va vu, est celle qui consiste a étudier 
la réponse des organes linéaires à un signal harmonique de pulsation w, d’où 
la fonction de transfert 7'(îw), puis de remplacer iw par la pulsation complexe p. 
Si les deux conditions suivantes se trouvent remplies, à savoir: 


1) la caractéristique f(w) de l’organe non linéaire est telle que dans la 
série de Fourier de f(cos wt) l’amplitude du fondamental est nettement plus 
élevée que celle des harmoniques de rangs supérieurs (condition réalisée, par 
exemple, par les relais, dont la réponse affecte la forme schématique de la 
fig. 9, et pour lesquels 

200 


f(u, cos wt) = > 


ITA) n 


1 — cos nz Wo 


Ò l 
sin (» ‘are COS 


cos nat ; 
Uy 
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2) ensemble des organes linéaires a une fonction de transfert telle que 
les harmoniques supérieurs de f(w) sont tous affaiblis de fagon notable, quelle 
que soit la pulsation fondamentale ©, alors la réponse de l’organe non linéaire 
au signal harmonique est pratiquement harmonique elle aussi, et réduite, dans 
Vexemple de la fig. 9, a 


4AM . Up 
vt — sin | are cos —] cos wt , 
TO) Uy 


de sorte que lon peut parler d’une « fonction de 
transfert », où interviennent la pulsation et Vampli- 
tude du signal d’entrée. Pour une amplitude don- 


P SCI seuil de 
nee, on pourra ecrire fonctionnement 


= X Big, 9 
1+T(p) + S(p)’ 


Z(p) 


et appliquer le critére de Hurwitz-Nyquist avec pour nouveau point critique 
celui d’affixe —(1-+ S(p)). Lorsqu’on fait varier u,, le point critique décrit 
une courbe du plan complexe, courbe qui comporte d’ailleurs un point d’arrét, 
correspondant au seuil wu): pour que le système soit stable, il faut et il suffit 
que cette courbe soit entièrement située & l’intérieur de celle qui représente 
la fonction de transfert de la partie linéaire. 

D’une facon plus générale, si l’on peut admettre l’existence d’une fonction 
de transfert T(x,, 0) de la chaîne, dépendant à la fois de la pulsation et de 
Vamplitude du signal d’entrée lorsque celui-ci est harmonique, l’affixe de 7° 
couvrira, lorsque les deux arguments prennent toutes les valeurs possibles, 
une région, généralement limitée, du plan complexe (ainsi, dans le cas pré- 
cédent, la réponse ne peut dépasser 40/7, et \ T| tend vers 0 lorque w, croit 
indéfiniment). Le système sera alors stable si le point — (1 + #0) est extérieur 
à cette région; au reste, dans le cas contraire, l’équation 


T(a,,@) =—1+%0 


donne Vamplitude et la pulsation auxquelles une auto-oscillation est possible. 
Cette généralisation du critere de Hurwitz-Nyquist est due a J. LOEB. 
Dans les considérations qui précédent, comme dans celles relatives aux 
systémes linéaires qui les ont inspirées, on voit que le rayon vecteur | T (ico) | 
joue un role important, encore qu’un peu vague: en particulier, pour un 
système dont la stabilité n’est pas simplement due au fait que l’affixe tourne 
autour du point critique un nombre de fois convenable dans les deux sens, 
mais & celui que le point critique est laissé a l’extérieur du lieu de l’affixe, il 
est clair que le système sera d’autant plus stable que la distance de l’affixe 
à Vorigine sera plus faible; de méme, sa réponse sera d’autant plus rapide 
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(au sens du n. 2) que cette distance tendra plus rapidement vers zéro au cours 
du temps. On peut partir de ce point de vue un peu imprécis, mais très général, 
pour étudier des systèmes non linéaires affranchis des hypothèses que néces- 
site le critére de LoEB; un mémoire récent de KANG et FETT donne, 4 cet 
égard, des résultats très intéressants, dont nous nous bornerons à indiquer 
sommairement le principe: 

On considère, pour un système régi par une équation d’ordre n, un espace 
de phase à n dimensions 


Ue — 08] da 


l’affixe du système est interprété comme un point matériel, dont le mouvement 
libre, régi par les terms linéaires, se fait suivant une famille déterminée de 
trajectoires. Les termes non linéaires du second membre jouent le ròle de 
forces appliquées, dont l’effet est de contraindre le point matériel a passer 
dune trajectoire sur l’autre è chaque instant. On peut alors introduire une 
métrique dans Vespace de phase en prenant pour fonction « distance » D(z, t) 
le lagrangien du système dynamique, ou mieux, compte tenu de la présence 
possible de termes oscillatoires, l’enveloppe de ce lagrangien, exprimé en co- 
ordonnées normales. Si D, est la fonction de distance en mouvement libre, 
et que F(t) soit ce que deviennent les termes non linéaires lorsqu’on y remplace 
les «; par leurs expressions en fonction de #, on trouve alors que 


dD dD, è 3 
ques nr + K(u) F(t) , 


où K est une forme linéaire homogène aux «,. Pour que D décroisse (cas d’un 
système stable), D, étant d’ailleurs décroissant a fortiori, il faut que les termes. 
non linéaires dépendent des uw; d’une facon convenable, et d’ailleurs différente 
suivant que l’affixe se trouve au-dessus ou au-dessous du plan A= 0; de 
plus, la valeur absolue de dD/dt donne une indication quantitative sur la 
rapidité du mécanisme, laquelle est naturellement différente pour chaque 
forme de signal d’entrée. 

Nous venons d’esquisser le tableau de état actuel en matière de stabilité 
des servomécanismes non linéaires. Si l’on examine maintenant le probleme 
du filtrage des messages statistiquement stationnaires, la présence d’organes 
non linéaires erée des difficultés d’une autre nature, qui tiennent surtout a 
la complication des calculs. Dans Vhypothése linéaire, la connaissance de la 
moyenne du second ordre permettait de définir la moyenne quadratique du 
signal de sortie, et, grace a l’intégrale de Duhamel, cette détermination était 
transferable vers l’amont: de sorte que finalement il suffisait de connaitre la 
covariance du signal entrée (supposé statistiquement stationnaire d’ordre 2) 
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pour connaitre le signal de sortie. Pour un mécanisme de transfert non linéaire, 
le principe de superposition, et par suite, l’intégrale de Duhamel, sont inap- 
plicables; on a, a leur place, une loi y = y(x), et si cette fonction est holo- 
morphe, soit 


(co) 


y(e) Do D A, ’ 


0 


il est clair que la détermination de y? exigera désormais la connaissance des 
moyennes de tous les ordres relatives au signal x(t). 

Heureusement, il est aisé d’obtenir toutes ces moyennes A la fois, si l’on 
connaît la fonction de répartition F(x(t)), ou la densité de probabilité p( x(t)), 
du signal: on sait, en effet, que la fonction caractéristique associée P,(w) est 
la transformée de Fourier de p(x), 


foo} 


E = I exp [tmx]: p(x) dx , 


—_ 


et que d’autre part le moment d’ordre n de x n’est autre chose que le coeffi- 
cient du développement de cette fonction 


Il est donc è prévoir que la considération de la fonction caractéristique associée 
permet de surmonter les difficultés dues a la non-linéarité du transmetteur. 

Effectivement, il suffit de supposer que la fonction y(x) qui définit Vorgane 
non linéaire possede une transformée de Fourier, pour que le signal de sortie 
se déduise en moyenne quadratique de la fonction caractéristique du signal 
d’entrée: car si i 

sd 
Y(iw)= | y(z) exp[— ina] de , 
Ò 

réciproquement 


n 


it 
ue) Y(iw) exp [iow] do, 


intégrale prise le long de l’axe réel, en contournant toutefois l’origine. Alors 


(oe) 


AR te) = y[ v(t) ] y[a(t + 1) ] at, 


-o 
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soit 


co 


r(Yy3 7) a di Il Y (iw) exp [iwa(t)| Y (io!) exp [tew'a(t + 7)] do do' = 


(ce) 


| Y (iw) ¥ (iw') P(w, 0‘; 7) dada’ , 


Je 


Ar? 


où P(0, o'; t) est, sous sa forme polaire, la fonetion caractéristique du signal 
composite w(t) + a(t + 7). k 

Si, comme au paragraphe 3, le signal d’entrée se compose d’une partie 
utile w(t) et d’un bruit parasite b(t), nous aurons 


1 rr 
r(y; t) = ia? I Y (iw) Y(iew') Py (@, @'3 7)P,(@, @'; t) dada’ . 


Les fonetions caractéristiques P, et P, peuvent théoriquement étre caleulées 
si Pon se donne les formes et les propriétés statistiques des deux signaux; il 
en est de méme des transformées des réponses Y(iw), lorsqu’on connaît la 
constitution de Porgane non linéaire: le probléme est done résolu en principe. 
De fait, Ricr et BLANG-LAPIERRE (qui ont établi, indépendamment Pun de 
l’autre, la théorie qui précède), et A leur suite MimoLetoN, ont effectué ces 
calculs, déjà assez pénibles dans certains cas simples intéressant la technique 
électronique: signaux harmoniques modulés en amplitude, bruit & moyenne 
nulle et è énergie moyenne donnée dans une bande de fréquence, transmetteur 
à détection y=|r| ou y= a"; le relevé expérimental des courbes d’auto- 
corrélation è la sortie a confirmé de fagon très satisfaisante ces prévisions 
théoriques (WEINBERG et KRAFT). On se trouve done en possession d’un outil 
dont le maniement peut s’avérer difficile, mais qui semble étre en principe 
adéquat pour traiter des problémes de filtrage non linéaire plus généraux. 


5. — Indications bibliographiques sommaires. 


Ad num. 1. On consultera sur la Cybernétique: 
N. WipnpRr: Cybernetics, or Control and Communication in the Animal and the 
Machine, « Actu. sei. ind. », n. 1053 (Paris, 1948). 
Sur la théorie de l’information, 
C. E. SHANNON: Bell System Tech. Journ., 27, 379-423 et 623-656 (1948), 


est le texte fondamental; on trouvera des compléments importants dans le 
recueil suivant: 
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Symposium on Information Theory, Report of Proceedings (London, Sept. 1950). 


Adnum.2: Les ouvrages sur la théorie générale des servomécanismes linéaires, 
et les descriptions techniques de différentes réalisations, sont innombrables. 
On pourra consulter 


H.M. James, N. B. Nicnors and R. S. Pamurps: Theory of servomechanisms , 
M.I.T. Radiation Lab. Series vol. 25 (New York, 1948). 


ou encore 


P. NasLin: Les Systémes asservis (Paris, 1951). 
Sur le critére du diagramme de phase, voir 


M. DEMONTVIGNIER et P. Lerbvre: Rev. gén. élect., 58, 263-279 (1949). 
P. Lerbvre: Mémorial arti. fr., 26 503-588 (1952). 
M. GoLoms et E. Uspin: Journ. Franklin Inst., 253, 29-57 (1952). 


Sur l’utilisation de l’écart quadratique moyen, 


P. R. AIGRAIN et E. M. WixLiams: Proc. I.R.E., 37, 873-879 (1949). 


On lira aussi avec profit les mémoires suivants, non utilisés dans le corps de 
l’exposé: 


J. KuNnTZMANN, J. DANIEL et M. Y. Ma: Rev. gén. d’élect., 61, 149-152, 532-538 
(1952). 


Ad num. 3: Les sources fondamentales sont 


N. WIENER: The extrapolation, interpolation and smoothing of stationary time series 
(New York, 1949). 

Y. W. Lew: Application of statistical Methods to commun ication problems, « Tech. 
Rep. 181», Res. Lab. of Electronics, M.I.T. (1950). 


Les travaux particuliers dont nous avons rendu compte sont les suivants: 


M. J. PÉLEGRIN: Bull. Soc. fr. élée. (7° 8.), 2, 595-600 (1952); These se. math. (Paris, 
1952). 

Y. W. Les, T. P. CaratHAM, J. B. WIESNER: Proc. I.R.E., 38, 1165-1171 (1950). 

A. Favre: Proc. VIIth Intern. Congress for appl. Mechanics, London, 2, 44-55 
(1948). 

J. P. Costas: Interference filtering, « Tech. Rep. 185 », Res. Lab. of Electronics, 
M.I.T. (1951). 


Ad num. 4: On consultera pour plus de détails les publications suivantes : 


J. R. DuriLn: Onde élec., 30, 438-445 (1950). 
R. J. KocHENBURGER: Blec. Eng., 69, 687-692 (1950). 
J. Lors: Ann. Télécomm., 6, 346-352 (1951). 
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(voir aussi, sur ce sujet: 
A. BLraquièRrE: Journ. de Phys., 13, 527-540 et 636-644 (1952)). 
Cc 


. L. Kane et G. H. Fert: Journ. appl. Phys., 24, 38-41 (1953). 

S. O. Rice: Bell. Syst. tech. Journ., 23, 282-332 (1944). 

A. BLANC-LAPIERRE: Compt. Rend. Ac. sc., 217, 73 (1943); 218, 924 (1944); Bull. 
Soc. fr. élect., (6° 8.), 5, 343-351 (1945). 

D. MippLETON: Qly. appl. maths., 5, 445-498 (1948). 

L. WEINBERG et L. G. Krarr: Experimental study of nonlinear devices by corre- 
lation methods, «Tech. Rep. 178», Res. Lab. of Electronics, M.I.T. (1951). 


On trouvera des résultats intéressants, qui n’ont pu trouver place dans notre 
exposé, dans 


F. H. Raymonp: Ann. Télécomm., 4, 234-240, 307,314 et 347-357 (1949). 
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